Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F224 I0 R1 1457 284.0 9290224 95.7% 8890744 150.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,488,966 G→C 100% A533P (GCC→CCC)  nirB → nitrite reductase, large subunit, NAD(P)H‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,488,9660GC97.1% 1077.4 / 19.4 273A533P (GCC→CCC) nirBnitrite reductase, large subunit, NAD(P)H‑binding
Reads supporting (aligned to +/- strand):  ref base G (3/5);  new base C (148/117);  total (151/122)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 278 total aligned reads displayed.

AACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTG  >  CP009273/3488819‑3489112
                                                                                                                                                   |                                                                                                                                                  
aacCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGccc                                                                                                                                                  <  1:2433667/150‑1 (MQ=255)
 acCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCcct                                                                                                                                                 <  2:679413/150‑1 (MQ=255)
       gCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCctgct                                                                                                                                           <  1:4148994/150‑1 (MQ=255)
         gAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTgg                                                                                                                                         >  2:3860585/1‑150 (MQ=255)
         gAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTgg                                                                                                                                         >  1:2813455/1‑150 (MQ=255)
         gAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTgg                                                                                                                                         >  1:232083/1‑150 (MQ=255)
         gAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTgg                                                                                                                                         >  1:1802697/1‑150 (MQ=255)
         gAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTgg                                                                                                                                         >  1:4098450/1‑150 (MQ=255)
         gAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTgg                                                                                                                                         >  1:1966740/1‑150 (MQ=255)
          aaCACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGa                                                                                                                                        >  1:3630844/1‑150 (MQ=255)
          aaCACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGa                                                                                                                                        >  2:4469591/1‑150 (MQ=255)
           acacTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGaa                                                                                                                                       <  1:402972/150‑1 (MQ=255)
           acacTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGATTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGaa                                                                                                                                       >  1:3579209/1‑150 (MQ=255)
            cacTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAAc                                                                                                                                      <  2:477866/150‑1 (MQ=255)
              cTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGa                                                                                                                                    >  2:20960/1‑150 (MQ=255)
               tttGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGaa                                                                                                                                   >  1:3684262/1‑150 (MQ=255)
                  gCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATAc                                                                                                                                <  1:1612013/150‑1 (MQ=255)
                   cTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAGGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACa                                                                                                                               <  1:2514851/150‑1 (MQ=255)
                   cTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACa                                                                                                                               <  2:713246/150‑1 (MQ=255)
                   cTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACa                                                                                                                               <  1:3332563/150‑1 (MQ=255)
                    ttATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACAt                                                                                                                              >  1:1506632/1‑150 (MQ=255)
                       ttCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCt                                                                                                                           <  1:833990/150‑1 (MQ=255)
                         cgcgTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGa                                                                                                                         <  1:919713/150‑1 (MQ=255)
                           cgTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAg                                                                                                                       <  2:4547386/150‑1 (MQ=255)
                           cgTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAg                                                                                                                       <  1:1072901/150‑1 (MQ=255)
                            gTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGc                                                                                                                      >  2:342167/1‑150 (MQ=255)
                            gTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGc                                                                                                                      <  1:3208146/150‑1 (MQ=255)
                                 gAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGaa                                                                                                                 >  2:1217043/1‑150 (MQ=255)
                                  aaCTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAAc                                                                                                                <  1:3172348/150‑1 (MQ=255)
                                  aaCTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAAc                                                                                                                <  2:826861/150‑1 (MQ=255)
                                      gTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATAc                                                                                                            <  1:837533/150‑1 (MQ=255)
                                      gTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATAc                                                                                                            >  1:447020/1‑150 (MQ=255)
                                       ttCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACt                                                                                                           <  1:1729099/150‑1 (MQ=255)
                                       ttCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACt                                                                                                           >  1:1203468/1‑150 (MQ=255)
                                        tCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTc                                                                                                          <  1:523654/150‑1 (MQ=255)
                                           aTTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGc                                                                                                       >  1:4449840/1‑150 (MQ=255)
                                            tttGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCt                                                                                                      <  1:2066892/150‑1 (MQ=255)
                                                      gTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCt                                                                                            >  1:3098682/1‑150 (MQ=255)
                                                         gAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAAc                                                                                         >  1:726614/1‑150 (MQ=255)
                                                           aGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGa                                                                                       >  1:86180/1‑150 (MQ=255)
                                                            ggCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGAc                                                                                      >  2:2648039/1‑150 (MQ=255)
                                                              cATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACaa                                                                                    >  1:1863662/1‑150 (MQ=255)
                                                              cATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACaa                                                                                    >  2:312015/1‑150 (MQ=255)
                                                                 tAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACtt                                                                                 >  2:1666229/1‑150 (MQ=255)
                                                                 tAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACtt                                                                                 <  2:4405814/150‑1 (MQ=255)
                                                                  aaaaCCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTc                                                                                >  2:2127927/1‑150 (MQ=255)
                                                                    aaCCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCt                                                                              <  2:3675124/150‑1 (MQ=255)
                                                                     aCCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTc                                                                             >  2:1364810/1‑150 (MQ=255)
                                                                       cTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGc                                                                           >  1:3747224/1‑150 (MQ=255)
                                                                        ttCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCt                                                                          <  2:1823363/150‑1 (MQ=255)
                                                                          cGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTaa                                                                        >  2:2278307/1‑150 (MQ=255)
                                                                            aagaaCTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACa                                                                      <  1:4365764/150‑1 (MQ=255)
                                                                             agaaCTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACAt                                                                     >  2:731549/1‑150 (MQ=255)
                                                                                    ctgGCGAAACACGGCATAGGCTACGGTTGTGTAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGaaa                                                              <  2:2435473/150‑1 (MQ=255)
                                                                                    ctgGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGaaa                                                              <  2:1752238/150‑1 (MQ=255)
                                                                                         gAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACgg                                                         >  2:1578391/1‑150 (MQ=255)
                                                                                          aaaCACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGc                                                        >  2:2685618/1‑150 (MQ=255)
                                                                                          aaaCACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGc                                                        >  1:4061513/1‑150 (MQ=255)
                                                                                          aaaCACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGc                                                        >  1:2176407/1‑150 (MQ=255)
                                                                                           aaCACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCa                                                       >  2:168466/1‑150 (MQ=255)
                                                                                            acacGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCAc                                                      <  2:4126391/150‑1 (MQ=255)
                                                                                                ggCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTAc                                                  >  2:1946021/1‑150 (MQ=255)
                                                                                                 gCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACt                                                 <  1:521543/150‑1 (MQ=255)
                                                                                                 gCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACt                                                 >  2:2917846/1‑150 (MQ=255)
                                                                                                   aaaGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCg                                               >  2:3942265/1‑150 (MQ=255)
                                                                                                      ggCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTg                                            <  2:2822220/150‑1 (MQ=255)
                                                                                                          aCGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATcc                                        >  2:3753216/1‑150 (MQ=255)
                                                                                                               tgtgAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGt                                   >  1:2135776/1‑150 (MQ=255)
                                                                                                               tgtgAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGt                                   >  2:1435729/1‑150 (MQ=255)
                                                                                                               tgtgAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGt                                   >  1:3659618/1‑150 (MQ=255)
                                                                                                                gtgAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGtt                                  <  1:2218873/150‑1 (MQ=255)
                                                                                                                gtgAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGtt                                  <  1:3864080/150‑1 (MQ=255)
                                                                                                                 tgAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTtc                                 <  1:1476109/150‑1 (MQ=255)
                                                                                                                 tgAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTtc                                 >  1:101173/1‑150 (MQ=255)
                                                                                                                  gAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTtct                                <  2:748187/150‑1 (MQ=255)
                                                                                                                    aGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATCCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTcc                              >  2:4050618/1‑150 (MQ=255)
                                                                                                                     gTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCg                             <  2:4529943/150‑1 (MQ=255)
                                                                                                                       ttGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCggg                           <  2:18938/150‑1 (MQ=255)
                                                                                                                          tAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCgg                        >  1:1384141/1‑150 (MQ=255)
                                                                                                                          tAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCgg                        >  1:3059146/1‑150 (MQ=255)
                                                                                                                           aaaCCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGt                       >  1:1824304/1‑150 (MQ=255)
                                                                                                                             accaaccGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGa                     <  2:1090881/150‑1 (MQ=255)
                                                                                                                             accaaccGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGa                     <  2:859259/150‑1 (MQ=255)
                                                                                                                             accaaccGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGa                     >  2:3034191/1‑150 (MQ=255)
                                                                                                                             accaaccGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACGTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGa                     <  2:4268086/150‑1 (MQ=255)
                                                                                                                                 accGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATc                 <  2:1431907/150‑1 (MQ=255)
                                                                                                                                 accGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATc                 <  2:147448/150‑1 (MQ=255)
                                                                                                                                  ccGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCa                >  2:1284641/1‑150 (MQ=255)
                                                                                                                                   cGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCAc               <  2:1966740/150‑1 (MQ=255)
                                                                                                                                    gTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCAcc              >  2:4528594/1‑150 (MQ=255)
                                                                                                                                      cGGTTCGCTGCTGCCCTCCTGCTGGAACGACTACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCAcccc            >  2:2690137/1‑150 (MQ=255)
                                                                                                                                      cGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCAcccc            <  2:3359866/150‑1 (MQ=255)
                                                                                                                                      cGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCAcccc            <  2:3630844/150‑1 (MQ=255)
                                                                                                                                        gTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCgg          >  2:261273/1‑150 (MQ=255)
                                                                                                                                            gctgctGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAgg      <  1:1257343/150‑1 (MQ=255)
                                                                                                                                            gctgctGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAgg      <  2:3592084/150‑1 (MQ=255)
                                                                                                                                              tgctgcCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGc    <  2:3101088/150‑1 (MQ=255)
                                                                                                                                              tgctgcCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGc    >  1:4224619/1‑150 (MQ=255)
                                                                                                                                               gctgCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCt   >  2:1596207/1‑150 (MQ=255)
                                                                                                                                                ctgCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTg  >  2:898033/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTG  >  CP009273/3488819‑3489112

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTGATG  >  CP009273/3488819‑3489115
                                                                                                                                                   |                                                                                                                                                     
AACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTTTTCCATTTGATCCGCGTTGAATGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGATGTGAAGTTTGTAAACCAACCGTCTGTTCGCTGCTGCCC                                                                                                                                                     >  SRR19973072.1927758/1‑150 (MQ=60)
AACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCC                                                                                                                                                     <  SRR19973072.2438939/150‑1 (MQ=60)
AACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCC                                                                                                                                                     >  SRR19973072.1927760/1‑150 (MQ=60)
 ACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCT                                                                                                                                                    >  SRR19973072.4128586/1‑150 (MQ=60)
 ACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCT                                                                                                                                                    >  SRR19973072.4128571/1‑150 (MQ=60)
    TGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCT                                                                                                                                                 <  SRR19973072.682467/150‑1 (MQ=60)
    TGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCT                                                                                                                                                 <  SRR19973072.4046062/150‑1 (MQ=60)
    TGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCT                                                                                                                                                 <  SRR19973072.4045685/150‑1 (MQ=60)
     GTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTG                                                                                                                                                >  SRR19973072.836324/1‑150 (MQ=60)
     GTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTG                                                                                                                                                >  SRR19973072.836331/1‑150 (MQ=60)
       GCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCT                                                                                                                                              <  SRR19973072.4157691/150‑1 (MQ=60)
         GAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGG                                                                                                                                            >  SRR19973072.1806635/1‑150 (MQ=60)
         GAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGG                                                                                                                                            >  SRR19973072.4107057/1‑150 (MQ=60)
         GAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGG                                                                                                                                            >  SRR19973072.2819481/1‑150 (MQ=60)
         GAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGG                                                                                                                                            >  SRR19973072.1970950/1‑150 (MQ=60)
         GAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGG                                                                                                                                            >  SRR19973072.232657/1‑150 (MQ=60)
          AACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGA                                                                                                                                           >  SRR19973072.3638700/1‑150 (MQ=60)
          AACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGA                                                                                                                                           >  SRR19973072.3638495/1‑150 (MQ=60)
           ACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAA                                                                                                                                          <  SRR19973072.403862/150‑1 (MQ=60)
           ACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGATTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAA                                                                                                                                          >  SRR19973072.3586747/1‑150 (MQ=60)
            CACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAAC                                                                                                                                         >  SRR19973072.72809/1‑150 (MQ=60)
              CTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGA                                                                                                                                       <  SRR19973072.2623792/150‑1 (MQ=60)
               TTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAA                                                                                                                                      >  SRR19973072.1747056/1‑150 (MQ=60)
               TTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAA                                                                                                                                      >  SRR19973072.3691952/1‑150 (MQ=60)
               TTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAA                                                                                                                                      >  SRR19973072.3692029/1‑150 (MQ=60)
                 TGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATA                                                                                                                                    >  SRR19973072.1756495/1‑150 (MQ=60)
                  GCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATAC                                                                                                                                   >  SRR19973072.3599641/1‑150 (MQ=60)
                  GCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATAC                                                                                                                                   <  SRR19973072.1615508/150‑1 (MQ=60)
                  GCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATAC                                                                                                                                   <  SRR19973072.1615505/150‑1 (MQ=60)
                  GCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATAC                                                                                                                                   <  SRR19973072.1615222/150‑1 (MQ=60)
                  GCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATAC                                                                                                                                   >  SRR19973072.3599663/1‑150 (MQ=60)
                   CTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAGGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACA                                                                                                                                  <  SRR19973072.2520314/150‑1 (MQ=60)
                   CTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACA                                                                                                                                  <  SRR19973072.786134/150‑1 (MQ=60)
                   CTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACA                                                                                                                                  <  SRR19973072.3339634/150‑1 (MQ=60)
                   CTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACA                                                                                                                                  <  SRR19973072.2011892/150‑1 (MQ=60)
                    TTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACAT                                                                                                                                 <  SRR19973072.953257/150‑1 (MQ=60)
                    TTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACAT                                                                                                                                 >  SRR19973072.1509840/1‑150 (MQ=60)
                      ATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTC                                                                                                                               >  SRR19973072.3107625/1‑150 (MQ=60)
                       TTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCT                                                                                                                              >  SRR19973072.922938/1‑150 (MQ=60)
                       TTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCT                                                                                                                              <  SRR19973072.835781/150‑1 (MQ=60)
                         CGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGA                                                                                                                            <  SRR19973072.921692/150‑1 (MQ=60)
                         CGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGA                                                                                                                            <  SRR19973072.921711/150‑1 (MQ=60)
                           CGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAG                                                                                                                          <  SRR19973072.1075188/150‑1 (MQ=60)
                            GTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGC                                                                                                                         <  SRR19973072.3214930/150‑1 (MQ=60)
                             TCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCC                                                                                                                        <  SRR19973072.492456/150‑1 (MQ=60)
                             TCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCC                                                                                                                        <  SRR19973072.492812/150‑1 (MQ=60)
                             TCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCC                                                                                                                        <  SRR19973072.3371881/150‑1 (MQ=60)
                              CAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCG                                                                                                                       >  SRR19973072.2216603/1‑150 (MQ=60)
                              CAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCG                                                                                                                       <  SRR19973072.1367479/150‑1 (MQ=60)
                               AGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGG                                                                                                                      >  SRR19973072.1819001/1‑150 (MQ=60)
                               AGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGG                                                                                                                      <  SRR19973072.2132542/150‑1 (MQ=60)
                                  AACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAAC                                                                                                                   >  SRR19973072.4224562/1‑150 (MQ=60)
                                  AACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAAC                                                                                                                   <  SRR19973072.3179041/150‑1 (MQ=60)
                                  AACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAAC                                                                                                                   >  SRR19973072.2145407/1‑150 (MQ=60)
                                  AACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAAC                                                                                                                   <  SRR19973072.167938/150‑1 (MQ=60)
                                  AACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAAC                                                                                                                   <  SRR19973072.167879/150‑1 (MQ=60)
                                  AACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAAC                                                                                                                   >  SRR19973072.2145316/1‑150 (MQ=60)
                                      GTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATAC                                                                                                               >  SRR19973072.448054/1‑150 (MQ=60)
                                      GTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATAC                                                                                                               <  SRR19973072.839395/150‑1 (MQ=60)
                                      GTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATAC                                                                                                               <  SRR19973072.839338/150‑1 (MQ=60)
                                       TTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACT                                                                                                              >  SRR19973072.1206033/1‑150 (MQ=60)
                                       TTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACT                                                                                                              <  SRR19973072.1732886/150‑1 (MQ=60)
                                       TTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACT                                                                                                              >  SRR19973072.1092141/1‑150 (MQ=60)
                                        TCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTC                                                                                                             <  SRR19973072.524902/150‑1 (MQ=60)
                                         CCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCC                                                                                                            <  SRR19973072.2392551/150‑1 (MQ=60)
                                           ATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGC                                                                                                          >  SRR19973072.4459222/1‑150 (MQ=60)
                                            TTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCT                                                                                                         <  SRR19973072.2071390/150‑1 (MQ=60)
                                                    GCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATT                                                                                                 <  SRR19973072.3711171/150‑1 (MQ=60)
                                                    GCGTTGAAGGCAGTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATT                                                                                                 >  SRR19973072.796615/1‑150 (MQ=60)
                                                      GTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCT                                                                                               >  SRR19973072.3105210/1‑150 (MQ=60)
                                                         GAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAAC                                                                                            >  SRR19973072.3074133/1‑150 (MQ=60)
                                                         GAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAAC                                                                                            >  SRR19973072.728181/1‑150 (MQ=60)
                                                           AGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGA                                                                                          >  SRR19973072.86397/1‑150 (MQ=60)
                                                              CATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAA                                                                                       >  SRR19973072.1867683/1‑150 (MQ=60)
                                                              CATTAAAACCTTCGAAGAACTGCTGGCGAAACAAGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGACGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACAGACTCCGCTGCAGGATTCTAACGACAA                                                                                       >  SRR19973072.1867679/1‑150 (MQ=60)
                                                                  AAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTC                                                                                   >  SRR19973072.3181525/1‑150 (MQ=60)
                                                                  AAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTC                                                                                   >  SRR19973072.3181379/1‑150 (MQ=60)
                                                                      CCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCG                                                                               >  SRR19973072.4214628/1‑150 (MQ=60)
                                                                       CTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGC                                                                              >  SRR19973072.3755103/1‑150 (MQ=60)
                                                                           GAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAAC                                                                          >  SRR19973072.4391524/1‑150 (MQ=60)
                                                                           GAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGACGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAAC                                                                          >  SRR19973072.2395721/1‑150 (MQ=60)
                                                                            AAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACA                                                                         <  SRR19973072.4374935/150‑1 (MQ=60)
                                                                            AAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACA                                                                         >  SRR19973072.3822511/1‑150 (MQ=60)
                                                                               AACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCC                                                                      >  SRR19973072.3892202/1‑150 (MQ=60)
                                                                                 CTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCTGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAG                                                                    >  SRR19973072.4356933/1‑150 (MQ=60)
                                                                                 CTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAG                                                                    >  SRR19973072.2023360/1‑150 (MQ=60)
                                                                                  TGATGGCGAAACACGGCAAAGGCTACGGTTTAGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTACCTCGCTAACATCCAGA                                                                   >  SRR19973072.2879300/1‑150 (MQ=60)
                                                                                     TGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAG                                                                <  SRR19973072.2254335/150‑1 (MQ=60)
                                                                                     TGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAG                                                                >  SRR19973072.3734344/1‑150 (MQ=60)
                                                                                        CGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACG                                                             >  SRR19973072.3504382/1‑150 (MQ=60)
                                                                                          AAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGC                                                           >  SRR19973072.2181149/1‑150 (MQ=60)
                                                                                          AAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGC                                                           >  SRR19973072.4070037/1‑150 (MQ=60)
                                                                                          AAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGC                                                           >  SRR19973072.4070030/1‑150 (MQ=60)
                                                                                           AACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCA                                                          >  SRR19973072.1486254/1‑150 (MQ=60)
                                                                                           AACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCA                                                          <  SRR19973072.1051039/150‑1 (MQ=60)
                                                                                            ACACGGCAAAGGCTACGGTTGTGAAGATTGTAAACCAACCCACGGTTCGCTGCTGCCCTCCTGCTGGAAAGAATACATTCTGAATCCGGAACATACTCCGCTGCAGGATACTAACTACAACTTACTCTCTAACATCCAGAAAGACGGCAC                                                         >  SRR19973072.1411134/1‑150 (MQ=60)
                                                                                              ACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCT                                                       <  SRR19973072.1199353/150‑1 (MQ=60)
                                                                                                 GCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTCACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACT                                                    <  SRR19973072.523070/150‑1 (MQ=60)
                                                                                                 GCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACT                                                    <  SRR19973072.522784/150‑1 (MQ=60)
                                                                                                            GGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCG                                         >  SRR19973072.2959890/1‑150 (MQ=60)
                                                                                                             GTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGC                                        <  SRR19973072.862808/150‑1 (MQ=60)
                                                                                                             GTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGC                                        >  SRR19973072.3803061/1‑150 (MQ=60)
                                                                                                               TGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGT                                      >  SRR19973072.1339012/1‑150 (MQ=60)
                                                                                                               TGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGT                                      >  SRR19973072.3667460/1‑150 (MQ=60)
                                                                                                               TGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGT                                      >  SRR19973072.2140433/1‑150 (MQ=60)
                                                                                                               TGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGT                                      >  SRR19973072.3667328/1‑150 (MQ=60)
                                                                                                                GTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTT                                     <  SRR19973072.3872118/150‑1 (MQ=60)
                                                                                                                GTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTT                                     <  SRR19973072.2223691/150‑1 (MQ=60)
                                                                                                                 TGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTC                                    <  SRR19973072.2689406/150‑1 (MQ=60)
                                                                                                                 TGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTC                                    >  SRR19973072.101431/1‑150 (MQ=60)
                                                                                                                 TGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTC                                    <  SRR19973072.1479362/150‑1 (MQ=60)
                                                                                                                 TGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTC                                    <  SRR19973072.1479267/150‑1 (MQ=60)
                                                                                                                  GAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCT                                   >  SRR19973072.2287463/1‑150 (MQ=60)
                                                                                                                   AAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTC                                  <  SRR19973072.2466809/150‑1 (MQ=60)
                                                                                                                     GTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCG                                >  SRR19973072.4478859/1‑150 (MQ=60)
                                                                                                                     GTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCG                                >  SRR19973072.4478967/1‑150 (MQ=60)
                                                                                                                          TAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACCACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGG                           >  SRR19973072.329439/1‑150 (MQ=60)
                                                                                                                          TAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGG                           >  SRR19973072.3066202/1‑150 (MQ=60)
                                                                                                                          TAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGG                           >  SRR19973072.1387159/1‑150 (MQ=60)
                                                                                                                          TAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGG                           >  SRR19973072.3065608/1‑150 (MQ=60)
                                                                                                                          TAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGG                           >  SRR19973072.123013/1‑150 (MQ=60)
                                                                                                                           AAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGT                          >  SRR19973072.1862549/1‑150 (MQ=60)
                                                                                                                           AAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGT                          >  SRR19973072.1828281/1‑150 (MQ=60)
                                                                                                                                 ACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATC                    <  SRR19973072.3259498/150‑1 (MQ=60)
                                                                                                                                      CGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCC               <  SRR19973072.4391016/150‑1 (MQ=60)
                                                                                                                                         TTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGA            >  SRR19973072.2200225/1‑150 (MQ=60)
                                                                                                                                           CGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAG          <  SRR19973072.1219641/150‑1 (MQ=60)
                                                                                                                                            GCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGG         <  SRR19973072.1260046/150‑1 (MQ=60)
                                                                                                                                             CTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGG        <  SRR19973072.4096459/150‑1 (MQ=60)
                                                                                                                                              TGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGC       >  SRR19973072.4233432/1‑150 (MQ=60)
                                                                                                                                                 TGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTGA    <  SRR19973072.3117494/150‑1 (MQ=60)
                                                                                                                                                  GCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTGAT   <  SRR19973072.4055021/150‑1 (MQ=60)
                                                                                                                                                  GCCCTCCTGCTGGAACGAATACAGTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTGAT   <  SRR19973072.4054854/150‑1 (MQ=60)
                                                                                                                                                   CCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTGATG  >  SRR19973072.1977060/1‑150 (MQ=60)
                                                                                                                                                   |                                                                                                                                                     
AACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTGATG  >  CP009273/3488819‑3489115

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: