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 279,800 G→T 86.3% E297* (GAA→TAA)  yagF → CP4‑6 prophage; dehydratase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273279,8000GT86.3% 231.2 / 36.8 102E297* (GAA→TAA) yagFCP4‑6 prophage; dehydratase family protein
Reads supporting (aligned to +/- strand):  ref base G (4/10);  new base T (44/44);  total (48/54)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 102 total aligned reads displayed.

TGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGT  >  CP009273/279657‑279942
                                                                                                                                               |                                                                                                                                              
tGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTaaaa                                                                                                                                            <  1:3094664/148‑1 (MQ=255)
  ggTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAAcg                                                                                                                                          <  2:1602551/148‑1 (MQ=255)
    tCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAAcgcg                                                                                                                                        <  2:1304813/148‑1 (MQ=255)
     cTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGa                                                                                                                                       <  2:820890/148‑1 (MQ=255)
     cTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGa                                                                                                                                       <  2:2837160/148‑1 (MQ=255)
        gcgcTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATgg                                                                                                                                    >  2:1149318/1‑148 (MQ=255)
          gcTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTg                                                                                                                                  <  1:2872085/148‑1 (MQ=255)
           cTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGa                                                                                                                                 <  2:1716000/148‑1 (MQ=255)
           cTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGa                                                                                                                                 <  1:1583119/148‑1 (MQ=255)
           cTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGa                                                                                                                                 <  1:3142343/148‑1 (MQ=255)
             gCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGTCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATc                                                                                                                               <  2:2670840/148‑1 (MQ=255)
             gCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATc                                                                                                                               <  2:2670836/148‑1 (MQ=255)
              ccGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATcc                                                                                                                              <  2:1866559/148‑1 (MQ=255)
               cGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCTATCTAAAACGCGATGGTGATCCa                                                                                                                             <  2:1385377/148‑1 (MQ=255)
                   cTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCAcgcg                                                                                                                         >  2:438653/1‑148 (MQ=255)
                    tCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACgcgg                                                                                                                        >  2:1224248/1‑148 (MQ=255)
                    tCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACgcgg                                                                                                                        >  1:2970160/1‑148 (MQ=255)
                    tCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACgcgg                                                                                                                        >  1:2866376/1‑148 (MQ=255)
                     ccGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACgcggc                                                                                                                       >  1:2009079/1‑148 (MQ=255)
                     ccGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACgcggc                                                                                                                       >  1:2008987/1‑148 (MQ=255)
                      cgcgCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACgcggcg                                                                                                                      <  2:1987840/148‑1 (MQ=255)
                       gcgcTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGt                                                                                                                     <  2:1379053/148‑1 (MQ=255)
                        cgcTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGtt                                                                                                                    <  2:1526682/148‑1 (MQ=255)
                            ggCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTcggc                                                                                                                >  1:1582074/1‑148 (MQ=255)
                              cgcCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTcggcgg                                                                                                              <  2:2356255/148‑1 (MQ=255)
                              cgcCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTcggcgg                                                                                                              <  2:2356635/148‑1 (MQ=255)
                              cgcCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTcggcgg                                                                                                              <  1:2352651/148‑1 (MQ=255)
                                           ggCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTa                                                                                                 >  1:2100596/1‑148 (MQ=255)
                                              ggTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTActg                                                                                              <  1:2449390/148‑1 (MQ=255)
                                              ggTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTActg                                                                                              <  2:1637847/148‑1 (MQ=255)
                                                     cTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGcaca                                                                                       <  1:1805971/148‑1 (MQ=255)
                                                        gagaTCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTc                                                                                    <  1:2635568/148‑1 (MQ=255)
                                                           aTCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCgg                                                                                 >  1:1970327/1‑148 (MQ=255)
                                                           aTCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCgg                                                                                 >  1:1969843/1‑148 (MQ=255)
                                                                 cgccAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCg                                                                           >  2:157254/1‑148 (MQ=255)
                                                                 cgccAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCg                                                                           >  2:157189/1‑148 (MQ=255)
                                                                  gccAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGc                                                                          <  1:2936242/148‑1 (MQ=255)
                                                                      gTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCAc                                                                      >  1:2036098/1‑148 (MQ=255)
                                                                       tCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCAcg                                                                     <  2:1689833/148‑1 (MQ=255)
                                                                       tCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCAcg                                                                     <  2:1689184/148‑1 (MQ=255)
                                                                            gcgcgcgGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACgcggcg                                                                <  2:2099245/148‑1 (MQ=255)
                                                                            gcgcgcgGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACgcggcg                                                                <  2:2100182/148‑1 (MQ=255)
                                                                            gcgcgcgGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACgcggcg                                                                <  2:2101490/148‑1 (MQ=255)
                                                                              gcgcgGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCggg                                                              <  1:1118897/148‑1 (MQ=255)
                                                                               cgcgGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGc                                                             <  1:2756625/148‑1 (MQ=255)
                                                                                 cgGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTg                                                           >  2:819971/1‑148 (MQ=255)
                                                                                 cgGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTg                                                           <  2:1658407/148‑1 (MQ=255)
                                                                                   gTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCa                                                         >  1:431183/1‑148 (MQ=255)
                                                                                      aGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGa                                                      <  1:2507861/148‑1 (MQ=255)
                                                                                      aGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGa                                                      <  1:3137279/148‑1 (MQ=255)
                                                                                        cGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCAGGCCATCGCCCACGCGGCGGGCTGCACGATc                                                    >  2:2311030/1‑148 (MQ=255)
                                                                                          aGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATccc                                                  >  2:912624/1‑148 (MQ=255)
                                                                                                 tAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGtt                                           >  2:1153016/1‑148 (MQ=255)
                                                                                                  aGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTg                                          <  1:1018405/148‑1 (MQ=255)
                                                                                                  aGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTg                                          <  1:2547674/148‑1 (MQ=255)
                                                                                                  aGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTg                                          <  1:1018402/148‑1 (MQ=255)
                                                                                                  aGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTg                                          <  1:1018388/148‑1 (MQ=255)
                                                                                                   gCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGa                                         >  2:497781/1‑148 (MQ=255)
                                                                                                    ccGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAg                                        <  1:3210230/148‑1 (MQ=255)
                                                                                                    ccGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAg                                        <  1:438653/148‑1 (MQ=255)
                                                                                                     cgcgGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGc                                       <  2:1666754/148‑1 (MQ=255)
                                                                                                         gCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTg                                   >  2:485566/1‑148 (MQ=255)
                                                                                                         gCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTg                                   >  2:485546/1‑148 (MQ=255)
                                                                                                            tCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAc                                >  1:3234055/1‑148 (MQ=255)
                                                                                                            tCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAc                                >  1:2399170/1‑148 (MQ=255)
                                                                                                             caccacGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAcg                               >  2:806803/1‑148 (MQ=255)
                                                                                                             caccacGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAcg                               >  2:2712940/1‑148 (MQ=255)
                                                                                                               ccacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAcgcg                             >  1:503742/1‑148 (MQ=255)
                                                                                                               ccacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAcgcg                             <  2:1657645/148‑1 (MQ=255)
                                                                                                               ccacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAcgcg                             >  1:503934/1‑148 (MQ=255)
                                                                                                               ccacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAcgcg                             <  2:3295176/148‑1 (MQ=255)
                                                                                                                cacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACgcgc                            >  1:2934967/1‑148 (MQ=255)
                                                                                                                cacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACgcgc                            >  1:2934960/1‑148 (MQ=255)
                                                                                                                cacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACgcgc                            >  1:2934849/1‑148 (MQ=255)
                                                                                                                 acGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCa                           <  1:486844/148‑1 (MQ=255)
                                                                                                                  cgcgGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCAt                          <  1:2529584/148‑1 (MQ=255)
                                                                                                                    cgGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCa                        >  2:1415131/1‑148 (MQ=255)
                                                                                                                      ggATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTGACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAAc                      <  1:1589061/148‑1 (MQ=255)
                                                                                                                          atCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTa                  >  2:2588082/1‑148 (MQ=255)
                                                                                                                          atCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTa                  <  2:1499312/148‑1 (MQ=255)
                                                                                                                          atCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTa                  <  2:1499169/148‑1 (MQ=255)
                                                                                                                           tCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTaa                 <  1:550625/148‑1 (MQ=255)
                                                                                                                           tCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTaa                 <  1:1925397/148‑1 (MQ=255)
                                                                                                                           tCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTaa                 >  1:2989664/1‑148 (MQ=255)
                                                                                                                           tCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTaa                 <  1:2631061/148‑1 (MQ=255)
                                                                                                                           tCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTaa                 <  2:3075994/148‑1 (MQ=255)
                                                                                                                            ccTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCCAGTTACTGCTGCACATTCCGGCCATAGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTaaa                >  1:2653051/1‑148 (MQ=255)
                                                                                                                            ccTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTaaa                >  1:1875963/1‑148 (MQ=255)
                                                                                                                              tctcCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGt              >  1:3146620/1‑148 (MQ=255)
                                                                                                                              tctcCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGt              >  1:1258014/1‑148 (MQ=255)
                                                                                                                               ctcCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTg             <  2:2808343/148‑1 (MQ=255)
                                                                                                                               ctcCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTg             >  2:2227609/1‑148 (MQ=255)
                                                                                                                                tcCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGc            >  2:332385/1‑148 (MQ=255)
                                                                                                                                tcCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGc            >  1:1866238/1‑148 (MQ=255)
                                                                                                                                    aTAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCcgcg        >  2:3183986/1‑148 (MQ=255)
                                                                                                                                     tAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGt       <  2:2202076/148‑1 (MQ=255)
                                                                                                                                      aaaGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTc      >  2:2514478/1‑148 (MQ=255)
                                                                                                                                       aaGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCt     >  1:524233/1‑148 (MQ=255)
                                                                                                                                         gCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTgg   >  1:2579148/1‑148 (MQ=255)
                                                                                                                                          ccATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGt  >  2:1749722/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                              
TGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGT  >  CP009273/279657‑279942

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
AGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGAGCGTG  >  CP009273/279651‑279949
                                                                                                                                                     |                                                                                                                                                     
AGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCT                                                                                                                                                       <  SRR19973070.2371405/150‑1 (MQ=60)
AGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCT                                                                                                                                                       <  SRR19973070.2900292/150‑1 (MQ=60)
      TGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACG                                                                                                                                                 <  SRR19973070.3114600/150‑1 (MQ=60)
                GCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGAT                                                                                                                                       <  SRR19973070.2890690/150‑1 (MQ=60)
                 CTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATC                                                                                                                                      <  SRR19973070.1594211/150‑1 (MQ=60)
                 CTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATC                                                                                                                                      <  SRR19973070.3162648/150‑1 (MQ=60)
                        ACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGG                                                                                                                               >  SRR19973070.2989240/1‑150 (MQ=60)
                        ACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGG                                                                                                                               >  SRR19973070.669847/1‑150 (MQ=60)
                        ACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGG                                                                                                                               >  SRR19973070.2884946/1‑150 (MQ=60)
                         CTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGC                                                                                                                              >  SRR19973070.2022587/1‑150 (MQ=60)
                         CTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGC                                                                                                                              >  SRR19973070.2022679/1‑150 (MQ=60)
                                CTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGC                                                                                                                       >  SRR19973070.1593160/1‑150 (MQ=60)
                                    CGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCT                                                                                                                   <  SRR19973070.2368242/150‑1 (MQ=60)
                                               CAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTA                                                                                                        >  SRR19973070.2114692/1‑150 (MQ=60)
                                                    GGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCT                                                                                                   <  SRR19973070.2465411/150‑1 (MQ=60)
                                                           CTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATT                                                                                            <  SRR19973070.1818437/150‑1 (MQ=60)
                                                              GAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCG                                                                                         <  SRR19973070.2652753/150‑1 (MQ=60)
                                                               AGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGG                                                                                        >  SRR19973070.1983152/1‑150 (MQ=60)
                                                               AGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGG                                                                                        >  SRR19973070.1983639/1‑150 (MQ=60)
                                                                        GCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCC                                                                               <  SRR19973070.2955170/150‑1 (MQ=60)
                                                                          CAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCAC                                                                             >  SRR19973070.2049872/1‑150 (MQ=60)
                                                                                    GCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCT                                                                   <  SRR19973070.1126538/150‑1 (MQ=60)
                                                                                     CGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTG                                                                  <  SRR19973070.2774314/150‑1 (MQ=60)
                                                                                       CGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCA                                                                >  SRR19973070.434318/1‑150 (MQ=60)
                                                                                            AGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATC                                                           <  SRR19973070.3157549/150‑1 (MQ=60)
                                                                                            AGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATC                                                           <  SRR19973070.2524164/150‑1 (MQ=60)
                                                                                                        AGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAG                                               <  SRR19973070.1025635/150‑1 (MQ=60)
                                                                                                        AGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAG                                               <  SRR19973070.1025621/150‑1 (MQ=60)
                                                                                                        AGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAG                                               <  SRR19973070.1025638/150‑1 (MQ=60)
                                                                                                        AGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAG                                               <  SRR19973070.2564336/150‑1 (MQ=60)
                                                                                                          CCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCA                                             <  SRR19973070.3230994/150‑1 (MQ=60)
                                                                                                          CCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCA                                             <  SRR19973070.441838/150‑1 (MQ=60)
                                                                                                                CATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAC                                       >  SRR19973070.3254930/1‑150 (MQ=60)
                                                                                                                CATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAC                                       >  SRR19973070.2415001/1‑150 (MQ=60)
                                                                                                                   CACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCG                                    >  SRR19973070.507427/1‑150 (MQ=60)
                                                                                                                   CACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCG                                    >  SRR19973070.507234/1‑150 (MQ=60)
                                                                                                                    ACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGC                                   >  SRR19973070.2953774/1‑150 (MQ=60)
                                                                                                                    ACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGC                                   >  SRR19973070.2953886/1‑150 (MQ=60)
                                                                                                                    ACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGC                                   >  SRR19973070.2953893/1‑150 (MQ=60)
                                                                                                                       ACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATC                                <  SRR19973070.490267/150‑1 (MQ=60)
                                                                                                                        CGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCA                               <  SRR19973070.2546109/150‑1 (MQ=60)
                                                                                                                            GGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTGACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCG                           <  SRR19973070.1600177/150‑1 (MQ=60)
                                                                                                                               TATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAA                        >  SRR19973070.3008818/1‑150 (MQ=60)
                                                                                                                                ATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCCAGTTACTGCTGCACATTCCGGCCATAGCCCACGCGCCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAAACGTAAA                       >  SRR19973070.2670313/1‑150 (MQ=60)
                                                                                                                                ATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAA                       >  SRR19973070.1888741/1‑150 (MQ=60)
                                                                                                                                 TCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAG                      <  SRR19973070.554279/150‑1 (MQ=60)
                                                                                                                                 TCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAG                      <  SRR19973070.2648220/150‑1 (MQ=60)
                                                                                                                                 TCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAG                      <  SRR19973070.1938372/150‑1 (MQ=60)
                                                                                                                                  CCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGT                     >  SRR19973070.3166954/1‑150 (MQ=60)
                                                                                                                                  CCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGT                     >  SRR19973070.1266831/1‑150 (MQ=60)
                                                                                                                                    TCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGC                   >  SRR19973070.1878985/1‑150 (MQ=60)
                                                                                                                                           TAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCT            >  SRR19973070.527788/1‑150 (MQ=60)
                                                                                                                                             AAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGG          >  SRR19973070.2596008/1‑150 (MQ=60)
                                                                                                                                                  ATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGAGC     >  SRR19973070.2967484/1‑150 (MQ=60)
                                                                                                                                                  ATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGAGC     >  SRR19973070.2077155/1‑150 (MQ=60)
                                                                                                                                                  ATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGAGC     >  SRR19973070.964431/1‑150 (MQ=60)
                                                                                                                                                  ATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGAGC     >  SRR19973070.964551/1‑150 (MQ=60)
                                                                                                                                                  ANCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGAGC     >  SRR19973070.2967762/1‑150 (MQ=60)
                                                                                                                                                   TCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGAGCG    <  SRR19973070.2512879/150‑1 (MQ=60)
                                                                                                                                                    CTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGAGCGT   >  SRR19973070.2208338/1‑150 (MQ=60)
                                                                                                                                                    CTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGAGCGT   >  SRR19973070.2208271/1‑150 (MQ=60)
                                                                                                                                                     TAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGAGCGTG  >  SRR19973070.2575524/1‑150 (MQ=60)
                                                                                                                                                     |                                                                                                                                                     
AGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGAGCGTG  >  CP009273/279651‑279949

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: