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 F204 I0 R1 1155 255.8 8356456 95.7% 7997128 150.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,204,284 C→A 89.9% D98Y (GAT→TAT)  tfaE ← e14 prophage; putative tail fiber assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,204,2840CA89.9% 321.8 / 33.1 119D98Y (GAT→TAT) tfaEe14 prophage; putative tail fiber assembly protein
Reads supporting (aligned to +/- strand):  ref base C (7/5);  new base A (51/56);  total (58/61)
Fisher's exact test for biased strand distribution p-value = 5.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

Only 100 of 119 total aligned reads displayed.

TTGCAGCATCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAAC  >  CP009273/1204137‑1204430
                                                                                                                                                   |                                                                                                                                                  
ttGCAGCATCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAAc                                                                                                                                                  <  2:2920600/150‑1 (MQ=255)
 tGCAGCATCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAAcc                                                                                                                                                 >  1:4039909/1‑150 (MQ=255)
   cagcaTCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCg                                                                                                                                               >  2:2977081/1‑150 (MQ=255)
       atcaatcaGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGc                                                                                                                                           <  1:286244/150‑1 (MQ=255)
       atcaatcaGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGGATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGc                                                                                                                                           >  2:3965498/1‑150 (MQ=255)
         caatcaGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCtt                                                                                                                                         >  1:2509508/1‑150 (MQ=255)
          aatcaGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCttt                                                                                                                                        <  1:2362467/150‑1 (MQ=255)
          aatcaGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCttt                                                                                                                                        <  2:1079767/150‑1 (MQ=255)
          aatcaGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCttt                                                                                                                                        <  2:2880409/150‑1 (MQ=255)
          aatcaGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCttt                                                                                                                                        <  2:3493959/150‑1 (MQ=255)
           atcaGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGGTGTATTTTCAGGGTAATAACCCGGAGCTTTg                                                                                                                                       >  2:2015961/1‑150 (MQ=255)
              aGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGAtt                                                                                                                                    >  1:734425/1‑150 (MQ=255)
              aGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGAtt                                                                                                                                    <  2:3456196/150‑1 (MQ=255)
              aGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGAtt                                                                                                                                    <  2:3456000/150‑1 (MQ=255)
              aGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGAtt                                                                                                                                    >  1:3489602/1‑150 (MQ=255)
                tGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTc                                                                                                                                  >  2:1370190/1‑150 (MQ=255)
                 gACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCt                                                                                                                                 >  1:3114156/1‑150 (MQ=255)
                  aCTGGCGCTGTGCTTGTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAATGGGGCTATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCtt                                                                                                                                <  1:3230106/150‑1 (MQ=255)
                  aCTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATACAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCTGAGCTTTGATTTCtt                                                                                                                                >  2:2803618/1‑150 (MQ=255)
                  aCTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCtt                                                                                                                                <  1:2856708/150‑1 (MQ=255)
                   cTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCttt                                                                                                                               <  2:2813398/150‑1 (MQ=255)
                       cgcTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGAt                                                                                                                           <  1:1934547/150‑1 (MQ=255)
                              cTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAg                                                                                                                    <  2:2963001/150‑1 (MQ=255)
                              cTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAg                                                                                                                    >  1:456584/1‑150 (MQ=255)
                              cTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAg                                                                                                                    >  1:1079288/1‑150 (MQ=255)
                              cTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAg                                                                                                                    >  1:1450434/1‑150 (MQ=255)
                                 cTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGttt                                                                                                                 >  1:173664/1‑150 (MQ=255)
                                   gCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTg                                                                                                               <  2:2155422/150‑1 (MQ=255)
                                   gCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTg                                                                                                               <  2:2155485/150‑1 (MQ=255)
                                    ccGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTgg                                                                                                              >  1:2833025/1‑150 (MQ=255)
                                    acgagTCTACTGCGGCGAAAAGCTGTGCTTCAGAATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAAGGGCGATAGTGGTATTATTTTCAGGGTAATAACCCGTAGCATTGATTACATTTGATTCTCAAGTTTTgg                                                                                                              >  2:1490090/5‑150 (MQ=25)
                                      gcgTCTACTGCGGCGCTATGCTGTGATTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGTGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTg                                                                                                            <  1:430977/150‑1 (MQ=255)
                                           tACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCtata                                                                                                       <  2:2577898/150‑1 (MQ=255)
                                           tACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCtata                                                                                                       >  1:490182/1‑150 (MQ=255)
                                           tACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCtata                                                                                                       >  1:3939708/1‑150 (MQ=255)
                                              tGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGAc                                                                                                    <  2:887544/150‑1 (MQ=255)
                                              tGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGAc                                                                                                    >  1:2626349/1‑150 (MQ=255)
                                                    gcTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTc                                                                                              >  1:1492800/1‑150 (MQ=255)
                                                               cTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGt                                                                                   <  2:1461668/150‑1 (MQ=255)
                                                               cTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGt                                                                                   <  2:1461955/150‑1 (MQ=255)
                                                                ttCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTc                                                                                  <  2:1740868/150‑1 (MQ=255)
                                                                 tCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCt                                                                                 <  1:4156145/150‑1 (MQ=255)
                                                                 tCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCt                                                                                 <  1:4156215/150‑1 (MQ=255)
                                                                   aGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTgg                                                                               <  1:2682643/150‑1 (MQ=255)
                                                                     tATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGAGTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGca                                                                             >  2:1223310/1‑150 (MQ=255)
                                                                      aTCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGcac                                                                            <  2:1725623/150‑1 (MQ=255)
                                                                       tCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGcaca                                                                           >  1:1587942/1‑150 (MQ=255)
                                                                        cGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACat                                                                          >  1:1805297/1‑150 (MQ=255)
                                                                        cGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACat                                                                          >  1:1804925/1‑150 (MQ=255)
                                                                        cGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTGGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGGTTCACCGCGATGGTCTGGcacct                                                                          >  2:2840837/1‑148 (MQ=255)
                                                                         ggTCACCCATTTCTCACAATCCCAATTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATTGACCGTTTCACAGCGATGGTCTGGCACata                                                                         >  2:3551468/1‑150 (MQ=255)
                                                                          gTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGTTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACatat                                                                        <  2:2009130/150‑1 (MQ=255)
                                                                           tCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCGCATAtt                                                                       >  1:2429618/1‑150 (MQ=255)
                                                                                     ctcACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGtt                                                             <  2:2245202/150‑1 (MQ=255)
                                                                                     ctcACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGtt                                                             <  2:2245194/150‑1 (MQ=255)
                                                                                     ctcACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGtt                                                             <  2:2245196/150‑1 (MQ=255)
                                                                                      tcACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTa                                                            >  2:1114967/1‑150 (MQ=255)
                                                                                          cATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAAt                                                        <  1:2473495/150‑1 (MQ=255)
                                                                                           aTCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATc                                                       <  1:3959671/150‑1 (MQ=255)
                                                                                           aTCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATc                                                       <  1:3959816/150‑1 (MQ=255)
                                                                                            tCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCt                                                      <  2:3467903/150‑1 (MQ=255)
                                                                                            tCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCt                                                      <  2:3467892/150‑1 (MQ=255)
                                                                                              ccATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGc                                                    <  1:1302748/150‑1 (MQ=255)
                                                                                               cATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCa                                                   >  2:2127534/1‑150 (MQ=255)
                                                                                               cATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCa                                                   >  2:2219134/1‑150 (MQ=255)
                                                                                                 tttATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGa                                                 <  1:3038726/150‑1 (MQ=255)
                                                                                                      cGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGc                                            >  1:163708/1‑150 (MQ=255)
                                                                                                           ggAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGAtt                                       >  1:2054212/1‑150 (MQ=255)
                                                                                                               aTAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCAt                                   <  1:213425/150‑1 (MQ=255)
                                                                                                               aTAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCAt                                   <  1:213452/150‑1 (MQ=255)
                                                                                                                 aaaGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATaa                                 >  2:809421/1‑150 (MQ=255)
                                                                                                                 aaaGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATaa                                 >  2:1849163/1‑150 (MQ=255)
                                                                                                                   aGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAAcc                               <  2:826915/150‑1 (MQ=255)
                                                                                                                     gggCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAg                             >  2:3150326/1‑150 (MQ=255)
                                                                                                                     gggCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAg                             >  2:3150314/1‑150 (MQ=255)
                                                                                                                      ggCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGc                            >  1:710680/1‑150 (MQ=255)
                                                                                                                       gCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGcc                           >  1:3270379/1‑150 (MQ=255)
                                                                                                                         taTAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCtt                         <  2:1667503/149‑1 (MQ=255)
                                                                                                                         gATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCtt                         >  1:2155424/1‑150 (MQ=255)
                                                                                                                         gATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCtt                         >  2:1790463/1‑150 (MQ=255)
                                                                                                                         gATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCtt                         >  2:1790394/1‑150 (MQ=255)
                                                                                                                           tAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTAt                       <  1:2268497/150‑1 (MQ=255)
                                                                                                                             gTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTAtgt                     <  1:3726442/150‑1 (MQ=255)
                                                                                                                             gTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTAtgt                     >  2:4003375/1‑150 (MQ=255)
                                                                                                                             gTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTAtgt                     <  1:3726463/150‑1 (MQ=255)
                                                                                                                             gTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTAtgt                     >  2:3551823/1‑150 (MQ=255)
                                                                                                                             gTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTAtgt                     <  1:3726459/150‑1 (MQ=255)
                                                                                                                                gTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAAACAGCCTTATGTGTa                  >  1:2148279/1‑150 (MQ=255)
                                                                                                                                  tGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTAcc                <  2:382696/150‑1 (MQ=255)
                                                                                                                                   gTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCAGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCa               >  2:2958049/1‑150 (MQ=255)
                                                                                                                                      ttttCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAggg            <  1:3381132/150‑1 (MQ=255)
                                                                                                                                         tCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCa         <  1:3086838/150‑1 (MQ=255)
                                                                                                                                          cAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCAt        >  2:2989417/1‑150 (MQ=255)
                                                                                                                                           aGGGTAATCACCCGGAGCTGTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTACCATGAGTTAAAATCTGCAGAATGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATc       >  2:1476022/1‑150 (MQ=255)
                                                                                                                                           aGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATc       >  2:1468094/1‑150 (MQ=255)
                                                                                                                                            gggTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCt      <  2:4063094/150‑1 (MQ=255)
                                                                                                                                            gggTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCt      <  2:557746/150‑1 (MQ=255)
                                                                                                                                             ggTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTa     <  1:2990435/150‑1 (MQ=255)
                                                                                                                                             ggTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTa     <  1:2990431/150‑1 (MQ=255)
                                                                                                                                                aaTAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAAc  <  1:918957/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTGCAGCATCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAAC  >  CP009273/1204137‑1204430

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
CATTGCAGCATCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGG  >  CP009273/1204135‑1204433
                                                                                                                                                     |                                                                                                                                                     
CATTGCAGCATCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATA                                                                                                                                                       >  SRR19973074.558843/1‑150 (MQ=60)
CATTGCAGCATCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATA                                                                                                                                                       >  SRR19973074.4070608/1‑150 (MQ=60)
CATTGCAGCATCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATA                                                                                                                                                       >  SRR19973074.4069723/1‑150 (MQ=60)
 ATTGCAGCATCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAA                                                                                                                                                      <  SRR19973074.205035/150‑1 (MQ=60)
 ATTGCAGCATCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAA                                                                                                                                                      <  SRR19973074.204495/150‑1 (MQ=60)
 ATTGCAGCATCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAA                                                                                                                                                      >  SRR19973074.3901061/1‑150 (MQ=60)
   TGCAGCATCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACC                                                                                                                                                    >  SRR19973074.4047377/1‑150 (MQ=60)
         ATCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGC                                                                                                                                              <  SRR19973074.286784/150‑1 (MQ=60)
           CAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTT                                                                                                                                            >  SRR19973074.2514181/1‑150 (MQ=60)
            AATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTT                                                                                                                                           <  SRR19973074.2366849/150‑1 (MQ=60)
            AATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTT                                                                                                                                           <  SRR19973074.2366864/150‑1 (MQ=60)
                AGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATT                                                                                                                                       >  SRR19973074.3496003/1‑150 (MQ=60)
                AGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATT                                                                                                                                       >  SRR19973074.735802/1‑150 (MQ=60)
                   GACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCT                                                                                                                                    >  SRR19973074.3119869/1‑150 (MQ=60)
                    ACTGGCGCTGTGCTTGTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAATGGGGCTATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTT                                                                                                                                   <  SRR19973074.3236080/150‑1 (MQ=60)
                    ACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTT                                                                                                                                   <  SRR19973074.2862010/150‑1 (MQ=60)
                         CGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGAT                                                                                                                              <  SRR19973074.1938244/150‑1 (MQ=60)
                                CTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAG                                                                                                                       >  SRR19973074.457478/1‑150 (MQ=60)
                                CTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAG                                                                                                                       >  SRR19973074.1081304/1‑150 (MQ=60)
                                CTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAG                                                                                                                       >  SRR19973074.1453149/1‑150 (MQ=60)
                                   CTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTT                                                                                                                    >  SRR19973074.858650/1‑150 (MQ=60)
                                   CTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTT                                                                                                                    >  SRR19973074.3314994/1‑150 (MQ=60)
                                   CTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTT                                                                                                                    >  SRR19973074.174045/1‑150 (MQ=60)
                                      CCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGG                                                                                                                 >  SRR19973074.2838266/1‑150 (MQ=60)
                                        GCGTCTACTGCGGCGCTATGCTGTGATTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGTGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTG                                                                                                               <  SRR19973074.431808/150‑1 (MQ=60)
                                             TACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATA                                                                                                          >  SRR19973074.3946967/1‑150 (MQ=60)
                                             TACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATA                                                                                                          >  SRR19973074.491167/1‑150 (MQ=60)
                                                TGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGAC                                                                                                       >  SRR19973074.2631262/1‑150 (MQ=60)
                                                      GCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTC                                                                                                 >  SRR19973074.1495611/1‑150 (MQ=60)
                                                                   TCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCT                                                                                    <  SRR19973074.4163865/150‑1 (MQ=60)
                                                                   TCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCT                                                                                    <  SRR19973074.3823481/150‑1 (MQ=60)
                                                                   TCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCT                                                                                    <  SRR19973074.4163795/150‑1 (MQ=60)
                                                                   TCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCT                                                                                    <  SRR19973074.4163890/150‑1 (MQ=60)
                                                                     AGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGG                                                                                  <  SRR19973074.2687628/150‑1 (MQ=60)
                                                                     AGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGG                                                                                  <  SRR19973074.2687716/150‑1 (MQ=60)
                                                                         TCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACA                                                                              >  SRR19973074.1590963/1‑150 (MQ=60)
                                                                          CGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACAT                                                                             >  SRR19973074.1808359/1‑150 (MQ=60)
                                                                          CGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACAT                                                                             >  SRR19973074.1808731/1‑150 (MQ=60)
                                                                             TCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCGCATATT                                                                          >  SRR19973074.2434112/1‑150 (MQ=60)
                                                                                            CATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAAT                                                           <  SRR19973074.2478084/150‑1 (MQ=60)
                                                                                             ATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATC                                                          <  SRR19973074.3967114/150‑1 (MQ=60)
                                                                                             ATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATC                                                          <  SRR19973074.3966969/150‑1 (MQ=60)
                                                                                                CCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGC                                                       <  SRR19973074.1305208/150‑1 (MQ=60)
                                                                                                   TTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGA                                                    <  SRR19973074.3044327/150‑1 (MQ=60)
                                                                                                   TTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGA                                                    <  SRR19973074.3044347/150‑1 (MQ=60)
                                                                                                        CGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGC                                               >  SRR19973074.164074/1‑150 (MQ=60)
                                                                                                        CGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGC                                               >  SRR19973074.164199/1‑150 (MQ=60)
                                                                                                             GGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATT                                          >  SRR19973074.2058074/1‑150 (MQ=60)
                                                                                                                 ATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCAT                                      <  SRR19973074.213866/150‑1 (MQ=60)
                                                                                                                 ATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCAT                                      <  SRR19973074.213893/150‑1 (MQ=60)
                                                                                                                        GGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGC                               >  SRR19973074.712001/1‑150 (MQ=60)
                                                                                                                         GCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCC                              >  SRR19973074.3276444/1‑150 (MQ=60)
                                                                                                                           GATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTT                            >  SRR19973074.2159481/1‑150 (MQ=60)
                                                                                                                             TAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTAT                          <  SRR19973074.2272757/150‑1 (MQ=60)
                                                                                                                               GTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGT                        <  SRR19973074.3733334/150‑1 (MQ=60)
                                                                                                                               GTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGT                        <  SRR19973074.3733330/150‑1 (MQ=60)
                                                                                                                               GTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGT                        <  SRR19973074.3733313/150‑1 (MQ=60)
                                                                                                                               GTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGT                        >  SRR19973074.2111786/1‑150 (MQ=60)
                                                                                                                                  GTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAAACAGCCTTATGTGTA                     >  SRR19973074.2152318/1‑150 (MQ=60)
                                                                                                                                        TTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGG               <  SRR19973074.3387394/150‑1 (MQ=60)
                                                                                                                                           TCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCA            <  SRR19973074.3092512/150‑1 (MQ=60)
                                                                                                                                               GGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTA        <  SRR19973074.2995990/150‑1 (MQ=60)
                                                                                                                                               GGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTA        <  SRR19973074.2995994/150‑1 (MQ=60)
                                                                                                                                                  AATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAAC     <  SRR19973074.920709/150‑1 (MQ=60)
                                                                                                                                                    TAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAG   >  SRR19973074.358600/1‑150 (MQ=60)
                                                                                                                                                    TAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAG   >  SRR19973074.2952972/1‑150 (MQ=60)
                                                                                                                                                     AACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGG  >  SRR19973074.2240532/1‑150 (MQ=60)
                                                                                                                                                     |                                                                                                                                                     
CATTGCAGCATCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGG  >  CP009273/1204135‑1204433

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