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 3,509,102 G→T 93.9% S57R (AGC→AGA)  dam ← DNA adenine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,509,1020GT93.9% 468.9 / 20.2 167S57R (AGC→AGA) damDNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (8/2);  new base T (83/73);  total (91/76)
Fisher's exact test for biased strand distribution p-value = 1.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 170 total aligned reads displayed.

GGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCAC  >  CP009273/3508957‑3509246
                                                                                                                                                 |                                                                                                                                                
ggaacggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACACTGTTATAGAGACTGATCAGGTctct                                                                                                                                                >  2:831423/1‑148 (MQ=255)
ggaacggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTctct                                                                                                                                                <  1:2312692/148‑1 (MQ=255)
  aacggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTAt                                                                                                                                              <  1:3410006/148‑1 (MQ=255)
    cggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTg                                                                                                                                            <  1:594920/148‑1 (MQ=255)
       aTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCCCAATGTTATAGAGACTGATCAGGTCTCTATTGata                                                                                                                                         >  1:416058/1‑148 (MQ=255)
       aTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGata                                                                                                                                         <  1:875430/148‑1 (MQ=255)
       aTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGata                                                                                                                                         >  1:173934/1‑148 (MQ=255)
       aTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGata                                                                                                                                         >  1:173982/1‑148 (MQ=255)
           tGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCgg                                                                                                                                     >  1:841280/1‑148 (MQ=255)
           tGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCgg                                                                                                                                     >  1:840111/1‑148 (MQ=255)
           tGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCgg                                                                                                                                     <  2:1429672/148‑1 (MQ=255)
           tGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCgg                                                                                                                                     <  2:3197989/148‑1 (MQ=255)
             gCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCa                                                                                                                                   >  1:3138768/1‑148 (MQ=255)
             gCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCa                                                                                                                                   >  1:3138583/1‑148 (MQ=255)
                   gttgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGata                                                                                                                             <  2:2928548/148‑1 (MQ=255)
                    ttgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGatat                                                                                                                            <  2:1115064/148‑1 (MQ=255)
                       aaCTCTTCGCGGAACTGATATTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAAc                                                                                                                         >  2:2497175/1‑148 (MQ=255)
                        aCTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACg                                                                                                                        <  1:1478786/148‑1 (MQ=255)
                         ctctTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACga                                                                                                                       >  2:1479064/1‑148 (MQ=255)
                          tctTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACgag                                                                                                                      >  1:2193900/1‑148 (MQ=255)
                            ttCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGaa                                                                                                                    <  1:2935518/148‑1 (MQ=255)
                             tCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGaaa                                                                                                                   >  2:239310/1‑148 (MQ=255)
                             tCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGaaa                                                                                                                   >  2:2033545/1‑148 (MQ=255)
                               gcgGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAg                                                                                                                 <  1:1318447/148‑1 (MQ=255)
                                   aaCTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCgg                                                                                                             <  1:1010503/148‑1 (MQ=255)
                                   aaCTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCgg                                                                                                             <  2:3350443/148‑1 (MQ=255)
                                   aaCTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCgg                                                                                                             <  2:421725/148‑1 (MQ=255)
                                     cTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGgtg                                                                                                           <  1:265079/148‑1 (MQ=255)
                                     cTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGgtg                                                                                                           <  2:1181628/148‑1 (MQ=255)
                                      tGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGgtgt                                                                                                          >  2:1767138/1‑148 (MQ=255)
                                      tGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGgtgt                                                                                                          >  2:1767087/1‑148 (MQ=255)
                                               aCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAaca                                                                                                 >  1:145925/1‑148 (MQ=255)
                                                ccTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAacac                                                                                                >  2:1159602/1‑148 (MQ=255)
                                                  tCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCg                                                                                              <  2:2732028/148‑1 (MQ=255)
                                                   cGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGa                                                                                             <  2:2181723/148‑1 (MQ=255)
                                                   cGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGa                                                                                             <  2:2181694/148‑1 (MQ=255)
                                                      cgcAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGAccc                                                                                          >  1:515075/1‑148 (MQ=255)
                                                       gcAATTTTTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACTTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTTATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCg                                                                                         <  2:3022394/148‑1 (MQ=255)
                                                         aaTTTGTTTCGTGAAAAAACAGCTCGCTTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACATTGTTATAGAGACTGATCAGGTCTATATTGATATCGGCATGTATATAACGATAAATGTCGGTGTTGAGAAACACCGACCCGGc                                                                                       <  2:1944813/148‑1 (MQ=255)
                                                         aaTTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGc                                                                                       >  1:201493/1‑148 (MQ=255)
                                                          aTTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCa                                                                                      >  1:656869/1‑148 (MQ=255)
                                                           tttGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCAc                                                                                     >  1:1077771/1‑148 (MQ=255)
                                                                     gAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAgg                                                                           >  2:3055408/1‑148 (MQ=255)
                                                                     gAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAgg                                                                           >  2:3055764/1‑148 (MQ=255)
                                                                       acaaacaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCt                                                                         >  1:1980946/1‑148 (MQ=255)
                                                                       acaaacaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCACCTTCACAATGTTATAGAGACGGATCAGGTCTCTATTGATATCGGCAAGGATATAACGACAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCt                                                                         <  2:97016/148‑1 (MQ=255)
                                                                        caaacaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTc                                                                        >  1:2773125/1‑148 (MQ=255)
                                                                          aacaGCTCGCGTGCGGCCTGTACTTACTAATCAGTACGCATCTTCACATTGTTATAGATACTTAACAGGTCTCTATTTATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCaa                                                                      <  2:1077771/148‑1 (MQ=255)
                                                                          aacaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCaa                                                                      >  1:472003/1‑148 (MQ=255)
                                                                               cTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGa                                                                 >  1:2309448/1‑148 (MQ=255)
                                                                                tCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGAc                                                                >  2:382216/1‑148 (MQ=255)
                                                                                tCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGAc                                                                >  2:382208/1‑148 (MQ=255)
                                                                                 cgcgTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACa                                                               >  1:325326/1‑148 (MQ=255)
                                                                                  gcgTGCGCCCTGTACGTACACATCAGTAAGCATCTTCACAATGTTATAGAGACTCAACAGGTCTCTATTGATATCGGCAAGGATATACCGATAAAAGTCGTTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACAt                                                              <  2:2291089/148‑1 (MQ=255)
                                                                                     tGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCTGGTCTCTTTTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCg                                                           <  2:1987964/148‑1 (MQ=255)
                                                                                          ccTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCtt                                                      >  2:3075709/1‑148 (MQ=255)
                                                                                                cgtacTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATTGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTCCAAAAGGCTCAACCAGACATTCGCCCTTGGGCaa                                                >  2:3281681/1‑148 (MQ=255)
                                                                                                cgtacTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCaa                                                >  1:697495/1‑148 (MQ=255)
                                                                                                 gtacTCATCATAACGAATCTTCACAATTTTATAGACACTTATCAGGTATATATTTATATCGGCAAGTATATAAGGATAAAAGTGGGTGTTGACAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCaaa                                               <  1:886670/148‑1 (MQ=255)
                                                                                                 gtacTCATCAGTACGCATCTTCACAATGTTATAGAGACTGAGCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCaaa                                               <  2:3397749/148‑1 (MQ=255)
                                                                                                  tacTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAAt                                              >  2:2579892/1‑148 (MQ=255)
                                                                                                    cTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGc                                            >  2:1128589/1‑148 (MQ=255)
                                                                                                     tcatcaGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGcc                                           >  1:2239109/1‑148 (MQ=255)
                                                                                                       atcaGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGt                                         >  1:3546593/1‑148 (MQ=255)
                                                                                                            tACGCATATTCAAAATGTTATAGAGACTGATCAGGTCTATATTGATAGCGGCAAGTATATACCGATAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAat                                    <  2:656869/148‑1 (MQ=255)
                                                                                                               gCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatc                                 >  1:2707850/1‑148 (MQ=255)
                                                                                                                cATCTTCACACTGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGACAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatca                                <  2:656074/148‑1 (MQ=255)
                                                                                                                 aTCTTCACAATGTTATAGAGAATTATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatcat                               <  2:874316/148‑1 (MQ=255)
                                                                                                                 aTCTCCACAATGTTATAGATACTTATCAGGGCGATTTTTATATCGGCACGGATATAACGAGAAATGTCGGTTTTGACAAACCCCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatcat                               <  1:2306551/148‑1 (MQ=255)
                                                                                                                   cTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATAtcatca                             >  2:18354/1‑148 (MQ=255)
                                                                                                                   cTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATAtcatca                             <  2:2680243/148‑1 (MQ=255)
                                                                                                                   cTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATAtcatca                             >  2:18252/1‑148 (MQ=255)
                                                                                                                        caATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAg                        >  2:928191/1‑148 (MQ=255)
                                                                                                                        caATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAg                        >  2:1018981/1‑148 (MQ=255)
                                                                                                                        caATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCATGGGCAAATGCCGTGTAATATCATCAAGCAg                        >  2:278960/1‑148 (MQ=255)
                                                                                                                         aaTGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAgg                       >  1:2459382/1‑148 (MQ=255)
                                                                                                                         aaTGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAgg                       >  1:2459485/1‑148 (MQ=255)
                                                                                                                          ttGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAggg                      <  2:785035/147‑1 (MQ=255)
                                                                                                                           tGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAgggg                     >  1:3457477/1‑148 (MQ=255)
                                                                                                                              tataGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTGGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATa                  <  2:3138768/148‑1 (MQ=255)
                                                                                                                              tataGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATa                  >  2:133046/1‑148 (MQ=255)
                                                                                                                                  gagaCTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTg              >  2:364825/1‑148 (MQ=255)
                                                                                                                                    gaCTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGcc            <  2:2831545/148‑1 (MQ=255)
                                                                                                                                    gaCTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGcc            >  2:1238556/1‑148 (MQ=255)
                                                                                                                                    gaCTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGcc            >  2:883166/1‑148 (MQ=255)
                                                                                                                                     aCTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGccc           <  2:840111/148‑1 (MQ=255)
                                                                                                                                     aCTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGccc           <  2:841282/148‑1 (MQ=255)
                                                                                                                                      cttaTCAGGTCTATATTTATCGCGGAATGGATATAACGACAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGcccc          <  1:2431359/145‑1 (MQ=255)
                                                                                                                                      cTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGcccc          >  2:2455562/1‑148 (MQ=255)
                                                                                                                                       tGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGccccc         >  2:3408694/1‑148 (MQ=255)
                                                                                                                                        gATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCt        >  2:2552236/1‑148 (MQ=255)
                                                                                                                                         aTCAGGTCTCTATTGATATCGGCCAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTg       >  1:435847/1‑148 (MQ=255)
                                                                                                                                         aTCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTg       >  2:2765507/1‑148 (MQ=255)
                                                                                                                                         aTCAGGGCCCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAACAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTg       <  2:2911921/148‑1 (MQ=255)
                                                                                                                                          tCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGc      >  2:1594306/1‑148 (MQ=255)
                                                                                                                                          tCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGc      >  2:102890/1‑148 (MQ=255)
                                                                                                                                           cAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGcc     >  1:2140679/1‑148 (MQ=255)
                                                                                                                                             ggTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCa   >  1:116970/1‑148 (MQ=255)
                                                                                                                                              gTCTCTATTGATATCGGCAAGTATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCAc  <  2:173982/148‑1 (MQ=255)
                                                                                                                                              gTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCAc  >  2:538630/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCAC  >  CP009273/3508957‑3509246

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCA  >  CP009273/3508955‑3509250
                                                                                                                                                   |                                                                                                                                                    
ACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCT                                                                                                                                                    <  SRR19973067.622812/150‑1 (MQ=60)
 CGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTA                                                                                                                                                   <  SRR19973067.3257866/150‑1 (MQ=60)
 CGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTA                                                                                                                                                   <  SRR19973067.3258029/150‑1 (MQ=60)
  GGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTAT                                                                                                                                                  <  SRR19973067.2327870/150‑1 (MQ=60)
  GGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTAT                                                                                                                                                  >  SRR19973067.2242801/1‑150 (MQ=60)
  GGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTAT                                                                                                                                                  <  SRR19973067.2327915/150‑1 (MQ=60)
    AACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTG                                                                                                                                                <  SRR19973067.3431715/150‑1 (MQ=60)
      CGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGAT                                                                                                                                              <  SRR19973067.598697/150‑1 (MQ=60)
      CGGATCCTGGCTTTTGTTGAACTCTTCGAGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGAT                                                                                                                                              >  SRR19973067.3377387/1‑150 (MQ=60)
       GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCCCAATGTTATAGAGACTGATCAGGTCTCTATTGATA                                                                                                                                             >  SRR19973067.418945/1‑150 (MQ=60)
       GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATA                                                                                                                                             >  SRR19973067.175384/1‑150 (MQ=60)
       GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATA                                                                                                                                             >  SRR19973067.175336/1‑150 (MQ=60)
       GGATCCTGGCTATTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATA                                                                                                                                             <  SRR19973067.1853639/150‑1 (MQ=60)
         ATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATC                                                                                                                                           <  SRR19973067.881273/150‑1 (MQ=60)
          TCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCG                                                                                                                                          <  SRR19973067.234631/150‑1 (MQ=60)
           CCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGG                                                                                                                                         >  SRR19973067.845809/1‑150 (MQ=60)
           CCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGG                                                                                                                                         >  SRR19973067.846987/1‑150 (MQ=60)
           CCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGG                                                                                                                                         >  SRR19973067.846989/1‑150 (MQ=60)
             TGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCA                                                                                                                                       >  SRR19973067.3158700/1‑150 (MQ=60)
             TGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCA                                                                                                                                       >  SRR19973067.3158513/1‑150 (MQ=60)
                        GAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACG                                                                                                                            <  SRR19973067.963777/150‑1 (MQ=60)
                          ACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAG                                                                                                                          <  SRR19973067.1488623/150‑1 (MQ=60)
                          ACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAG                                                                                                                          <  SRR19973067.1488611/150‑1 (MQ=60)
                          ACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAG                                                                                                                          >  SRR19973067.2208537/1‑150 (MQ=60)
                              TTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAA                                                                                                                      <  SRR19973067.2954142/150‑1 (MQ=60)
                                CGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGT                                                                                                                    <  SRR19973067.418158/150‑1 (MQ=60)
                                 GCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTC                                                                                                                   <  SRR19973067.1327447/150‑1 (MQ=60)
                                     AACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTG                                                                                                               <  SRR19973067.1017445/150‑1 (MQ=60)
                                       CTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTT                                                                                                             <  SRR19973067.266876/150‑1 (MQ=60)
                                          ATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAG                                                                                                          >  SRR19973067.2720795/1‑150 (MQ=60)
                                          ATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAG                                                                                                          >  SRR19973067.2584974/1‑150 (MQ=60)
                                               AAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACA                                                                                                     >  SRR19973067.147193/1‑150 (MQ=60)
                                                     CGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACC                                                                                               <  SRR19973067.1472886/150‑1 (MQ=60)
                                                      GGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCC                                                                                              >  SRR19973067.518571/1‑150 (MQ=60)
                                                         GCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGC                                                                                           >  SRR19973067.203029/1‑150 (MQ=60)
                                                         GCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGC                                                                                           >  SRR19973067.2011861/1‑150 (MQ=60)
                                                          CAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCA                                                                                          >  SRR19973067.661333/1‑150 (MQ=60)
                                                           AATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCAC                                                                                         >  SRR19973067.1085155/1‑150 (MQ=60)
                                                               TGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTAC                                                                                     >  SRR19973067.507523/1‑150 (MQ=60)
                                                                       GAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCT                                                                             >  SRR19973067.1994354/1‑150 (MQ=60)
                                                                        AACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTC                                                                            >  SRR19973067.2791035/1‑150 (MQ=60)
                                                                          CAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAA                                                                          >  SRR19973067.475265/1‑150 (MQ=60)
                                                                               AGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGA                                                                     >  SRR19973067.2324660/1‑150 (MQ=60)
                                                                                 CTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACA                                                                   >  SRR19973067.327546/1‑150 (MQ=60)
                                                                                      GTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGC                                                              >  SRR19973067.1292850/1‑150 (MQ=60)
                                                                                                TACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAA                                                    >  SRR19973067.702311/1‑150 (MQ=60)
                                                                                                TACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAA                                                    >  SRR19973067.702258/1‑150 (MQ=60)
                                                                                                   TTACTCATCAGTACGAATCTACACAATGTTATAGAGACTGACCTGGTCTATATTTATAGCGGCAAGGATATAACGACAAAAGTGGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATG                                                 <  SRR19973067.892435/150‑1 (MQ=60)
                                                                                                   GTACTCATCATAACGAATCTTCACAATTTTATAGACACTTATCAGGTATATATTTATATCGGCAAGTATATAAGGATAAAAGTGGGTGTTGACAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATG                                                 <  SRR19973067.892555/150‑1 (MQ=60)
                                                                                                     ACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCC                                               >  SRR19973067.2254003/1‑150 (MQ=60)
                                                                                                       TCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGT                                             >  SRR19973067.3568971/1‑150 (MQ=60)
                                                                                                              TACGCATCTACAAAATGTTATAGAGACTGAACAGGTCTATATTGATATCGCACAGGATATAACGATAAACGTGGGTGTTGACAACCCCCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATAT                                      <  SRR19973067.819884/150‑1 (MQ=60)
                                                                                                               ACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATC                                     >  SRR19973067.593571/1‑150 (MQ=60)
                                                                                                               ACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATC                                     >  SRR19973067.2725288/1‑150 (MQ=60)
                                                                                                                   ATCTCCACAATGTTATAGATACTTATCAGGGCGATTTTTATATCGGCACGGATATAACGAGAAATGTCGGTTTTGACAAACCCCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCA                                 <  SRR19973067.2321745/150‑1 (MQ=60)
                                                                                                                        CACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAG                            >  SRR19973067.419731/1‑150 (MQ=60)
                                                                                                                         ACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGG                           >  SRR19973067.2475550/1‑150 (MQ=60)
                                                                                                                         ACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGG                           >  SRR19973067.2475571/1‑150 (MQ=60)
                                                                                                                         ACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGG                           >  SRR19973067.2475653/1‑150 (MQ=60)
                                                                                                                           AATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGG                         >  SRR19973067.3479520/1‑150 (MQ=60)
                                                                                                                           AATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGG                         >  SRR19973067.3479461/1‑150 (MQ=60)
                                                                                                                           AATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGG                         >  SRR19973067.310414/1‑150 (MQ=60)
                                                                                                                                   AGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGC                 >  SRR19973067.203103/1‑150 (MQ=60)
                                                                                                                                        CTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAATGTGGGTGTTGACAAACCCCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCT            <  SRR19973067.1143484/150‑1 (MQ=60)
                                                                                                                                         TGATCAGGTCTCTATTGATATCGGCCAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTG           >  SRR19973067.438877/1‑150 (MQ=60)
                                                                                                                                           ATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCC         >  SRR19973067.2154989/1‑150 (MQ=60)
                                                                                                                                           ATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCC         >  SRR19973067.2888160/1‑150 (MQ=60)
                                                                                                                                             CAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCA       >  SRR19973067.118087/1‑150 (MQ=60)
                                                                                                                                             CAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCA       >  SRR19973067.2268920/1‑150 (MQ=60)
                                                                                                                                               GGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACT     >  SRR19973067.1791341/1‑150 (MQ=60)
                                                                                                                                               GGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACT     >  SRR19973067.2830434/1‑150 (MQ=60)
                                                                                                                                                  CTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTcgatcggaagagcacacg  >  SRR19973067.3127551/1‑132 (MQ=60)
                                                                                                                                                  CTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTagatcggaagagcacacg  >  SRR19973067.125103/1‑132 (MQ=60)
                                                                                                                                                   |                                                                                                                                                    
ACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCA  >  CP009273/3508955‑3509250

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