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 F309 I0 R1 1697 197.2 6627756 95.7% 6342762 147.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 217,887 C→T 91.8% W74* (TGG→TGA metI ← DL‑methionine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273217,8870CT91.8% 495.0 / 28.8 159W74* (TGG→TGAmetIDL‑methionine transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (3/10);  new base T (78/68);  total (81/78)
Fisher's exact test for biased strand distribution p-value = 4.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 162 total aligned reads displayed.

AAGCTTCAATTAACCCGGTTGGGATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACAATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAACCGGCAGACCAATCACAAAGCCAAAA  >  CP009273/217742‑218034
                                                                                                                                                 |                                                                                                                                                    
aaGCTTCAATTAACCCGGTTGGGATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCa                                                                                                                                                    >  2:2558473/1‑148 (MQ=255)
  gCTTCAATTAACCCGGTTGGGATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATa                                                                                                                                                  <  2:493087/148‑1 (MQ=255)
  gCTTCAATTAACCCGGTTGGGATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATa                                                                                                                                                  >  2:1218621/1‑148 (MQ=255)
    ttCAATTAACCCGGTTGGGATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATAca                                                                                                                                                <  1:2201313/148‑1 (MQ=255)
     tCAATTAACCCGGTTGGGATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATAcaa                                                                                                                                               >  2:101953/1‑148 (MQ=255)
      cAATTAACCCGGTTGGGATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATAcaag                                                                                                                                              <  1:1856394/148‑1 (MQ=255)
        aTTAACCCGGTTGGGATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTACACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATAcaagca                                                                                                                                            <  1:1571948/148‑1 (MQ=255)
          tAACCCGGTTGGGATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATAcaagcaag                                                                                                                                          <  2:1159988/148‑1 (MQ=255)
          tAACCCGGTTGGGATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATAcaagcaag                                                                                                                                          <  2:1592934/148‑1 (MQ=255)
          tAACCCGGTTGGGATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATAcaagcaag                                                                                                                                          <  2:2361129/148‑1 (MQ=255)
          tAACCCGGTTGGGATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATAcaagcaag                                                                                                                                          <  2:2420101/148‑1 (MQ=255)
             cccGGTTGGGATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATa                                                                                                                                       <  2:2455771/148‑1 (MQ=255)
                 gTTGGGATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGa                                                                                                                                   <  2:2022854/148‑1 (MQ=255)
                    gggATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACg                                                                                                                                <  2:860033/148‑1 (MQ=255)
                      gATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACggg                                                                                                                              <  1:1131099/148‑1 (MQ=255)
                       aTCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGa                                                                                                                             <  1:128550/148‑1 (MQ=255)
                        tctcCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGAt                                                                                                                            >  1:2397108/1‑148 (MQ=255)
                         ctcCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATg                                                                                                                           <  2:806790/148‑1 (MQ=255)
                         ctcCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATg                                                                                                                           <  2:806406/148‑1 (MQ=255)
                              gcagAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACg                                                                                                                      >  2:592325/1‑148 (MQ=255)
                              gcagAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACg                                                                                                                      >  2:1768549/1‑148 (MQ=255)
                                 gagCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGaa                                                                                                                   >  2:806575/1‑148 (MQ=255)
                                   gCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGaaaa                                                                                                                 <  2:2344909/148‑1 (MQ=255)
                                    cGTTCTCGACCATACCGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGACTCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAAt                                                                                                                >  1:2450790/1‑148 (MQ=255)
                                      ttCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGt                                                                                                              <  1:1765887/148‑1 (MQ=255)
                                       tctcGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGtt                                                                                                             <  2:1250597/148‑1 (MQ=255)
                                             ccATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCAt                                                                                                       <  1:2213257/148‑1 (MQ=255)
                                               aTACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCg                                                                                                     >  1:1517895/1‑148 (MQ=255)
                                                tACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGc                                                                                                    >  2:518965/1‑148 (MQ=255)
                                                   gggCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGa                                                                                                 >  1:1674388/1‑148 (MQ=255)
                                                    ggCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGaa                                                                                                >  2:1640148/1‑148 (MQ=255)
                                                    ggCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGaa                                                                                                <  1:674430/148‑1 (MQ=255)
                                                    ggCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACAATCGCAGaa                                                                                                <  2:1509463/148‑1 (MQ=255)
                                                         tAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGt                                                                                           >  1:897628/1‑148 (MQ=255)
                                                         tAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGt                                                                                           >  1:897635/1‑148 (MQ=255)
                                                         tAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGt                                                                                           <  1:779969/148‑1 (MQ=255)
                                                         tAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAACCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGt                                                                                           >  2:3287683/1‑148 (MQ=255)
                                                          aaaCGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTa                                                                                          <  1:2093066/148‑1 (MQ=255)
                                                          aaaCGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTa                                                                                          >  2:1505809/1‑148 (MQ=255)
                                                          aaaCGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTa                                                                                          <  1:2093058/148‑1 (MQ=255)
                                                             cGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGa                                                                                       >  1:620193/1‑148 (MQ=255)
                                                             cGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGa                                                                                       >  1:1530343/1‑148 (MQ=255)
                                                              ggTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGAt                                                                                      <  1:1060033/148‑1 (MQ=255)
                                                                tgctgcACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATAc                                                                                    <  2:1486901/148‑1 (MQ=255)
                                                                tgctgcACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATAc                                                                                    >  1:2803277/1‑148 (MQ=255)
                                                                tgctgcACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATAc                                                                                    <  2:2938105/148‑1 (MQ=255)
                                                                tgctgcACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACAATCGCAGAAACGGTACGATAc                                                                                    <  1:2807374/148‑1 (MQ=255)
                                                                   tgcACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACAATCGCAGAAACGGTACGATACAGc                                                                                 <  1:1999999/148‑1 (MQ=255)
                                                                     cACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCtt                                                                               <  1:489136/148‑1 (MQ=255)
                                                                     cACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCtt                                                                               >  1:2871814/1‑148 (MQ=255)
                                                                       ccAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTcg                                                                             >  1:2370903/1‑148 (MQ=255)
                                                                       ccAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTcg                                                                             <  1:69010/148‑1 (MQ=255)
                                                                        cAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTcgc                                                                            >  1:2533770/1‑148 (MQ=255)
                                                                        cAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTcgc                                                                            >  2:1124758/1‑148 (MQ=255)
                                                                          aCGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGt                                                                          >  2:140512/1‑148 (MQ=255)
                                                                          aCGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGt                                                                          >  2:136876/1‑148 (MQ=255)
                                                                            ggTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTa                                                                        <  2:2077544/148‑1 (MQ=255)
                                                                               tAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCa                                                                    >  1:996682/1‑148 (MQ=255)
                                                                                 aCGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCaat                                                                  >  1:1484987/1‑148 (MQ=255)
                                                                                 aCGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCaat                                                                  >  2:1375203/1‑148 (MQ=255)
                                                                                  cGGAACAATCGCTGCCTGCAAACCCATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCaata                                                                 <  2:1200434/148‑1 (MQ=255)
                                                                                  cGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCaata                                                                 <  1:2117929/148‑1 (MQ=255)
                                                                                   ggAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCaataa                                                                <  1:1141062/148‑1 (MQ=255)
                                                                                    gAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCaataat                                                               >  2:2848047/1‑148 (MQ=255)
                                                                                     aaCAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAAtt                                                              >  2:2935034/1‑148 (MQ=255)
                                                                                     aaCAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAAtt                                                              >  2:575996/1‑148 (MQ=255)
                                                                                      aCAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAAATAAtt                                                              >  1:1215908/1‑148 (MQ=255)
                                                                                      aCAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAAttt                                                             >  1:2463086/1‑148 (MQ=255)
                                                                                       cAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCACGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTg                                                            >  2:3192219/1‑148 (MQ=255)
                                                                                             ctgcctgcAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCgg                                                      <  1:2558473/148‑1 (MQ=255)
                                                                                             ctgcctgcAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCgg                                                      >  2:2362817/1‑148 (MQ=255)
                                                                                             ctgcctgcAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCgg                                                      >  2:2360822/1‑148 (MQ=255)
                                                                                               gcctgcAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGAc                                                    >  1:2875319/1‑148 (MQ=255)
                                                                                                   gcAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGt                                                >  1:1303450/1‑148 (MQ=255)
                                                                                                     aaaCCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGa                                              <  1:1094065/148‑1 (MQ=255)
                                                                                                       aCCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACa                                            <  2:2728881/148‑1 (MQ=255)
                                                                                                        ccAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAATCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACAt                                           >  2:20898/1‑148 (MQ=255)
                                                                                                         cAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATa                                          <  2:3100009/148‑1 (MQ=255)
                                                                                                         cAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATa                                          <  2:3099649/148‑1 (MQ=255)
                                                                                                           aTCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATaaa                                        >  2:2004751/1‑148 (MQ=255)
                                                                                                            tCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAg                                       <  2:446860/148‑1 (MQ=255)
                                                                                                              gATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACAATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCa                                     <  1:252235/148‑1 (MQ=255)
                                                                                                                tGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGa                                   >  1:871404/1‑148 (MQ=255)
                                                                                                                 gTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGaa                                  <  1:1249073/148‑1 (MQ=255)
                                                                                                                     cGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGcc                              <  1:2732905/148‑1 (MQ=255)
                                                                                                                      gACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCa                             >  1:3210379/1‑148 (MQ=255)
                                                                                                                       aCAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCaa                            >  2:2284263/1‑148 (MQ=255)
                                                                                                                       aCAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCaa                            >  1:64518/1‑148 (MQ=255)
                                                                                                                        cAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAAc                           >  1:1234388/1‑148 (MQ=255)
                                                                                                                            aaCGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAACCGGc                       >  1:3155823/1‑148 (MQ=255)
                                                                                                                              cgcgGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAACCGGCAg                     >  1:807615/1‑148 (MQ=255)
                                                                                                                                    tAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAACCGGCAGACCAAt               <  1:2787832/148‑1 (MQ=255)
                                                                                                                                     aaaCGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAACCGGCAGACCAATc              <  1:1354271/148‑1 (MQ=255)
                                                                                                                                      aaCGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAACCGGCAGACCAATca             <  1:3287683/148‑1 (MQ=255)
                                                                                                                                       aCGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAACCGGCAGACCAATcac            >  2:192765/1‑148 (MQ=255)
                                                                                                                                          gAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAAGAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAACCGGCAGACCAATCACaaa         >  1:2704865/1‑148 (MQ=255)
                                                                                                                                           aaTCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAACCGGCAGACCAATCACAAAg        >  2:3196621/1‑148 (MQ=255)
                                                                                                                                           aaTCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAACCGGCAGACCAATCACAAAg        >  2:2382326/1‑148 (MQ=255)
                                                                                                                                           aaTCATCCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACAATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAACCGGCAGACCAATCACAAAg        <  1:2339570/148‑1 (MQ=255)
                                                                                                                                            atcattCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAACCGGCAGACCAATCACAAAGc       >  1:3102503/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
AAGCTTCAATTAACCCGGTTGGGATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACAATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAACCGGCAGACCAATCACAAAGCCAAAA  >  CP009273/217742‑218034

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
AGCTTCAATTAACCCGGTTGGGATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACAATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAACCGGCAGACCAATCACAAAGCCAAAA  >  CP009273/217743‑218034
                                                                                                                                                |                                                                                                                                                    
AGCTTCAATTAACCCGGTTGGGATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATAC                                                                                                                                                 >  SRR19973070.2816264/1‑150 (MQ=60)
   TTCAATTAACCCGGTTGGGATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAG                                                                                                                                              <  SRR19973070.2215823/150‑1 (MQ=60)
     CAATTAACCCGGTTGGGATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCA                                                                                                                                            <  SRR19973070.1869107/150‑1 (MQ=60)
       ATTAACCCGGTTGGGATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTACACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAG                                                                                                                                          <  SRR19973070.1582988/150‑1 (MQ=60)
         TAACCCGGTTGGGATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGAT                                                                                                                                        >  SRR19973070.1381153/1‑150 (MQ=60)
                 TTGGGATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACG                                                                                                                                <  SRR19973070.1561232/150‑1 (MQ=60)
                     GATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTTCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGAT                                                                                                                            >  SRR19973070.878389/1‑150 (MQ=60)
                     GATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGAT                                                                                                                            >  SRR19973070.2412929/1‑150 (MQ=60)
                     GATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGAT                                                                                                                            <  SRR19973070.1138905/150‑1 (MQ=60)
                      ATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATG                                                                                                                           <  SRR19973070.129682/150‑1 (MQ=60)
                      ATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATG                                                                                                                           <  SRR19973070.129673/150‑1 (MQ=60)
                         TCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAA                                                                                                                        <  SRR19973070.958899/150‑1 (MQ=60)
                          CCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATGAACGGGATGGAAC                                                                                                                       <  SRR19973070.1227161/150‑1 (MQ=60)
                                GAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAA                                                                                                                 >  SRR19973070.2988888/1‑150 (MQ=60)
                                 AGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAAT                                                                                                                >  SRR19973070.1539791/1‑150 (MQ=60)
                                 AGCGTTCTCGACCATACCGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGACTCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAAT                                                                                                                >  SRR19973070.2466826/1‑150 (MQ=60)
                                  GCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATG                                                                                                               <  SRR19973070.1887978/150‑1 (MQ=60)
                                     TTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTC                                                                                                            <  SRR19973070.1778102/150‑1 (MQ=60)
                                         CGACCATACGCGCAATAAACGGAGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTCCACCA                                                                                                        <  SRR19973070.1749369/150‑1 (MQ=60)
                                            CCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCG                                                                                                     >  SRR19973070.1528533/1‑150 (MQ=60)
                                            CCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCG                                                                                                     <  SRR19973070.2227804/150‑1 (MQ=60)
                                                ACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGA                                                                                                 >  SRR19973070.1685856/1‑150 (MQ=60)
                                                  GGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAA                                                                                               >  SRR19973070.2129398/1‑150 (MQ=60)
                                                   GGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAAC                                                                                              <  SRR19973070.679238/150‑1 (MQ=60)
                                                     CAATAAACGGTGCTGCACCAACAGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGG                                                                                            >  SRR19973070.2851262/1‑150 (MQ=60)
                                                      AATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGT                                                                                           >  SRR19973070.903986/1‑150 (MQ=60)
                                                      AATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGT                                                                                           >  SRR19973070.903993/1‑150 (MQ=60)
                                                        TAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTAC                                                                                         <  SRR19973070.785418/150‑1 (MQ=60)
                                                         AAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACG                                                                                        <  SRR19973070.2106004/150‑1 (MQ=60)
                                                         AAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACG                                                                                        <  SRR19973070.2107130/150‑1 (MQ=60)
                                                         AAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACG                                                                                        <  SRR19973070.2107138/150‑1 (MQ=60)
                                                          AACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGA                                                                                       >  SRR19973070.1541092/1‑150 (MQ=60)
                                                          AACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGA                                                                                       >  SRR19973070.624558/1‑150 (MQ=60)
                                                          AACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGAAGTACCGACCAGAACGCGGGTAAACGGACTCATTCATACAAGCACGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGA                                                                                       >  SRR19973070.625257/1‑150 (MQ=60)
                                                             GGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATAC                                                                                    <  SRR19973070.1067467/150‑1 (MQ=60)
                                                             GGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATAC                                                                                    >  SRR19973070.2821390/1‑150 (MQ=60)
                                                               TGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACAATCGCAGAAACGGTACGATACAG                                                                                  <  SRR19973070.2825519/150‑1 (MQ=60)
                                                               TGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACAATCGCAGAAACGGTACGATACAG                                                                                  <  SRR19973070.2784239/150‑1 (MQ=60)
                                                               TGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACAATCGCAGAAACGGTACGATACAG                                                                                  <  SRR19973070.1166841/150‑1 (MQ=60)
                                                                  TGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTT                                                                               >  SRR19973070.2890417/1‑150 (MQ=60)
                                                                  TGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACAATCGCAGAAACGGTACGATACAGCTT                                                                               <  SRR19973070.2013544/150‑1 (MQ=60)
                                                                    CACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCG                                                                             >  SRR19973070.2386589/1‑150 (MQ=60)
                                                                    CACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCG                                                                             >  SRR19973070.2386591/1‑150 (MQ=60)
                                                                    CACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCG                                                                             <  SRR19973070.2166590/150‑1 (MQ=60)
                                                                    CACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCG                                                                             <  SRR19973070.492565/150‑1 (MQ=60)
                                                                     ACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGC                                                                            >  SRR19973070.2550330/1‑150 (MQ=60)
                                                                     ACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGC                                                                            >  SRR19973070.2550441/1‑150 (MQ=60)
                                                                      CCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCG                                                                           <  SRR19973070.69775/150‑1 (MQ=60)
                                                                            GTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCA                                                                    >  SRR19973070.1003787/1‑150 (MQ=60)
                                                                              TAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAAT                                                                  >  SRR19973070.1495383/1‑150 (MQ=60)
                                                                                 CGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAAT                                                               <  SRR19973070.2132123/150‑1 (MQ=60)
                                                                                  GGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATT                                                              <  SRR19973070.1148979/150‑1 (MQ=60)
                                                                                   GAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAAATAATT                                                              >  SRR19973070.1224430/1‑150 (MQ=60)
                                                                                   GAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTT                                                             >  SRR19973070.2479166/1‑150 (MQ=60)
                                                                                    AACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTG                                                            >  SRR19973070.1320662/1‑150 (MQ=60)
                                                                                            CTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGAC                                                    >  SRR19973070.2893942/1‑150 (MQ=60)
                                                                                            CTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGAC                                                    <  SRR19973070.2575188/150‑1 (MQ=60)
                                                                                                CTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGT                                                >  SRR19973070.1312511/1‑150 (MQ=60)
                                                                                                  GCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGA                                              <  SRR19973070.3254627/150‑1 (MQ=60)
                                                                                                    AAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACA                                            <  SRR19973070.1101615/150‑1 (MQ=60)
                                                                                                    AAACAAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGTATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACA                                            <  SRR19973070.1101619/150‑1 (MQ=60)
                                                                                                            CGATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAG                                    <  SRR19973070.2501680/150‑1 (MQ=60)
                                                                                                             GATGTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGA                                   >  SRR19973070.877508/1‑150 (MQ=60)
                                                                                                             GATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACAATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGA                                   <  SRR19973070.253951/150‑1 (MQ=60)
                                                                                                               TGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACAATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAAC                                 >  SRR19973070.81923/1‑150 (MQ=60)
                                                                                                                GTACCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACG                                <  SRR19973070.1257832/150‑1 (MQ=60)
                                                                                                                   CCGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCA                             >  SRR19973070.3231145/1‑150 (MQ=60)
                                                                                                                    CGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAA                            <  SRR19973070.2750521/150‑1 (MQ=60)
                                                                                                                    CGACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAA                            >  SRR19973070.65242/1‑150 (MQ=60)
                                                                                                                     GACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAAC                           >  SRR19973070.3311896/1‑150 (MQ=60)
                                                                                                                     GACAATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAAC                           >  SRR19973070.1243029/1‑150 (MQ=60)
                                                                                                                         ATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAACCGGC                       >  SRR19973070.3175689/1‑150 (MQ=60)
                                                                                                                         ATAACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAACCGGC                       >  SRR19973070.3176220/1‑150 (MQ=60)
                                                                                                                           AACGCGGGTAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAACCGGCAG                     >  SRR19973070.813156/1‑150 (MQ=60)
                                                                                                                                   TAAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAACCGGCAGACCAATCA             <  SRR19973070.2805780/150‑1 (MQ=60)
                                                                                                                                    AAACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAACCGGCAGACCAATCAC            <  SRR19973070.1363532/150‑1 (MQ=60)
                                                                                                                                     AACGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAACCGGCAGACCAATCACA           <  SRR19973070.3308785/150‑1 (MQ=60)
                                                                                                                                       CGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAACCGGCAGACCAATCACAAA         <  SRR19973070.695124/150‑1 (MQ=60)
                                                                                                                                       CGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAACCGGCAGACCAATCACAAA         <  SRR19973070.695006/150‑1 (MQ=60)
                                                                                                                                       CGGAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAAGAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAACCGGCAGACCAATCACAAA         >  SRR19973070.2722356/1‑150 (MQ=60)
                                                                                                                                         GAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAACCGGCAGACCAATCACAAAGC       >  SRR19973070.2416041/1‑150 (MQ=60)
                                                                                                                                         GAATCATTCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACCATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAACCGGCAGACCAATCACAAAGC       >  SRR19973070.3122509/1‑150 (MQ=60)
                                                                                                                                          AATCATCCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACAATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAACCGGCAGACCAATCACAAAGCC      <  SRR19973070.2355096/150‑1 (MQ=60)
                                                                                                                                          AATCATCCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACAATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAACCGGCAGACCAATCACAAAGCC      <  SRR19973070.2354940/150‑1 (MQ=60)
                                                                                                                                              ATCCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACAATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAACCGGCAGACCAATCACAAAGCCAAAA  >  SRR19973070.2624637/1‑150 (MQ=60)
                                                                                                                                                |                                                                                                                                                    
AGCTTCAATTAACCCGGTTGGGATCTCCAGCAGAGCGTTCTCGACCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATGAACGGGATGGAACGGAAAATGTTCACAATCGCAGAAACGGTACGATACAGCTTCGCGTTAGCAATAATTTGCCCCGGACGCGTGACATAAAGCAGAACGCCAACCGGCAGACCAATCACAAAGCCAAAA  >  CP009273/217743‑218034

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: