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 1,820,888 C→A 88.6% C98F (TGC→TTC)  astE ← succinylglutamate desuccinylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,820,8880CA88.6% 410.4 / 48.5 166C98F (TGC→TTC) astEsuccinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base C (9/10);  new base A (71/76);  total (80/86)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 170 total aligned reads displayed.

TAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTG  >  CP009273/1820743‑1821035
                                                                                                                                                 |                                                                                                                                                   
tAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGaaa                                                                                                                                                   <  2:1275860/148‑1 (MQ=255)
 aGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAAt                                                                                                                                                  <  2:2539733/148‑1 (MQ=255)
 aGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAAt                                                                                                                                                  <  2:1719655/148‑1 (MQ=255)
  gATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATa                                                                                                                                                 <  1:3006547/148‑1 (MQ=255)
    tCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAAc                                                                                                                                               <  1:1737778/148‑1 (MQ=255)
    tCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAAc                                                                                                                                               >  2:23467/1‑148 (MQ=255)
     cAATGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACg                                                                                                                                              >  1:1078695/1‑148 (MQ=255)
     cAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGANAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACg                                                                                                                                              <  1:1683242/148‑1 (MQ=255)
     cAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACGGTTCCTGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACg                                                                                                                                              <  1:1152266/148‑1 (MQ=255)
      aaGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGt                                                                                                                                             <  1:2309758/148‑1 (MQ=255)
      aaGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGt                                                                                                                                             <  1:951878/148‑1 (MQ=255)
      aaGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGt                                                                                                                                             <  1:1318668/148‑1 (MQ=255)
       aGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGtt                                                                                                                                            >  2:1376755/1‑148 (MQ=255)
       aGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGtt                                                                                                                                            >  1:3153219/1‑148 (MQ=255)
        ggTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGttt                                                                                                                                           >  1:1754112/1‑148 (MQ=255)
        ggTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGttt                                                                                                                                           >  1:1754113/1‑148 (MQ=255)
          tGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTcc                                                                                                                                         <  2:1917901/148‑1 (MQ=255)
           gCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTccc                                                                                                                                        <  2:419603/148‑1 (MQ=255)
              aGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTg                                                                                                                                     >  1:1497988/1‑148 (MQ=255)
              aGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTg                                                                                                                                     <  1:2463741/148‑1 (MQ=255)
               gcgcACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCACTAACGTTTCCCTTGc                                                                                                                                    >  2:2776639/1‑148 (MQ=255)
               gcgcACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGc                                                                                                                                    <  2:1400037/148‑1 (MQ=255)
                 gcACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCtt                                                                                                                                  <  2:699622/148‑1 (MQ=255)
                 gcACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCtt                                                                                                                                  <  2:699550/148‑1 (MQ=255)
                 gcACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCtt                                                                                                                                  <  2:2435173/148‑1 (MQ=255)
                  cacaGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTc                                                                                                                                 <  2:1718137/148‑1 (MQ=255)
                       aTTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAgcgc                                                                                                                            >  1:942940/1‑148 (MQ=255)
                       aTTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAgcgc                                                                                                                            >  1:1656230/1‑148 (MQ=255)
                        ttCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCa                                                                                                                           >  2:2798976/1‑148 (MQ=255)
                        ttCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCa                                                                                                                           <  2:466770/148‑1 (MQ=255)
                        ttCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCa                                                                                                                           >  1:2725453/1‑148 (MQ=255)
                        ttCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCa                                                                                                                           <  2:466557/148‑1 (MQ=255)
                        ttCTTTGCCCGGGTAATAAAAATCTTCCAGGCTCTGTTCCAGTTCGCGCGCCCTACAGGTTTCTCCGCTTTCAGCTATTAGCTGCCAACGACCGCCATACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCa                                                                                                                           <  2:466156/148‑1 (MQ=255)
                                ccTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGAtt                                                                                                                   <  1:632742/148‑1 (MQ=255)
                                ccTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGAtt                                                                                                                   <  2:203124/148‑1 (MQ=255)
                                 cTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTc                                                                                                                  <  1:347723/148‑1 (MQ=255)
                                  tGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTcc                                                                                                                 >  1:223060/1‑148 (MQ=255)
                                      cATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAg                                                                                                             >  2:2319033/1‑148 (MQ=255)
                                       aTAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAgg                                                                                                            <  2:2725068/148‑1 (MQ=255)
                                       aTAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAgg                                                                                                            <  1:1582328/148‑1 (MQ=255)
                                        tAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGa                                                                                                           <  1:1792700/148‑1 (MQ=255)
                                             aTCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCAcc                                                                                                      <  2:1728609/148‑1 (MQ=255)
                                                   cAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAAc                                                                                                >  2:831562/1‑148 (MQ=255)
                                                    aGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAccc                                                                                               >  2:2065414/1‑148 (MQ=255)
                                                    aGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAAcc                                                                                               >  2:1523815/1‑148 (MQ=255)
                                                      gCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGc                                                                                             >  2:3289754/1‑148 (MQ=255)
                                                       cACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGcc                                                                                            >  1:95/1‑148 (MQ=255)
                                                       cACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGcc                                                                                            >  2:1067289/1‑148 (MQ=255)
                                                        aCTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCa                                                                                           <  1:1032427/148‑1 (MQ=255)
                                                        aCTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCa                                                                                           <  2:347448/148‑1 (MQ=255)
                                                        aCTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCa                                                                                           <  1:1032402/148‑1 (MQ=255)
                                                             tCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTa                                                                                      >  1:2592875/1‑148 (MQ=255)
                                                              ccAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTaa                                                                                     <  2:1624246/148‑1 (MQ=255)
                                                                aGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACg                                                                                   >  2:1074697/1‑148 (MQ=255)
                                                                  ttCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACggg                                                                                 >  1:833158/1‑148 (MQ=255)
                                                                  ttCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACggg                                                                                 >  2:1534555/1‑148 (MQ=255)
                                                                  ttCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACggg                                                                                 >  1:861411/1‑148 (MQ=255)
                                                                   tCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGa                                                                                <  1:506097/148‑1 (MQ=255)
                                                                   tCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGa                                                                                >  2:1028710/1‑148 (MQ=255)
                                                                    cgcgcgCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGAt                                                                               >  2:317006/1‑148 (MQ=255)
                                                                     gcgcgcCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGAtc                                                                              <  1:2458303/148‑1 (MQ=255)
                                                                     gcgcgcCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGAtc                                                                              <  1:2457947/148‑1 (MQ=255)
                                                                      cgcgcCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGAtct                                                                             <  1:58214/148‑1 (MQ=255)
                                                                      cgcgcCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGAtct                                                                             <  2:2707369/148‑1 (MQ=255)
                                                                        cgcCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCg                                                                           >  2:2203410/1‑148 (MQ=255)
                                                                              aCAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTga                                                                     >  2:2817132/1‑148 (MQ=255)
                                                                               cAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTgag                                                                    <  2:1239030/148‑1 (MQ=255)
                                                                                  gTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTACCGGGATCTCGCCGTGAGata                                                                 >  2:1840127/1‑148 (MQ=255)
                                                                                  gTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGata                                                                 >  2:3276223/1‑148 (MQ=255)
                                                                                         cGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGcc                                                          >  2:1457696/1‑148 (MQ=255)
                                                                                           cTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCaa                                                        >  2:3212152/1‑148 (MQ=255)
                                                                                            tttCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAg                                                       >  2:1044826/1‑148 (MQ=255)
                                                                                             ttCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGc                                                      >  1:3130841/1‑148 (MQ=255)
                                                                                               cAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCaa                                                    <  2:1576748/148‑1 (MQ=255)
                                                                                               cAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCaa                                                    >  1:2155467/1‑148 (MQ=255)
                                                                                                  cAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAAcgc                                                 >  2:523475/1‑148 (MQ=255)
                                                                                                    aaTAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGt                                               <  2:147249/148‑1 (MQ=255)
                                                                                                        gCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAg                                           >  1:1885558/1‑148 (MQ=255)
                                                                                                         cTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGc                                          >  2:1798149/1‑148 (MQ=255)
                                                                                                          tGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCa                                         >  2:981152/1‑148 (MQ=255)
                                                                                                          tGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCa                                         >  2:2635306/1‑148 (MQ=255)
                                                                                                           gCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCAt                                        <  2:294670/148‑1 (MQ=255)
                                                                                                           gCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCAt                                        <  1:1074697/148‑1 (MQ=255)
                                                                                                           gCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCAt                                        <  1:2376612/148‑1 (MQ=255)
                                                                                                              aaCGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCAtctc                                     <  1:1961506/148‑1 (MQ=255)
                                                                                                                  aCCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCcaca                                 <  2:158005/148‑1 (MQ=255)
                                                                                                                  aCCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCcaca                                 <  2:158119/148‑1 (MQ=255)
                                                                                                                         aaCATTCGATTCATATCGCTATGGAAATTTGGGTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCg                          >  2:3155839/1‑148 (MQ=255)
                                                                                                                         aaCATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCg                          >  2:896203/1‑148 (MQ=255)
                                                                                                                          aCATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGt                         >  2:3078655/1‑148 (MQ=255)
                                                                                                                          aCATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGt                         >  2:3078437/1‑148 (MQ=255)
                                                                                                                                  ttCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTAc                 >  1:2674593/1‑148 (MQ=255)
                                                                                                                                  ttCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTAc                 >  1:671453/1‑148 (MQ=255)
                                                                                                                                     atatCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCgt              <  1:2666774/148‑1 (MQ=255)
                                                                                                                                     atatCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCgt              <  2:121607/148‑1 (MQ=255)
                                                                                                                                       atCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCgtgt            <  1:1044826/148‑1 (MQ=255)
                                                                                                                                         cGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTAt          <  1:1556198/148‑1 (MQ=255)
                                                                                                                                         cGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAAGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTAt          >  2:745143/1‑148 (MQ=255)
                                                                                                                                             aTGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTccc      >  2:2614031/1‑148 (MQ=255)
                                                                                                                                                 cAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTg  >  2:3224902/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTG  >  CP009273/1820743‑1821035

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
ATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTC  >  CP009273/1820740‑1821029
                                                                                                                                                    |                                                                                                                                             
ATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAA                                                                                                                                              >  SRR19973070.1674727/1‑150 (MQ=60)
     GATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAAC                                                                                                                                         <  SRR19973070.3025755/150‑1 (MQ=60)
      ATCAATGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACG                                                                                                                                        >  SRR19973070.1086199/1‑150 (MQ=60)
      ATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACG                                                                                                                                        <  SRR19973070.985019/150‑1 (MQ=60)
       TCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGT                                                                                                                                       <  SRR19973070.1749815/150‑1 (MQ=60)
        CAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGANAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTT                                                                                                                                      <  SRR19973070.1694784/150‑1 (MQ=60)
        CAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTT                                                                                                                                      >  SRR19973070.3173595/1‑150 (MQ=60)
        CAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACGGTTCCTGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTT                                                                                                                                      <  SRR19973070.1160306/150‑1 (MQ=60)
         ACGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTT                                                                                                                                     <  SRR19973070.1327792/150‑1 (MQ=60)
         AAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTT                                                                                                                                     >  SRR19973070.1766267/1‑150 (MQ=60)
         AAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTT                                                                                                                                     >  SRR19973070.1766266/1‑150 (MQ=60)
         AAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTT                                                                                                                                     <  SRR19973070.958667/150‑1 (MQ=60)
         AAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTT                                                                                                                                     <  SRR19973070.2325110/150‑1 (MQ=60)
         AAGGTGCCAGCGCACAGATTCTTTGCCCTGGAAAAAAAAATCTTCCAGCCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGATTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTT                                                                                                                                     <  SRR19973070.956713/150‑1 (MQ=60)
          AGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTC                                                                                                                                    >  SRR19973070.3321666/1‑150 (MQ=60)
           ggcgcccgcccccagccccttgGCCCTCGTCATAAAAATCTACCCGGCACTGTTCCAGTTCGCGCCCCCGCCAGGTTTCTCCGCTTTCACCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCC                                                                                                                                   <  SRR19973070.2887589/128‑1 (MQ=60)
              GCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTT                                                                                                                                <  SRR19973070.1759428/150‑1 (MQ=60)
               CCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTG                                                                                                                               >  SRR19973070.1508475/1‑150 (MQ=60)
                 AGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCT                                                                                                                             <  SRR19973070.2479823/150‑1 (MQ=60)
                  GCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTT                                                                                                                            <  SRR19973070.2631492/150‑1 (MQ=60)
                     CACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAG                                                                                                                         >  SRR19973070.1702474/1‑150 (MQ=60)
                        AGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGC                                                                                                                      >  SRR19973070.1667642/1‑150 (MQ=60)
                        AGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGC                                                                                                                      >  SRR19973070.949655/1‑150 (MQ=60)
                         GATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCA                                                                                                                     >  SRR19973070.2743034/1‑150 (MQ=60)
                             CTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAG                                                                                                                 >  SRR19973070.122690/1‑150 (MQ=60)
                                   CCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCC                                                                                                           <  SRR19973070.637233/150‑1 (MQ=60)
                                   CCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCC                                                                                                           >  SRR19973070.224667/1‑150 (MQ=60)
                                    CTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACTGTTTTCTCCGATTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCC                                                                                                          <  SRR19973070.350249/150‑1 (MQ=60)
                                       GTAATAAAAATCTTCCAGGCACTGTTCCTGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAG                                                                                                       <  SRR19973070.309828/150‑1 (MQ=60)
                                          ATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGAT                                                                                                    <  SRR19973070.1593414/150‑1 (MQ=60)
                                           TAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATC                                                                                                   <  SRR19973070.1805087/150‑1 (MQ=60)
                                           TAAAAATCTTCCAGGCACTGGTCCTGTTCGCGCGCCCGACTGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATC                                                                                                   <  SRR19973070.1805112/150‑1 (MQ=60)
                                             AAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCAC                                                                                                 <  SRR19973070.59179/150‑1 (MQ=60)
                                              AAATCTTCCCGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCGGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACC                                                                                                <  SRR19973070.370203/150‑1 (MQ=60)
                                                    TCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAAC                                                                                          <  SRR19973070.2736677/150‑1 (MQ=60)
                                                    TCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAAC                                                                                          <  SRR19973070.2414604/150‑1 (MQ=60)
                                                        GGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCC                                                                                      >  SRR19973070.96/1‑150 (MQ=60)
                                                           ACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAAC                                                                                   <  SRR19973070.1039731/150‑1 (MQ=60)
                                                           ACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAAC                                                                                   <  SRR19973070.1039706/150‑1 (MQ=60)
                                                              GTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTA                                                                                >  SRR19973070.2609819/1‑150 (MQ=60)
                                                              GTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTA                                                                                >  SRR19973070.1590856/1‑150 (MQ=60)
                                                                TCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAAC                                                                              <  SRR19973070.922743/150‑1 (MQ=60)
                                                                 CCATTTCGCGCGCCCGACAGGTTTCTCCGTTTTCAGAAAATAGCTGCAAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACG                                                                             <  SRR19973070.2801177/150‑1 (MQ=60)
                                                                   AGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGG                                                                           >  SRR19973070.867395/1‑150 (MQ=60)
                                                                   AGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGG                                                                           >  SRR19973070.838800/1‑150 (MQ=60)
                                                                   AGTTCGCGCGCCCGACAGGTTTCTCCGCTGTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGG                                                                           >  SRR19973070.867389/1‑150 (MQ=60)
                                                                    GTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGA                                                                          <  SRR19973070.2497130/150‑1 (MQ=60)
                                                                      TCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATC                                                                        <  SRR19973070.509595/150‑1 (MQ=60)
                                                                       CGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCT                                                                       <  SRR19973070.1604030/150‑1 (MQ=60)
                                                                        GCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTC                                                                      <  SRR19973070.2474340/150‑1 (MQ=60)
                                                                        GCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTC                                                                      <  SRR19973070.2474366/150‑1 (MQ=60)
                                                                        GCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTC                                                                      <  SRR19973070.2474007/150‑1 (MQ=60)
                                                                         CGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCG                                                                     <  SRR19973070.58889/150‑1 (MQ=60)
                                                                                     GTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATC                                                         >  SRR19973070.1797292/1‑150 (MQ=60)
                                                                                              CTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGC                                                >  SRR19973070.3151067/1‑150 (MQ=60)
                                                                                                TTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAA                                              >  SRR19973070.2169706/1‑150 (MQ=60)
                                                                                                TTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAA                                              >  SRR19973070.2169803/1‑150 (MQ=60)
                                                                                                TTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAA                                              >  SRR19973070.2170089/1‑150 (MQ=60)
                                                                                                         TAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAG                                     >  SRR19973070.1898362/1‑150 (MQ=60)
                                                                                                              GCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCT                                <  SRR19973070.1082186/150‑1 (MQ=60)
                                                                                                              GCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCT                                <  SRR19973070.2392317/150‑1 (MQ=60)
                                                                                                                CCACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCC                              <  SRR19973070.2256267/150‑1 (MQ=60)
                                                                                                                 AACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCA                             <  SRR19973070.1974727/150‑1 (MQ=60)
                                                                                                                               ATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCA               <  SRR19973070.1122433/150‑1 (MQ=60)
                                                                                                                                   GATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTAC           >  SRR19973070.2691958/1‑150 (MQ=60)
                                                                                                                                   GATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTAC           >  SRR19973070.675836/1‑150 (MQ=60)
                                                                                                                                   GATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTAC           >  SRR19973070.676247/1‑150 (MQ=60)
                                                                                                                                        ATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGT      <  SRR19973070.2684095/150‑1 (MQ=60)
                                                                                                                                          ATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTAT    <  SRR19973070.1052189/150‑1 (MQ=60)
                                                                                                                                          ATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTAT    <  SRR19973070.1051623/150‑1 (MQ=60)
                                                                                                                                           TCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATT   <  SRR19973070.1077160/150‑1 (MQ=60)
                                                                                                                                            CGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTC  <  SRR19973070.1567131/150‑1 (MQ=60)
                                                                                                                                            CGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTC  <  SRR19973070.1567142/150‑1 (MQ=60)
                                                                                                                                                    |                                                                                                                                             
ATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTC  >  CP009273/1820740‑1821029

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