Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 985,690 Δ1 bp 100% intergenic (+29/‑55) betB → / → betA betaine‑aldehyde dehydrogenase/choline dehydrogenase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034985,6900C→.100.0% 54.2 232intergenic (+29/‑55)betB/betAbetaine‑aldehyde dehydrogenase/choline dehydrogenase
Reads supporting (aligned to +/- strand):  new base (82/129):  ref base (18/0):  total (103/129)
Fisher's exact test strand distribution p-value = 1.30e-07
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.38e-15

Only 100 of 285 total aligned reads displayed.

GCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGCCGGTCCCCCCCCGGACCCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCGGGCG  >  NC_014034/985590‑985784
                                                                                                    |                                                                                              
gCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGAc                                                                                                <  2:5838247/101‑1 (MQ=255)
 ccGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGAcc                                                                                               <  1:842555/101‑1 (MQ=255)
   gcgATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGGCCCCCCCCCGGcccc                                                                                             >  2:2375680/1‑101 (MQ=255)
   gcgATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGGCCCCCCCCCGCcccc                                                                                             >  1:3295505/1‑101 (MQ=255)
    cgATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGAccccc                                                                                            >  2:2955046/1‑101 (MQ=255)
    cgATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGAccccc                                                                                            <  1:3150177/101‑1 (MQ=255)
    cgATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGGCCCCCCCCGGAccccc                                                                                            >  2:6699541/1‑101 (MQ=255)
    cgATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGGCCCCCCCCGGAccccc                                                                                            >  2:6591203/1‑101 (MQ=255)
     gATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCTGGCCGGGGGGTCCCCCCCCGGAcccccc                                                                                           <  1:5067986/101‑1 (MQ=255)
     gATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGCcccccc                                                                                           >  2:4813542/1‑101 (MQ=255)
     gATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGAcccccc                                                                                           <  1:2435087/101‑1 (MQ=255)
      aTCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGAccccccg                                                                                          <  2:3490675/101‑2 (MQ=255)
         gAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGAccccccggcg                                                                                       <  2:1897641/101‑5 (MQ=255)
         gAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGGCCCCCCCCGGAccccccggcg                                                                                       >  1:6986476/1‑97 (MQ=255)
            cACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGACGGGGGGGCCCCCCCCGGAccccccggcggat                                                                                    >  2:5254035/1‑94 (MQ=255)
             aCTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAACCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTAtt                                                                                  <  2:4611187/101‑1 (MQ=255)
              cTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCCGACCCCCCCgcgtgttt                                                                                  >  2:6569902/1‑94 (MQ=255)
               tataCGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTg                                                                                <  2:6752582/101‑1 (MQ=255)
               tataCGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTg                                                                                <  1:1710696/101‑1 (MQ=255)
                ataCGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGGCCCCCCCCGGACCCCCCC.GCGGATTTgc                                                                               >  1:6487747/1‑101 (MQ=255)
                 taCGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGGCCCCCCCCGGACCCCCCCgcggactttcg                                                                               >  2:610294/1‑91 (MQ=255)
                    gCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCa                                                                           <  1:19902/101‑1 (MQ=255)
                    gCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCa                                                                           <  2:5534647/101‑1 (MQ=255)
                    gCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGGCCCCCCCCGG.CCCCCCCGGCGTATTTGCGCCa                                                                           >  2:2236005/1‑101 (MQ=255)
                     cAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCaa                                                                          <  2:4275310/101‑1 (MQ=255)
                     cAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGGCCCCCCCCCGACCCCCCC.GCGGACTTGCGCCaa                                                                          >  1:4794300/1‑101 (MQ=255)
                      aGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAg                                                                         <  2:5391381/101‑1 (MQ=255)
                        gTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGaa                                                                       <  2:4117207/101‑1 (MQ=255)
                         tGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATGTGCGCCAAGaaa                                                                      <  1:6401733/101‑1 (MQ=255)
                         tGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGGCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCCAGaaa                                                                      >  1:2817806/1‑101 (MQ=38)
                         tGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGGCCCCCCCCCGGCCCCCCCGCCGATTTTGGCCCAAAaa                                                                       >  1:2762313/1‑101 (MQ=255)
                          gAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGaaaa                                                                     >  2:6249583/1‑101 (MQ=255)
                           aaGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGGCCCCCCC.GCGTATTTGCGCCAAGaaaaa                                                                    >  2:89207/1‑101 (MQ=255)
                           aaGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGaaaaa                                                                    <  1:119861/101‑1 (MQ=255)
                            aGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGaaaaaa                                                                   <  1:3788305/101‑1 (MQ=255)
                            aGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGaaaaaa                                                                   <  2:4444981/101‑1 (MQ=255)
                            aGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGaaaaaa                                                                   <  1:610314/101‑1 (MQ=255)
                             ggCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAc                                                                  <  1:2846215/101‑1 (MQ=255)
                             ggCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGGCTTTGGGGCAAGaaaaaaa                                                                  >  2:360942/1‑100 (MQ=255)
                             ggCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGGCCCCCCCCGGACCCCCCC.GCGTAATTGCGCCCAGAAAAAAc                                                                  >  2:4204117/1‑101 (MQ=255)
                                 gTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGGCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAggg                                                              >  1:6896299/1‑101 (MQ=38)
                                  tCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGt                                                             <  2:2601291/101‑1 (MQ=255)
                                   cTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGGATTTGCGCCAAGAAAAAAAAggggt                                                            >  1:4486468/1‑99 (MQ=255)
                                    tATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTc                                                           <  1:4167946/101‑1 (MQ=255)
                                    tATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGGTc                                                           >  2:6926397/1‑101 (MQ=38)
                                    tATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGGCCCCCCCCCGACCCCCCC.GCGTATTTGCGCCAAGAAAAAACAGGGTTc                                                           >  2:2875437/1‑101 (MQ=38)
                                     aTGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCa                                                          <  2:376965/101‑1 (MQ=255)
                                      tGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAg                                                         <  1:6725882/101‑1 (MQ=255)
                                        tCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGGCCCCCCCCCGACCCCCCCGG.GGATTTGCGCCCAGAAAAAACAGGGGTCAAGc                                                       >  2:2421102/1‑101 (MQ=255)
                                         cGGCATGGGGCCGGTCGAGACCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCa                                                      <  2:3950862/101‑1 (MQ=38)
                                          ggCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCag                                                     <  2:6487417/101‑1 (MQ=255)
                                          ggCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCag                                                     >  1:1217351/1‑101 (MQ=255)
                                          ggCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGG.GCCCCCCGGCGTATTTGCGCCAAGCAAAAACAGGGTTCCGGCag                                                     >  2:5116959/1‑101 (MQ=255)
                                            cATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGGCCCCCCCCCGA.CCCCCCGGCGGATTTGCGCCCAGAAAAAACAGGGGTCAGGCggag                                                   >  1:1440156/1‑101 (MQ=255)
                                             aTGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGACCCCCCCCGG.GCCCCCCGGGGACTTTTCCCCCCGAAAAAACAGGGGCCGAGCAGAgg                                                  >  1:7129506/1‑101 (MQ=255)
                                                 ggCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTcg                                              >  1:4191650/1‑101 (MQ=38)
                                                 ggCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTcg                                              <  1:4915039/101‑1 (MQ=255)
                                                     ggTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTcgcgcg                                          <  1:1138851/101‑1 (MQ=255)
                                                     ggTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTcgcgcg                                          >  1:1016013/1‑101 (MQ=37)
                                                      gTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGGCCCCCCCCCGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGGTCAGGCAGAGGGGcgcgcgg                                         >  1:1548107/1‑100 (MQ=37)
                                                       tCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCgat                                        >  2:5443943/1‑101 (MQ=38)
                                                         gagaGCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGG.GACCCCCGGCGGATTTGCGCCAAGAAAAAACAGGGGCCAGGGAGAGGGGCGCGCgatgg                                      >  2:1375840/1‑101 (MQ=25)
                                                          agagCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGCTCAGGCAGAGGGTCGCGCgatgga                                     >  2:7128078/1‑101 (MQ=38)
                                                          agagCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGGATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCgatgga                                     >  2:2075771/1‑101 (MQ=38)
                                                           gagCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCgatggat                                    <  1:5445871/101‑1 (MQ=255)
                                                           gagCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCgatggat                                    <  2:4375908/101‑1 (MQ=255)
                                                            agCCCCTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGGCGAGGGTcgccggcaggatg                                   >  2:6950386/1‑91 (MQ=25)
                                                             gCCCCTATTGAAACGCAGGCCGGGGGGGCCCCCCCCCGACCCCCCC.GCGGATTTGCGCCAAGAAAAAACAAGGGTCCGGCAGAGGGGCGCGCgatgggtgg                                  >  1:107758/1‑100 (MQ=255)
                                                              ccccTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCg                                 <  2:5320024/101‑1 (MQ=255)
                                                                ccTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGa                               <  2:1432639/101‑1 (MQ=255)
                                                                ccTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGa                               <  2:1452390/101‑1 (MQ=255)
                                                                ccTATTGAAACGCAGGCCGGGGGGGCCCCCCCCCGACCCCCCC.GCGGATTTTGGCCCAGAAAAAACAGGGGTCAGGCAGAGGGTGGCGCGATGGATGCGGa                               >  1:4999950/1‑101 (MQ=255)
                                                                 cTATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGGTCAGGCAGAGGGTCGCGCGATGGATGCGGAt                              >  1:6556326/1‑101 (MQ=37)
                                                                  tATTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGAtt                             <  2:2762211/101‑1 (MQ=255)
                                                                   aTTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTa                            <  1:2578766/101‑1 (MQ=255)
                                                                   aTTGAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGGCGCGCGATGGATGCGGATTa                            >  1:6356043/1‑101 (MQ=37)
                                                                      gAAACGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTAtgt                         >  2:1745270/1‑101 (MQ=38)
                                                                       aaaCGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTAtgtg                        <  1:1906971/101‑1 (MQ=255)
                                                                       aaaCGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTAtgtg                        <  1:7050515/101‑1 (MQ=255)
                                                                        aaCGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCCAGAAAAAACCGGGTTCCGGCCGAGGGTCGCGCGATGGATGCGGATTAtgtgg                       >  2:1014344/1‑100 (MQ=37)
                                                                         aCGCAGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGAt                      <  2:3236784/101‑1 (MQ=255)
                                                                          cGCAGGCCGGGGGGTCCCGCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATc                     <  2:1267213/101‑1 (MQ=38)
                                                                           gCAGGCCGGCGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCg                    <  2:3295389/101‑1 (MQ=255)
                                                                             cggccggggGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATAGATGCGGATTATGTGATCGtt                  <  1:4281027/100‑1 (MQ=38)
                                                                             aGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGtt                  >  2:1588168/1‑101 (MQ=255)
                                                                             aGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGtt                  <  1:2237263/101‑1 (MQ=255)
                                                                             aGGCCGGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGtt                  >  1:2703061/1‑101 (MQ=255)
                                                                              ggccggggGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTg                 >  2:5010506/1‑101 (MQ=255)
                                                                              ggccggggGGGCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTg                 <  1:209389/101‑1 (MQ=38)
                                                                                 cggggggTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGGTCAGGCAGAGGGTCGCGCGATGGATGCGGGTTATGTGATCGTTgggg              >  2:6655886/6‑101 (MQ=38)
                                                                                  ggggggTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGc             >  2:1034136/5‑101 (MQ=255)
                                                                                  ggggggTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGc             <  1:5291749/97‑1 (MQ=255)
                                                                                  gggaggTCCCCCCCCGGA.CCCCCCGGCGTACTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCAGATTATGTGATCGTTGGGGc             >  1:5618037/5‑101 (MQ=38)
                                                                                    ggggTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAg           <  1:1305204/99‑1 (MQ=255)
                                                                                       gTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGcc        >  2:572969/1‑100 (MQ=255)
                                                                                       gTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCt        <  1:6208070/101‑1 (MQ=255)
                                                                                         ccccccccGGA.CCCCCCGGCGTATTTTCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGCTGCGGATTATGTGATCGTTGGGGCAGGCTCg      >  1:395258/1‑101 (MQ=255)
                                                                                            cccccGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCGggc   <  1:4013274/101‑1 (MQ=255)
                                                                                            cccccGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCGggc   <  1:1048417/101‑1 (MQ=255)
                                                                                             ccccGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCGggcg  >  1:3740633/1‑101 (MQ=255)
                                                                                                    |                                                                                              
GCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGCCGGTCCCCCCCCGGACCCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCGGGCG  >  NC_014034/985590‑985784

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG