Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,725,840 Δ1 bp 5.0% intergenic (‑112/+85) cysT ← / ← cysK1 sulfate ABC transporter permease CysT/cysteine synthase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,725,8400C→. 5.0% 17.8 ‑1.9 ‑0.2159intergenic (‑112/+85)cysT/cysK1sulfate ABC transporter permease CysT/cysteine synthase
Reads supporting (aligned to +/- strand):  new base (3/5):  ref base (97/53):  total (101/58)
Fisher's exact test strand distribution p-value = 1.44e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.56e-01

Only 100 of 206 total aligned reads displayed.

TCCTTCCGTGAGCTTCGTCATTTCCAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGCGGGGCCTTTGAAAGACG.CGGGGAAACCGATCAGGCCGGGGCGCCG  >  NC_014034/2725740‑2725940
                                                                                                    |                                                                                                       
tccttccGTGAGCTTCGTCATTTCCAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTc                                                                                                         >  2:1627749/1‑101 (MQ=255)
 ccttccGTGAGCTTCGTCATTTCCAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTcc                                                                                                        <  2:983512/101‑1 (MQ=255)
 ccttccGTGAGCTTCGTCATTTCCAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTcc                                                                                                        <  1:4195122/101‑1 (MQ=255)
 ccttccGTGAGCTTCGTCATTTCCAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTcc                                                                                                        <  1:1338906/101‑1 (MQ=255)
  cttccGTGAGCTTCGTCATTTCCAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTccc                                                                                                       <  2:837191/101‑1 (MQ=255)
   ttccGTGAGCTTCGTCATTTCCAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTcccc                                                                                                      <  1:1566093/101‑1 (MQ=255)
      cGTGAGCTTCGTCATTTCCAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTccccccc                                                                                                   >  2:1970636/1‑101 (MQ=255)
      cGTGAGCTTCGTCATTTCCAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTccccccc                                                                                                   <  1:2261120/101‑1 (MQ=255)
      cGTGAGCTTCGTCATTTACAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTccccccc                                                                                                   >  2:498380/1‑101 (MQ=255)
       gTGAGCTTCGTCATTTCCAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTcccccccc                                                                                                  <  2:725669/101‑1 (MQ=255)
       gTGAGCTTCGTCATTTCCAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTcccccccc                                                                                                  <  2:435972/101‑1 (MQ=255)
       gTGAGCTTCGTCATTTCCAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTcccccccc                                                                                                  >  2:2993728/1‑101 (MQ=255)
       gTGAGCTTCGTCATTTCCAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTcccccccc                                                                                                  <  1:1296912/101‑1 (MQ=255)
         gAGCTTCGTCATTTCCAAAATATCCACTGTGTCAGTCGAGATTTAAAGGGATAAATCTCTCCATTCTTGGTCGTGAATATAGAATTTTATTcccccccccc                                                                                                <  2:3982217/101‑1 (MQ=255)
         gAGCTTCGTCATTTCCAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTcccccccccc                                                                                                <  1:109911/101‑1 (MQ=255)
          aGCTTCGTCATTTCCAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCa                                                                                               >  2:2531986/1‑101 (MQ=255)
          aGCTTCGTCATTTCCAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCa                                                                                               >  1:2981973/1‑101 (MQ=255)
           gCTTCGTCATTTCCAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCAc                                                                                              <  1:195104/101‑1 (MQ=255)
             ttCGTCATTTCCAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTc                                                                                            >  1:3410462/1‑101 (MQ=255)
             ttCGTCATTTCCAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATT.CCCCCCCCCACTcc                                                                                           >  1:3624644/1‑101 (MQ=255)
             ttCGTCATTTCCAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATT.CCCCCCCCCACTcc                                                                                           <  1:2832167/101‑1 (MQ=255)
             ctCGTCATTTCCAAAATATTCACTGTGTCAGACGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTAAAAATAGAATTTTATTCCCCCCCCCCACTc                                                                                            <  2:1506053/100‑1 (MQ=255)
              tCGTCATTTCCAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTcc                                                                                           <  2:2623208/101‑1 (MQ=255)
              tCGTCATTTCCAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTcc                                                                                           <  1:4021413/101‑1 (MQ=255)
                gTCATTTCCAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCCCTcccc                                                                                         >  1:2135486/1‑101 (MQ=255)
                 tCATTTCCAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTccccc                                                                                        <  1:1439080/101‑1 (MQ=255)
                  cATTTCCAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCg                                                                                       <  1:3070429/101‑1 (MQ=255)
                  cATTTCCAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCg                                                                                       >  2:2091182/1‑101 (MQ=255)
                  cATTTCCAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCg                                                                                       >  2:2651286/1‑101 (MQ=255)
                   aTTTCCAAAATATCCCCTGTGGCAGTCGAGATTGAAAGGGATAAATCCCCCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGt                                                                                      <  1:4216020/101‑1 (MQ=255)
                   aTTTCCAAAATATCCCCTGAGTCAGTAGAGATTGAAAGGGAAAAATCCCGCCATTCTTGGTAGTGAATATAAAATTTTATTCCCCCCCCCCACTCCCCCGt                                                                                      <  1:435893/101‑1 (MQ=255)
                    tttCCAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGtc                                                                                     >  2:3382103/1‑100 (MQ=255)
                       ccAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCg                                                                                  >  1:4418222/1‑101 (MQ=255)
                        cAAAATATCCCCTGTGGCAATCCAAATTGAAAGGGATAAATCCCGCCCTTCTTGGGCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCgg                                                                                 <  1:3488883/101‑1 (MQ=255)
                        cAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTCCCgg                                                                                 >  1:182538/1‑101 (MQ=255)
                         aaaaTATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCCCTCCCCCGTCCCggg                                                                                >  2:2146062/1‑101 (MQ=255)
                         aaaaTATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCggg                                                                                <  1:742057/101‑1 (MQ=255)
                             tatCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGgcgc                                                                            >  1:1862137/1‑101 (MQ=255)
                             tatCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGgcgc                                                                            >  1:3111734/1‑101 (MQ=255)
                             tatCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGgcgc                                                                            >  1:2740138/1‑101 (MQ=255)
                             tatCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGgcgc                                                                            >  2:602637/1‑101 (MQ=255)
                              atCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGcc                                                                           >  1:3503993/1‑101 (MQ=255)
                                 cACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAgc                                                                        >  1:944500/1‑100 (MQ=255)
                                 cACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGa                                                                        >  2:4177753/1‑101 (MQ=255)
                                  ccTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAc                                                                       <  2:1454969/100‑1 (MQ=255)
                                  aCTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAc                                                                       <  2:889631/101‑1 (MQ=255)
                                   cTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.g                                                                     <  1:1970765/101‑1 (MQ=255)
                                        tCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAg                                                                >  2:955682/1‑101 (MQ=255)
                                        tCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAg                                                                <  1:2630660/101‑1 (MQ=255)
                                           gttgaaaTTTAAAGGGATAACACTCTCCATTCCTGGTCGTGAATATAGAAATTTTTTCCCCCCCCCCCCGCCCCCGTGCCGGGGCGCCAGAC.GAAAAGCcc                                                             <  2:4293819/95‑1 (MQ=255)
                                             cGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTTCCCCGTCCCGGGGCGCCCGAC.GAAAAGGcccc                                                           >  2:1925507/1‑101 (MQ=255)
                                             cGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGcccc                                                           <  1:619372/101‑1 (MQ=255)
                                              gagaTTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCccg                                                          >  1:996866/1‑101 (MQ=255)
                                              gagaTTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCccg                                                          <  1:880014/101‑1 (MQ=255)
                                               agaTTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCC.CCCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCccg.cc                                                       <  1:3776662/101‑1 (MQ=255)
                                                gaTTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCccg.cc                                                       >  1:2151019/1‑101 (MQ=255)
                                                  ttGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCgg                                                     >  2:3211654/1‑101 (MQ=255)
                                                  ttGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCAGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCgg                                                     >  1:3836832/1‑101 (MQ=255)
                                                     aaaGGGATAAATCTCGCCGTTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCCGAC.GAAAAGGCCCCG.CCggggg                                                  >  1:2383714/1‑101 (MQ=255)
                                                     aaaGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCggggg                                                  >  1:2763097/1‑101 (MQ=255)
                                                     aaaGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCggggg                                                  >  2:1827216/1‑101 (MQ=255)
                                                     aaaGGGAAAAAACCCCCCATTCTTGGGCGTTAATATAAAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCggggg                                                  <  1:1019493/101‑1 (MQ=255)
                                                      aaGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCgggggg                                                 >  1:1207451/1‑101 (MQ=255)
                                                      aaGGGATAAATCTCGCCATTCTCGGTCGCGAATATAGAATTTTATTTCCCCCCCCCACCCCCCCGTGCCGGGGCGCCAGAC.GGAAGGGCCCCG.CCgcgggg                                                 >  2:1140098/1‑101 (MQ=255)
                                                        gggATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGggggcg                                               <  2:1731736/101‑1 (MQ=255)
                                                           aTAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCAACTCCCCCGT.CCCGGGCGCCCGAA.GGAAAGGCCCCGCCCGGggggcgggc                                            >  2:234287/1‑100 (MQ=255)
                                                             aaaTCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGCGGGGcc                                          >  1:3320505/1‑101 (MQ=255)
                                                             aaaTCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGGGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGCGGGGcc                                          >  2:1531417/1‑101 (MQ=255)
                                                               aTCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGCGGGGCCtt                                        >  1:424917/1‑101 (MQ=255)
                                                                   cGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGCGGGGCCTTTGaa                                    >  1:2775796/1‑101 (MQ=255)
                                                                    gCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGCGGGGCCTTTGaaa                                   >  2:3379028/1‑101 (MQ=255)
                                                                    gCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGCGGGGCCTTTGaaa                                   >  2:2998806/1‑101 (MQ=255)
                                                                         tCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGCGGGGCCTTTGAAAGAcg.c                             >  1:415981/1‑101 (MQ=255)
                                                                            tGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGCGGGGCCTTTGAAAGACG.Cggg                          >  2:1998330/1‑101 (MQ=255)
                                                                             gggcGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGCGGGGCCTTTGAAAGACG.Cgggg                         <  2:3806359/98‑1 (MQ=255)
                                                                             ggTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGCGGGTCCTTTGAAAGACG.Cgggg                         <  2:1207358/101‑1 (MQ=255)
                                                                                tgtaaaTATAGAATTTTATCCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGCGGGGCCTTTGAAAGACG.CGGGGaaa                      <  2:3836579/97‑1 (MQ=255)
                                                                                cGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGCGGGGCCTTTGAAAGACG.CGGGGaaa                      <  2:2868163/101‑1 (MQ=255)
                                                                                cGTGAATATAGAATTTTATTCCCCCCCCCC.CCCCCCCGTGCCCGGGCGCCCGACGGAAAAGGCCCCC.CCGGGGGGGGGGGCCTTGGAAA.ACG.CGGGGGAAc                     >  2:3521053/1‑101 (MQ=255)
                                                                                 gTGAATATAGAATTTTATCCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGCGGGGCCTTTGAAAGACG.CGGGGAAAc                     <  2:415941/101‑1 (MQ=255)
                                                                                  tGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGCGGGGCCTTTGAAAGACG.CGGGGAAAcc                    >  2:2794401/1‑101 (MQ=255)
                                                                                   gAATATGAAATTTTACCCCCCCCCCCCGCTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGCGGGGCCTTTGAAAGACG.CGGGGAAACCg                   <  1:2197313/101‑1 (MQ=255)
                                                                                   gAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGCGGGGCCTTTGAAAGACG.CGGGGAAACCg                   >  2:359046/1‑101 (MQ=255)
                                                                                    aaTATAGAATTTTATTCCCCCCCCCCCCCCCCCCG.GCCCGGGGGCCAGGA.GAAAAGGCCCCC.CCCGGGGGCGGGGCCTTTGAAAGACG.CGGGGGAAcccga                 >  2:1798570/1‑98 (MQ=255)
                                                                                    aaTATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGGGGGGCCTTTGAAAGACC.CGGGGAAACCGa                  >  2:2378210/1‑101 (MQ=255)
                                                                                    aaTATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGCGGGGCCTTTGAAAGACG.CGGGGAACCCGa                  >  1:3463073/1‑101 (MQ=255)
                                                                                         aGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGCGGGGCCTTTGAAAGACG.CGGGGAAACCGATCAgg             >  2:4003047/1‑101 (MQ=255)
                                                                                         aGAATTTTATT.CCCCCCCCCACCCCCCCGGGCCGGGGCGCCCGACGGAAAAGGCCCCC.CCGGGGGGGGGGGCCTTTG.AAGACG.CCGGGGAACCGAACCGGc            >  2:74978/1‑101 (MQ=255)
                                                                                          gAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGCGGGGCCTTTGAAAGACG.CGGGGAAACCGATCAGGc            >  2:2842135/1‑101 (MQ=255)
                                                                                          gAATTTTATGCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGCGGGGCCTTTGAAAGACG.CGGGGAAACCGATCAGGc            <  1:2930041/101‑1 (MQ=255)
                                                                                           aaTTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGCGGGGCCTTTGAAAGACG.CGGGGAAACCGATCAGGcc           <  2:3305209/101‑1 (MQ=255)
                                                                                            aTTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGCGGGGCCTTTGAAAGACG.CGGGGAAACCGATCAGGCCg          >  1:114142/1‑101 (MQ=255)
                                                                                            aTTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGCGGGGCCTTTGAAAGACG.CGGGGAAACCGAGCAGGCCg          >  2:1275280/1‑101 (MQ=255)
                                                                                             ttttATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCC.CCGGGGGGCGGGGCCTTTGAAAGACC.CGGGGAAACCCATCAGGCCgg         >  1:3140328/1‑101 (MQ=255)
                                                                                               ttATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGCGGGGCCTTTGAAAGACG.CGGGGAAACCGATCAGGCCgggg       >  2:15708/1‑101 (MQ=255)
                                                                                                tATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGCGGGGCCTTTGAAAGACG.CGGGGAAACCGATCAGGCCGGGgc      <  2:3111510/101‑1 (MQ=255)
                                                                                                 aTTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCCGAC.GAAAAGGCCCCC.CCGGGGGGCGGGGCCTTTGAAAGACC.GGGGGCAACCGGTCGGGCCGGGgca     >  1:3866006/1‑100 (MQ=255)
                                                                                                   tCCCCCCCCCCCCCCCCCCTGGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGGGGGGCCTTTGAAA.ACGCCGGGGAAACCGACCGGGCCGGGGCCcc   >  2:325087/1‑101 (MQ=255)
                                                                                                    ccccccccccCCTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGCGGGGCCTTTGAAAGACG.CGGGGAAACCGATCAGGCCGGGGCGCCg  <  2:3607800/101‑1 (MQ=255)
                                                                                                    ccccccccccACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGCGGGGCCTTTGAAAGACG.CGGGGGAACCGATCAGGCCGGGGCGCCg  >  1:3214721/1‑101 (MQ=255)
                                                                                                    |                                                                                                       
TCCTTCCGTGAGCTTCGTCATTTCCAAAATATCCACTGTGTCAGTCGAGATTGAAAGGGATAAATCTCGCCATTCTTGGTCGTGAATATAGAATTTTATTCCCCCCCCCCACTCCCCCGTGCCGGGGCGCCAGAC.GAAAAGGCCCCG.CCGGGGGGCGGGGCCTTTGAAAGACG.CGGGGAAACCGATCAGGCCGGGGCGCCG  >  NC_014034/2725740‑2725940

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