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 F2 I199 R1 220 57.5 3163680 97.0% 3068769 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,898,828 (T)5→4 coding (130/1878 nt) bglF ← PTS beta‑glucoside transporter subunit IIABC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,898,8240T.98.8% 346.4 / ‑1.9 83coding (134/1878 nt)bglFPTS beta‑glucoside transporter subunit IIABC
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (50/32);  total (51/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.85e-01

TAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTT  >  NZ_CP009273/3898735‑3898909
                                                                                         |                                                                                     
tAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCt                                                                                       >  2:1461937/1‑90 (MQ=255)
tAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCt                                                                                       <  2:48147/90‑1 (MQ=255)
tAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCt                                                                                       >  1:936376/1‑90 (MQ=255)
     gCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCttttca                                                                                  >  1:1544322/1‑88 (MQ=255)
      ccAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCttttcag                                                                                 <  2:952099/90‑4 (MQ=255)
      ccAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCttttcag                                                                                 >  2:614713/1‑87 (MQ=255)
      ccAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCttttcag                                                                                 >  1:92839/1‑87 (MQ=255)
       cAGGAAGACATCGGCCACATGGTTACCTATGACCCCCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCttttcagt                                                                                <  1:324436/90‑5 (MQ=255)
       cAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATCCCGGGGGTCttttcagt                                                                                <  2:92839/90‑5 (MQ=255)
       cAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCttttcagt                                                                                <  2:607139/90‑5 (MQ=255)
        aGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTa                                                                              >  1:1201373/1‑90 (MQ=255)
        aGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTa                                                                              >  2:1053598/1‑90 (MQ=255)
        aGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTa                                                                              >  2:461647/1‑90 (MQ=255)
          gAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTAcc                                                                            <  1:59736/90‑1 (MQ=255)
                       acaTGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCgctt                                                               <  2:1460776/90‑1 (MQ=255)
                       acaTGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCgctt                                                               <  1:1404288/90‑1 (MQ=255)
                              tACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTc                                                        >  2:1542595/1‑90 (MQ=255)
                              tACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTc                                                        >  1:1371584/1‑90 (MQ=255)
                                ccTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCAt                                                      <  2:964167/90‑1 (MQ=255)
                                  tATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATcc                                                    >  2:1109538/1‑90 (MQ=255)
                                  tATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATcc                                                    >  2:113409/1‑90 (MQ=255)
                                  tATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATcc                                                    >  1:155912/1‑90 (MQ=255)
                                    tGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCtt                                                  >  1:509208/1‑90 (MQ=255)
                                     gACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCttt                                                 >  2:1150534/1‑90 (MQ=255)
                                     gACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGACTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCttt                                                 >  1:537643/1‑90 (MQ=255)
                                         accTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCtttaatt                                             >  1:617961/1‑90 (MQ=255)
                                         accTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCtttaatt                                             >  1:906594/1‑90 (MQ=255)
                                         accTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCtttaatt                                             >  2:330222/1‑90 (MQ=255)
                                          ccTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACct                                                                           <  2:331776/59‑1 (MQ=255)
                                          ccTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACct                                                                           >  1:331776/1‑59 (MQ=255)
                                          ccTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCtttaattt                                            >  1:1299281/1‑90 (MQ=255)
                                          ccTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCtttaattt                                            >  1:1164463/1‑90 (MQ=255)
                                          ccTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCtttaattt                                            >  2:344248/1‑90 (MQ=255)
                                            tGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCtttaatttaa                                          <  1:330222/90‑1 (MQ=255)
                                               aaCTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAAc                                       <  1:113409/90‑1 (MQ=255)
                                                   ggCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTaa                                   >  2:1180246/1‑90 (MQ=255)
                                                   ggCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTaa                                   >  1:1364058/1‑90 (MQ=255)
                                                    gCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAAt                                  >  2:623789/1‑90 (MQ=255)
                                                    gCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAAt                                  >  2:1377596/1‑90 (MQ=255)
                                                    gCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAAt                                  >  1:864087/1‑90 (MQ=255)
                                                    gCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAAt                                  <  1:1121677/90‑1 (MQ=255)
                                                    gCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAAt                                  >  1:926986/1‑90 (MQ=255)
                                                    gCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAAt                                  >  1:1179739/1‑90 (MQ=255)
                                                    gCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAAt                                  >  1:476975/1‑90 (MQ=255)
                                                    gCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAAt                                  <  1:189074/90‑1 (MQ=255)
                                                    gCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAAt                                  <  2:136853/90‑1 (MQ=255)
                                                    gCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAAt                                  >  2:647840/1‑90 (MQ=255)
                                                    gCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAAt                                  >  2:1286787/1‑90 (MQ=255)
                                                        ccGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATcgcg                              >  2:1139699/1‑90 (MQ=255)
                                                         cGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGt                             <  2:314365/90‑1 (MQ=255)
                                                         cGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGt                             <  1:461647/90‑1 (MQ=255)
                                                         cGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGt                             <  1:373360/90‑1 (MQ=255)
                                                         cGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGt                             <  1:1266762/90‑1 (MQ=255)
                                                              tCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTgcgc                        <  1:1508817/90‑1 (MQ=255)
                                                              tCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTgcgc                        <  2:1310863/90‑1 (MQ=255)
                                                              tCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTgcgc                        <  1:1118572/90‑1 (MQ=255)
                                                                caccaccaTAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCaa                      <  2:497650/90‑1 (MQ=255)
                                                                caccaccaTAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCaa                      <  2:524936/90‑1 (MQ=255)
                                                                caccaccaTAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCaa                      <  1:1542595/90‑1 (MQ=255)
                                                                    accaTAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGc                   >  1:1349178/1‑90 (MQ=255)
                                                                     ccaTAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCaa                      >  1:1280597/1‑85 (MQ=255)
                                                                     ccaTAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCaa                      <  2:1280597/85‑1 (MQ=255)
                                                                     ccaTAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCAt                 >  1:157673/1‑90 (MQ=255)
                                                                     ccaTAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCAt                 >  2:1260731/1‑90 (MQ=255)
                                                                       ataataataCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCa               >  1:316259/1‑90 (MQ=255)
                                                                       ataataataCCGGGGGTCTTTTCAGTACCCTGCTTGCGCTTTGCTTTCATCCttt                                                 >  1:112521/1‑55 (MQ=255)
                                                                       ataataataCCGGGGGTCTTTTCAGTACCCTGCTTGCGCTTTGCTTTCATCCttt                                                 <  2:112521/55‑1 (MQ=255)
                                                                       ataataataCCGGGGGACTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTTCGCAATGCATCa               >  1:1540099/1‑90 (MQ=255)
                                                                        taataataCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAg              <  2:1234621/90‑1 (MQ=255)
                                                                        taataataCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAg              <  2:672451/90‑1 (MQ=255)
                                                                             ataCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTca         <  1:1150534/90‑1 (MQ=255)
                                                                             ataCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTca         <  1:890628/90‑1 (MQ=255)
                                                                             ataCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTca         <  1:1009024/90‑1 (MQ=255)
                                                                             ataCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTca         <  1:205318/90‑1 (MQ=255)
                                                                              taCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTcac        <  1:1354748/90‑1 (MQ=255)
                                                                               aCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTcaca       >  2:747792/1‑90 (MQ=255)
                                                                               aCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTcaca       >  1:698505/1‑90 (MQ=255)
                                                                               aCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTcaca       <  1:639533/90‑1 (MQ=255)
                                                                               aCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTcaca       <  1:1180246/90‑1 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGGAATGCATCAGACTCACAATg    >  1:894500/1‑90 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAg              >  1:895426/1‑80 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAg              <  2:895426/80‑1 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg    >  2:385862/1‑90 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg    >  2:1087629/1‑90 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg    >  1:671709/1‑90 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg    >  1:632134/1‑90 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg    >  1:5283/1‑90 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg    >  1:127050/1‑90 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg    >  1:1226256/1‑90 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg    >  1:1067050/1‑90 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg    >  1:1023884/1‑90 (MQ=255)
                                                                                   ggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGt   >  1:1497804/1‑90 (MQ=255)
                                                                                    gggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGtt  >  2:375327/1‑90 (MQ=255)
                                                                                         |                                                                                     
TAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTT  >  NZ_CP009273/3898735‑3898909

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCC  >  NZ_CP009273/3898729‑3898927
                                                                                                   |                                                                                                   
CACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTC                                                                                                     <  SRR3722078.1477507/100‑1 (MQ=60)
  CTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAG                                                                                                  >  SRR3722078.93928/1‑100 (MQ=60)
    GTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTA                                                                                                >  SRR3722078.325207/1‑100 (MQ=60)
    GTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTA                                                                                                >  SRR3722078.1220217/1‑100 (MQ=60)
     TTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTAC                                                                                               <  SRR3722078.1339481/100‑1 (MQ=60)
            ccatcgtcggcagcgtcagatgtgtataagagacagCCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCT                                                                                        <  SRR3722078.1381736/64‑1 (MQ=60)
             CAGGAAGACATCGGCCACATGGTTACCTATGACCCCCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTT                                                                                       <  SRR3722078.328414/100‑1 (MQ=60)
                GAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCG                                                                                    <  SRR3722078.60441/100‑1 (MQ=60)
                          GCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTC                                                                          >  SRR3722078.1393673/1‑100 (MQ=60)
                             ACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATC                                                                       <  SRR3722078.1427046/100‑1 (MQ=60)
                              CATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCC                                                                      >  SRR3722078.157731/1‑100 (MQ=60)
                                TGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTT                                                                    >  SRR3722078.515234/1‑100 (MQ=60)
                                 GGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGACTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTT                                                                   >  SRR3722078.543997/1‑100 (MQ=60)
                                     ACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATT                                                               >  SRR3722078.625085/1‑100 (MQ=60)
                                     ACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATT                                                               >  SRR3722078.918791/1‑100 (MQ=60)
                                      CCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCctgtctcttatacacatctga                                                              >  SRR3722078.335833/1‑79 (MQ=60)
                                      CCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTT                                                              >  SRR3722078.1319971/1‑100 (MQ=60)
                                      CCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTT                                                              >  SRR3722078.1182537/1‑100 (MQ=60)
                                               ACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAA                                                     >  SRR3722078.1386000/1‑100 (MQ=60)
                                                CCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAAT                                                    >  SRR3722078.875656/1‑100 (MQ=60)
                                                CCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAAT                                                    >  SRR3722078.939723/1‑100 (MQ=60)
                                                CCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAAT                                                    >  SRR3722078.1198155/1‑100 (MQ=60)
                                                CCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAAT                                                    >  SRR3722078.482634/1‑100 (MQ=60)
                                                  TGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCG                                                  <  SRR3722078.334263/100‑1 (MQ=60)
                                                     AACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGT                                               <  SRR3722078.114745/100‑1 (MQ=60)
                                                          GCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGC                                          <  SRR3722078.191318/100‑1 (MQ=60)
                                                          GCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGC                                          <  SRR3722078.1138800/100‑1 (MQ=60)
                                                               CGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGC                                     <  SRR3722078.467163/100‑1 (MQ=60)
                                                               CGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGC                                     <  SRR3722078.1286863/100‑1 (MQ=60)
                                                               CGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGC                                     <  SRR3722078.377925/100‑1 (MQ=60)
                                                                GCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGC                                     >  SRR3722078.1370832/1‑100 (MQ=60)
                                                                 CTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCAT                                   >  SRR3722078.159512/1‑100 (MQ=60)
                                                                 CTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCAT                                   >  SRR3722078.1300966/1‑100 (MQ=60)
                                                                   TTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCA                                 >  SRR3722078.320131/1‑100 (MQ=60)
                                                                   TTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACctgtctcttatacacatctgacgct                                >  SRR3722078.113851/1‑75 (MQ=60)
                                                                   TTCCACCACCATAATAATACCGGGGGACTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTTCGCAATGCATCA                                 >  SRR3722078.1565207/1‑100 (MQ=60)
                                                                    TCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAG                                <  SRR3722078.1533422/100‑1 (MQ=60)
                                                                    TCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAG                                <  SRR3722078.1135652/100‑1 (MQ=60)
                                                                      CACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGAC                              <  SRR3722078.1567744/100‑1 (MQ=60)
                                                                           CCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACA                         >  SRR3722078.706895/1‑100 (MQ=60)
                                                                              TAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGGAATGCATCAGACTCACAATG                      >  SRR3722078.906523/1‑100 (MQ=60)
                                                                              TAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATG                      >  SRR3722078.639404/1‑100 (MQ=60)
                                                                              TAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATG                      >  SRR3722078.679376/1‑100 (MQ=60)
                                                                              TAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATG                      >  SRR3722078.907458/1‑100 (MQ=60)
                                                                              TAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATG                      >  SRR3722078.5339/1‑100 (MQ=60)
                                                                              TAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATG                      >  SRR3722078.1039054/1‑100 (MQ=60)
                                                                              TAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATG                      >  SRR3722078.1083144/1‑100 (MQ=60)
                                                                              TAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATG                      >  SRR3722078.128517/1‑100 (MQ=60)
                                                                              TAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATG                      >  SRR3722078.1245570/1‑100 (MQ=60)
                                                                               TATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGT                     >  SRR3722078.1522222/1‑100 (MQ=60)
                                                                                   ATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATC                 <  SRR3722078.1023845/100‑1 (MQ=60)
                                                                                   ATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATC                 <  SRR3722078.207795/100‑1 (MQ=60)
                                                                                   ATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATC                 <  SRR3722078.902606/100‑1 (MQ=60)
                                                                                   ATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATC                 <  SRR3722078.1168276/100‑1 (MQ=60)
                                                                                    gtataagagacagTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATC                 <  SRR3722078.988777/87‑1 (MQ=60)
                                                                                    TACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCT                >  SRR3722078.1350768/1‑100 (MQ=60)
                                                                                    TACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCC                <  SRR3722078.1376494/100‑1 (MQ=60)
                                                                                     ACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTgtctcttatacacatctgacgctgccgacgatggcgggcgtgtaggtttcggtg               >  SRR3722078.515577/1‑46 (MQ=60)
                                                                                     ACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTG               <  SRR3722078.646889/100‑1 (MQ=60)
                                                                                     ACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTG               <  SRR3722078.1198676/100‑1 (MQ=60)
                                                                                           GGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCC         >  SRR3722078.514908/1‑100 (MQ=60)
                                                                                                 TTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACT    >  SRR3722078.380714/1‑100 (MQ=60)
                                                                                                   TCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCC  <  SRR3722078.757308/100‑1 (MQ=60)
                                                                                                   TCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCC  <  SRR3722078.381087/100‑1 (MQ=60)
                                                                                                   |                                                                                                   
CACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCC  >  NZ_CP009273/3898729‑3898927

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