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 F368 I0 R1 1472 214.4 7116750 95.9% 6824963 147.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 518,820 C→A 100% A35E (GCG→GAG)  rhsD → Rhs family putative polymorphic toxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273518,8200CA97.1% 534.9 / 6.2 173A35E (GCG→GAG) rhsDRhs family putative polymorphic toxin
Reads supporting (aligned to +/- strand):  ref base C (0/5);  new base A (99/69);  total (99/74)
Fisher's exact test for biased strand distribution p-value = 1.32e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 177 total aligned reads displayed.

TAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGG  >  CP009273/518677‑518962
                                                                                                                                               |                                                                                                                                              
tAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGTCGCTTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTACAGGGTTCGTCAGGTGTAAGAATTTGAGCGCCCAACGGCGTGGAGTGc                                                                                                                                            >  2:832199/1‑148 (MQ=255)
tAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGc                                                                                                                                            >  2:856765/1‑148 (MQ=255)
     caGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGgt                                                                                                                                       >  1:978309/1‑148 (MQ=255)
     caGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGgt                                                                                                                                       >  2:3089434/1‑148 (MQ=255)
      aGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGgtg                                                                                                                                      <  2:2570433/148‑1 (MQ=255)
        ggATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGgtgtg                                                                                                                                    >  1:3358690/1‑148 (MQ=255)
        ggATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGgtgtg                                                                                                                                    <  2:2876866/148‑1 (MQ=255)
         gATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGgtgtgt                                                                                                                                   <  2:1103463/148‑1 (MQ=255)
            atGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCcg                                                                                                                                >  2:1435175/1‑148 (MQ=255)
              gAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCcggg                                                                                                                              >  2:2123932/1‑148 (MQ=255)
               aaCAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCcgggc                                                                                                                             >  2:1524654/1‑148 (MQ=255)
               aaCAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCcgggc                                                                                                                             >  2:1524933/1‑148 (MQ=255)
                  aaTAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCcgggcggg                                                                                                                          <  2:1445904/148‑1 (MQ=255)
                  aaTAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCcgggcggg                                                                                                                          <  1:1379180/148‑1 (MQ=255)
                   aTAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGa                                                                                                                         <  2:3041080/148‑1 (MQ=255)
                    tAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGAt                                                                                                                        <  1:2757978/148‑1 (MQ=255)
                    tAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGAt                                                                                                                        >  1:350683/1‑148 (MQ=255)
                     aaCTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGTGCGGGATg                                                                                                                       <  2:2889887/148‑1 (MQ=255)
                     aaCTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATg                                                                                                                       >  1:2534205/1‑148 (MQ=255)
                     aaCTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATg                                                                                                                       >  1:2534235/1‑148 (MQ=255)
                       cTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGAc                                                                                                                     >  2:314862/1‑148 (MQ=255)
                       cTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGAc                                                                                                                     >  2:682599/1‑148 (MQ=255)
                       cTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGAc                                                                                                                     <  1:948401/148‑1 (MQ=255)
                         aaaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACtt                                                                                                                   <  1:2142285/148‑1 (MQ=255)
                         aaaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACtt                                                                                                                   <  1:1807228/148‑1 (MQ=255)
                          aaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGAAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTc                                                                                                                  >  2:985581/1‑148 (MQ=255)
                           aaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCg                                                                                                                 >  2:3118414/1‑148 (MQ=255)
                           aaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCg                                                                                                                 >  2:1306437/1‑148 (MQ=255)
                            aGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCgg                                                                                                                >  1:536239/1‑148 (MQ=255)
                                 aCTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAAc                                                                                                           >  2:474283/1‑148 (MQ=255)
                                        atGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGTATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTaa                                                                                                    <  2:2866458/148‑1 (MQ=255)
                                        atGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTaa                                                                                                    <  2:2866552/148‑1 (MQ=255)
                                        atGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTaa                                                                                                    >  2:712511/1‑148 (MQ=255)
                                        atGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTaa                                                                                                    >  2:1160094/1‑148 (MQ=255)
                                          gAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAAt                                                                                                  >  2:3335658/1‑148 (MQ=255)
                                          gAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAAt                                                                                                  >  2:3336140/1‑148 (MQ=255)
                                           aGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATc                                                                                                 <  1:2975206/148‑1 (MQ=255)
                                               gAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCgct                                                                                             >  1:413715/1‑148 (MQ=255)
                                                 aaaCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCgctgc                                                                                           >  2:1991409/1‑148 (MQ=255)
                                                  aaCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCgctgct                                                                                          >  1:363018/1‑148 (MQ=255)
                                                    ccAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTgg                                                                                        >  1:216515/1‑148 (MQ=255)
                                                     cAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTggg                                                                                       <  2:243942/148‑1 (MQ=255)
                                                     cAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTggg                                                                                       >  2:1651893/1‑148 (MQ=255)
                                                     cAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTggg                                                                                       <  2:1463762/148‑1 (MQ=255)
                                                           cgcgTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGaa                                                                                 <  1:841398/148‑1 (MQ=255)
                                                              gTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGt                                                                              <  1:943395/148‑1 (MQ=255)
                                                               tCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGtg                                                                             <  2:489606/148‑1 (MQ=255)
                                                                cAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGtgc                                                                            <  2:3358690/148‑1 (MQ=255)
                                                                  gggAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGtgctg                                                                          <  1:2815006/148‑1 (MQ=255)
                                                                  gggAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGtgctg                                                                          <  1:2365128/148‑1 (MQ=255)
                                                                   ggAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGtgctgc                                                                         >  1:1592912/1‑148 (MQ=255)
                                                                   ggAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGtgctgc                                                                         <  2:920520/148‑1 (MQ=255)
                                                                     agaTATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGccc                                                                       <  2:1025368/148‑1 (MQ=255)
                                                                     agaTATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGccc                                                                       <  2:1024962/148‑1 (MQ=255)
                                                                       atatGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCgg                                                                     >  2:1391631/1‑148 (MQ=255)
                                                                         atGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCg                                                                   >  2:1257896/1‑148 (MQ=255)
                                                                         atGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCg                                                                   <  2:356967/148‑1 (MQ=255)
                                                                         atGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCg                                                                   >  2:1594096/1‑148 (MQ=255)
                                                                           gACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCgag                                                                 >  1:569452/1‑148 (MQ=255)
                                                                            aCTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCgaga                                                                >  1:108751/1‑148 (MQ=255)
                                                                            aCTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCgaga                                                                >  1:2578456/1‑148 (MQ=255)
                                                                            aCTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCgaga                                                                >  1:1686105/1‑148 (MQ=255)
                                                                             cTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGAc                                                               <  2:1863754/148‑1 (MQ=255)
                                                                              tCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACg                                                              >  1:2098406/1‑148 (MQ=255)
                                                                               cAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACgg                                                             >  1:1447386/1‑148 (MQ=255)
                                                                                  tATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGAcc                                                          >  1:370078/1‑148 (MQ=255)
                                                                                  tATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGAcc                                                          >  1:2520117/1‑148 (MQ=255)
                                                                                  tATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGAcc                                                          >  1:3496521/1‑148 (MQ=255)
                                                                                     ggCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTg                                                       >  1:2575118/1‑148 (MQ=255)
                                                                                       cGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTgcg                                                     <  2:1371683/148‑1 (MQ=255)
                                                                                        ggTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTgcgc                                                    <  2:967264/148‑1 (MQ=255)
                                                                                        ggTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTgcgc                                                    <  2:2674468/148‑1 (MQ=255)
                                                                                             cATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgccc                                               >  1:2402746/1‑148 (MQ=255)
                                                                                             cATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgccc                                               >  2:139745/1‑148 (MQ=255)
                                                                                                tGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggc                                            >  1:1208274/1‑148 (MQ=255)
                                                                                                tGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggc                                            <  1:712511/148‑1 (MQ=255)
                                                                                                   ccAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggcccg                                         <  1:3028337/148‑1 (MQ=255)
                                                                                                   ccAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggcccg                                         <  1:1862372/148‑1 (MQ=255)
                                                                                                    cAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGc                                        <  2:2245373/148‑1 (MQ=255)
                                                                                                     aGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCt                                       >  2:2483775/1‑148 (MQ=255)
                                                                                                     aGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCt                                       >  2:2292704/1‑148 (MQ=255)
                                                                                                       ggTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGc                                     >  2:2134333/1‑148 (MQ=255)
                                                                                                             gCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCa                               <  1:2894042/148‑1 (MQ=255)
                                                                                                              cAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCAt                              <  1:1435355/148‑1 (MQ=255)
                                                                                                                ggTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATtc                            <  1:1840835/148‑1 (MQ=255)
                                                                                                                 gtgtAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATtct                           >  1:451254/1‑148 (MQ=255)
                                                                                                                   gtAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATtctct                         >  2:2215765/1‑148 (MQ=255)
                                                                                                                     aaGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTcc                       <  1:1870934/148‑1 (MQ=255)
                                                                                                                     aaGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTcc                       <  1:70266/148‑1 (MQ=255)
                                                                                                                        aaTTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGc                    >  1:3312170/1‑148 (MQ=255)
                                                                                                                         aTTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCa                   >  2:1687970/1‑148 (MQ=255)
                                                                                                                          ttGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCAc                  <  2:1236223/148‑1 (MQ=255)
                                                                                                                          ttGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCAc                  >  1:882305/1‑148 (MQ=255)
                                                                                                                           tGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCAcc                 >  1:407648/1‑148 (MQ=255)
                                                                                                                           tGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTAATTCTCTCCCGCAcc                 <  1:89944/148‑1 (MQ=255)
                                                                                                                                    cACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACagcagc        <  2:1592912/148‑1 (MQ=255)
                                                                                                                                    cACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACagcagc        >  1:131816/1‑148 (MQ=255)
                                                                                                                                    cACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACagcagc        >  1:2949423/1‑148 (MQ=255)
                                                                                                                                         gcgtggagtgCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCg   >  2:1682004/1‑148 (MQ=255)
                                                                                                                                         gcgtggagtgCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCg   >  2:3439703/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                              
TAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGG  >  CP009273/518677‑518962

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
TGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACG  >  CP009273/518671‑518965
                                                                                                                                                     |                                                                                                                                                 
TGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGA                                                                                                                                                   >  SRR19973067.856821/1‑150 (MQ=60)
  TTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGT                                                                                                                                                 <  SRR19973067.942909/150‑1 (MQ=60)
       AACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCG                                                                                                                                            >  SRR19973067.1876468/1‑150 (MQ=60)
         CACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGT                                                                                                                                          >  SRR19973067.985020/1‑150 (MQ=60)
         CACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGT                                                                                                                                          >  SRR19973067.985003/1‑150 (MQ=60)
           CAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGT                                                                                                                                        >  SRR19973067.2381462/1‑150 (MQ=60)
            AGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTG                                                                                                                                       >  SRR19973067.3380013/1‑150 (MQ=60)
               GATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCC                                                                                                                                    <  SRR19973067.3339785/150‑1 (MQ=60)
                        AATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGAT                                                                                                                           <  SRR19973067.1388583/150‑1 (MQ=60)
                        AATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGAT                                                                                                                           >  SRR19973067.852461/1‑150 (MQ=60)
                        AATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGAT                                                                                                                           >  SRR19973067.353108/1‑150 (MQ=60)
                         ATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATG                                                                                                                          >  SRR19973067.2550839/1‑150 (MQ=60)
                         ATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATG                                                                                                                          >  SRR19973067.2550824/1‑150 (MQ=60)
                         ATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATG                                                                                                                          >  SRR19973067.2550809/1‑150 (MQ=60)
                          TAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGA                                                                                                                         <  SRR19973067.2775801/150‑1 (MQ=60)
                             CTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTT                                                                                                                      <  SRR19973067.954889/150‑1 (MQ=60)
                             CTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGAGCTCCGTGTGTCCGGGCGGGATGACTT                                                                                                                      >  SRR19973067.3088725/1‑150 (MQ=60)
                              TAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTC                                                                                                                     >  SRR19973067.2554948/1‑150 (MQ=60)
                               AAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTNCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCG                                                                                                                    <  SRR19973067.1819427/150‑1 (MQ=60)
                               AAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCG                                                                                                                    <  SRR19973067.2156606/150‑1 (MQ=60)
                                AAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGG                                                                                                                   >  SRR19973067.539822/1‑150 (MQ=60)
                                AAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGG                                                                                                                   >  SRR19973067.539863/1‑150 (MQ=60)
                                               TGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATC                                                                                                    <  SRR19973067.525156/150‑1 (MQ=60)
                                                 AGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCG                                                                                                  >  SRR19973067.359442/1‑150 (MQ=60)
                                                 AGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCG                                                                                                  <  SRR19973067.2993961/150‑1 (MQ=60)
                                                 AGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCG                                                                                                  <  SRR19973067.2993949/150‑1 (MQ=60)
                                                   CGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCT                                                                                                >  SRR19973067.416590/1‑150 (MQ=60)
                                                    GGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTG                                                                                               <  SRR19973067.346537/150‑1 (MQ=60)
                                                      AAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCT                                                                                             >  SRR19973067.365541/1‑150 (MQ=60)
                                                      AAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCT                                                                                             >  SRR19973067.365558/1‑150 (MQ=60)
                                                       AAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTG                                                                                            >  SRR19973067.2509491/1‑150 (MQ=60)
                                                        AACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGG                                                                                           >  SRR19973067.218113/1‑150 (MQ=60)
                                                        AACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGG                                                                                           >  SRR19973067.2950705/1‑150 (MQ=60)
                                                                 CGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGG                                                                                  <  SRR19973067.847107/150‑1 (MQ=60)
                                                                  GCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGT                                                                                 <  SRR19973067.2846524/150‑1 (MQ=60)
                                                                   CGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTG                                                                                >  SRR19973067.1797782/1‑150 (MQ=60)
                                                                    GTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGC                                                                               <  SRR19973067.950746/150‑1 (MQ=60)
                                                                    GTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGC                                                                               <  SRR19973067.949831/150‑1 (MQ=60)
                                                                       AGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGC                                                                            >  SRR19973067.1603828/1‑150 (MQ=60)
                                                                       AGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGC                                                                            >  SRR19973067.1603832/1‑150 (MQ=60)
                                                                        GGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCC                                                                           <  SRR19973067.2833136/150‑1 (MQ=60)
                                                                        GGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCC                                                                           <  SRR19973067.2380532/150‑1 (MQ=60)
                                                                          GAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCG                                                                         <  SRR19973067.2135109/150‑1 (MQ=60)
                                                                           AGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGG                                                                        >  SRR19973067.2811486/1‑150 (MQ=60)
                                                                               ATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAG                                                                    <  SRR19973067.122366/150‑1 (MQ=60)
                                                                               ATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAG                                                                    >  SRR19973067.573148/1‑150 (MQ=60)
                                                                                TGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGA                                                                   >  SRR19973067.2595170/1‑150 (MQ=60)
                                                                                TGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGA                                                                   >  SRR19973067.109811/1‑150 (MQ=60)
                                                                                TGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGA                                                                   >  SRR19973067.1697724/1‑150 (MQ=60)
                                                                                TGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGA                                                                   >  SRR19973067.2595257/1‑150 (MQ=60)
                                                                                  ACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACG                                                                 <  SRR19973067.626843/150‑1 (MQ=60)
                                                                                  ACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACG                                                                 >  SRR19973067.2112310/1‑150 (MQ=60)
                                                                                   CTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGG                                                                >  SRR19973067.1457097/1‑150 (MQ=60)
                                                                                      AGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTCAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACC                                                             >  SRR19973067.3518702/1‑150 (MQ=60)
                                                                                      AGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACC                                                             >  SRR19973067.2536642/1‑150 (MQ=60)
                                                                                      AGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACC                                                             >  SRR19973067.372655/1‑150 (MQ=60)
                                                                                      AGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACC                                                             >  SRR19973067.372657/1‑150 (MQ=60)
                                                                                         ATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTG                                                          >  SRR19973067.2591911/1‑150 (MQ=60)
                                                                                                 CCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCC                                                  >  SRR19973067.2418432/1‑150 (MQ=60)
                                                                                                    ATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGC                                               >  SRR19973067.1216215/1‑150 (MQ=60)
                                                                                                      TGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCC                                             <  SRR19973067.717447/150‑1 (MQ=60)
                                                                                                       GTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCG                                            <  SRR19973067.1572461/150‑1 (MQ=60)
                                                                                                         CCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCT                                          <  SRR19973067.1875067/150‑1 (MQ=60)
                                                                                                         CCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCT                                          <  SRR19973067.3047617/150‑1 (MQ=60)
                                                                                                         CCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCCGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCT                                          <  SRR19973067.100434/150‑1 (MQ=60)
                                                                                                           AGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGC                                        <  SRR19973067.833209/150‑1 (MQ=60)
                                                                                                             GGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCG                                      <  SRR19973067.2207918/150‑1 (MQ=60)
                                                                                                                   GCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATT                                <  SRR19973067.2912522/150‑1 (MQ=60)
                                                                                                                    CAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTC                               <  SRR19973067.1444837/150‑1 (MQ=60)
                                                                                                                    CAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTC                               <  SRR19973067.1445018/150‑1 (MQ=60)
                                                                                                                     AGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCT                              >  SRR19973067.453010/1‑150 (MQ=60)
                                                                                                                     AGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCT                              >  SRR19973067.454385/1‑150 (MQ=60)
                                                                                                                      GGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTC                             <  SRR19973067.1853393/150‑1 (MQ=60)
                                                                                                                      GGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTC                             <  SRR19973067.1853847/150‑1 (MQ=60)
                                                                                                                           AAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCG                        <  SRR19973067.71083/150‑1 (MQ=60)
                                                                                                                           AAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCG                        <  SRR19973067.1883695/150‑1 (MQ=60)
                                                                                                                            AGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGC                       >  SRR19973067.3333096/1‑150 (MQ=60)
                                                                                                                              AATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCAC                     >  SRR19973067.888180/1‑150 (MQ=60)
                                                                                                                               ATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCCCACC                    >  SRR19973067.410487/1‑150 (MQ=60)
                                                                                                                                 TGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTA                  <  SRR19973067.1375969/150‑1 (MQ=60)
                                                                                                                                 TGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTAATTCTCTCCCGCACCTA                  <  SRR19973067.90872/150‑1 (MQ=60)
                                                                                                                                 TGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTAATTCTCTCCCGCACCTA                  <  SRR19973067.90907/150‑1 (MQ=60)
                                                                                                                                        CCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGC           >  SRR19973067.133005/1‑150 (MQ=60)
                                                                                                                                        CCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGC           >  SRR19973067.1941898/1‑150 (MQ=60)
                                                                                                                                        CCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGC           >  SRR19973067.2968086/1‑150 (MQ=60)
                                                                                                                                                 GTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACG  >  SRR19973067.1563066/1‑150 (MQ=60)
                                                                                                                                                 GTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACG  >  SRR19973067.1395085/1‑150 (MQ=60)
                                                                                                                                                 GTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACG  >  SRR19973067.1901442/1‑150 (MQ=60)
                                                                                                                                                     |                                                                                                                                                 
TGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACG  >  CP009273/518671‑518965

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: ‑