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 F1 I2 R1 1873 183.1 9376148 98.2% 9207377 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,998,039 A→G Y124H (TAC→CAC)  fliY ← putative periplasmic binding transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,998,0390AG99.2% 380.2 / NA 122Y124H (TAC→CAC) fliYputative periplasmic binding transport protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (64/57);  total (64/58)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

CCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGA  >  NC_000913/1997949‑1998126
                                                                                          |                                                                                       
cccccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCg                                                                                       >  1:3991489/1‑93 (MQ=255)
 ccccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCgg                                                                                      >  1:296106/1‑93 (MQ=255)
 ccccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCgg                                                                                      >  2:3356826/1‑93 (MQ=255)
  cccTCTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCggg                                                                                     <  1:1727990/93‑1 (MQ=255)
  cccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCggg                                                                                     >  2:2829047/1‑93 (MQ=255)
  cccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCggg                                                                                     >  1:870728/1‑93 (MQ=255)
   ccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCgggg                                                                                    <  1:4367975/93‑1 (MQ=255)
   ccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCgggg                                                                                    >  2:2487565/1‑93 (MQ=255)
     aCTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGtt                                                                                  >  1:1527881/1‑93 (MQ=255)
     aCTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGtt                                                                                  <  2:1098801/93‑1 (MQ=255)
      cTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTTGTTACCTTTTTTCACCAGCGGCCGAATACCTGAAATGGTGCGCGGGGCTg                                                                                 >  2:1225133/1‑93 (MQ=255)
       ttttttGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTga                                                                                >  1:2501142/1‑93 (MQ=255)
         ttttGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTgaga                                                                              >  1:1571566/1‑93 (MQ=255)
           ttGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGaaa                                                                            <  2:970311/93‑1 (MQ=255)
           ttGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGaaa                                                                            >  2:965884/1‑93 (MQ=255)
             gCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATc                                                                          <  2:1883156/93‑1 (MQ=255)
             gCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATc                                                                          <  1:1207083/93‑1 (MQ=255)
                tttCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTa                                                                       >  2:3112864/1‑93 (MQ=255)
                tttCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGGTGAGAAATCGTa                                                                       >  1:2261651/1‑93 (MQ=255)
                 ttCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTAt                                                                      <  2:2658851/93‑1 (MQ=255)
                    aGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTAtttt                                                                   <  2:1850789/93‑1 (MQ=255)
                    aGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTAtttt                                                                   <  1:1437213/93‑1 (MQ=255)
                       tcatcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTc                                                                >  1:2562213/1‑93 (MQ=255)
                        catcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCt                                                               >  2:4521152/1‑93 (MQ=255)
                        catcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCt                                                               <  1:3849358/93‑1 (MQ=255)
                        catcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCt                                                               >  2:4269485/1‑93 (MQ=255)
                          tcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTg                                                             >  2:3151699/1‑93 (MQ=255)
                          tcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTg                                                             >  2:2740650/1‑93 (MQ=255)
                          tcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTg                                                             >  2:2739028/1‑93 (MQ=255)
                          tcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTg                                                             >  2:738029/1‑93 (MQ=255)
                          tcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTg                                                             >  1:3026705/1‑93 (MQ=255)
                           cGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTgc                                                            <  1:3165768/93‑1 (MQ=255)
                             gCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTgcgc                                                          >  1:1698456/1‑93 (MQ=255)
                             gCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTgcgc                                                          >  1:4204097/1‑93 (MQ=255)
                               tGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtc                                                        >  1:565557/1‑93 (MQ=255)
                               tGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtc                                                        >  1:2420038/1‑93 (MQ=255)
                                gTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtca                                                       <  2:4040562/93‑1 (MQ=255)
                                gTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtca                                                       >  2:1125860/1‑93 (MQ=255)
                                gTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtca                                                       <  2:2161830/93‑1 (MQ=255)
                                gTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtca                                                       <  2:4040561/93‑1 (MQ=255)
                                 ttttAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtcat                                                      <  2:584254/93‑1 (MQ=255)
                                  tttAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtcatc                                                     <  1:1577080/93‑1 (MQ=255)
                                   ttAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAATTCGTATTTTTTCTTGCGCtcatca                                                    <  2:4518312/93‑1 (MQ=255)
                                   ttAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtcatca                                                    <  2:3296868/93‑1 (MQ=255)
                                    tAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAg                                                   <  1:3049078/93‑1 (MQ=255)
                                      aTGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGaa                                                 <  1:822193/93‑1 (MQ=255)
                                      aTGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGaa                                                 >  1:383797/1‑93 (MQ=255)
                                       tGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGaaa                                                >  1:210065/1‑93 (MQ=255)
                                       tGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGaaa                                                >  1:2761554/1‑93 (MQ=255)
                                       tGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGaaa                                                >  2:2580084/1‑93 (MQ=255)
                                       tGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGCGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGaaa                                                <  1:2507254/93‑1 (MQ=255)
                                       tGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATAGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGaaa                                                <  2:3795981/93‑1 (MQ=255)
                                        ggTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAAt                                               >  1:2083974/1‑93 (MQ=255)
                                         gTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATg                                              >  2:2364074/1‑93 (MQ=255)
                                           gCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGt                                            <  1:2056424/93‑1 (MQ=255)
                                            ccTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTg                                           >  2:459237/1‑93 (MQ=255)
                                             cTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGa                                          <  2:3042617/93‑1 (MQ=255)
                                             cTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGa                                          >  2:3238739/1‑93 (MQ=255)
                                              ttCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGAc                                         >  1:3392628/1‑93 (MQ=255)
                                               tCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGAcc                                        <  2:769425/93‑1 (MQ=255)
                                               tCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGAcc                                        <  1:1469068/93‑1 (MQ=255)
                                               tCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGAcc                                        >  1:562417/1‑93 (MQ=255)
                                                cGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCt                                       >  1:521775/1‑93 (MQ=255)
                                                cGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCt                                       <  1:216285/93‑1 (MQ=255)
                                                 gTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTg                                      >  2:1132605/1‑93 (MQ=255)
                                                 gTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTg                                      >  1:1747263/1‑93 (MQ=255)
                                                  ttcccTTTTTTCGCCAGCGCCTGAATACCAGAAATGGTGTGCGGTGTTGAGAAATCGTATTTTTTCTGGCGCTCATCAGAAATGGTGACCTGa                                     <  2:2044379/90‑1 (MQ=255)
                                                  ttACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGGAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGa                                     >  1:746308/1‑93 (MQ=255)
                                                  ttACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGa                                     <  1:3381215/93‑1 (MQ=255)
                                                  ttACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGa                                     <  1:2406350/93‑1 (MQ=255)
                                                   tACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGAt                                    <  1:3040802/93‑1 (MQ=255)
                                                   tACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGAt                                    <  2:4674138/93‑1 (MQ=255)
                                                   tACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCCGAt                                    <  2:4063884/93‑1 (MQ=255)
                                                     ccTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTa                                  >  2:3602549/1‑93 (MQ=255)
                                                      cTTTTTTCACCAGCGCGTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTaa                                 <  2:1569766/93‑1 (MQ=255)
                                                       ttttttCCCCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAAt                                <  1:729721/93‑1 (MQ=255)
                                                       ttttttCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAAt                                <  2:2670966/93‑1 (MQ=255)
                                                       ttttttCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAAt                                >  1:2232700/1‑93 (MQ=255)
                                                       ttttttCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGGTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAAt                                >  2:1933929/1‑93 (MQ=255)
                                                        tttttCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATc                               <  1:543764/93‑1 (MQ=255)
                                                        tttttCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATc                               <  1:4180337/93‑1 (MQ=255)
                                                          tttCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATcac                             >  2:867667/1‑93 (MQ=255)
                                                          tttCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATcac                             <  1:279379/93‑1 (MQ=255)
                                                              aCCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACcaca                         >  1:1430559/1‑93 (MQ=255)
                                                              aCCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACcaca                         >  1:181417/1‑93 (MQ=255)
                                                                cAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATc                       <  2:4603167/93‑1 (MQ=255)
                                                                  gcgcCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCaa                     <  2:3739501/93‑1 (MQ=255)
                                                                       tGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGt                >  2:947450/1‑93 (MQ=255)
                                                                       tGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGt                >  1:4472508/1‑93 (MQ=255)
                                                                          aTACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTa             >  2:3320805/1‑93 (MQ=255)
                                                                           tACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTag            >  1:1599306/1‑93 (MQ=255)
                                                                             ccAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTagag          <  1:2865799/93‑1 (MQ=255)
                                                                             ccAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTagag          <  2:3023767/93‑1 (MQ=255)
                                                                               aGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTc        <  2:2821170/93‑1 (MQ=255)
                                                                               aGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTc        <  2:3832486/93‑1 (MQ=255)
                                                                                gAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTcc       >  2:2587933/1‑93 (MQ=255)
                                                                                gAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTcc       <  1:3846927/93‑1 (MQ=255)
                                                                                 aaaTGTTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCa      <  1:2940804/93‑1 (MQ=255)
                                                                                     ggTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCGGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGa  >  1:2512474/1‑93 (MQ=255)
                                                                                     ggTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGa  <  1:357791/93‑1 (MQ=255)
                                                                                          |                                                                                       
CCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGA  >  NC_000913/1997949‑1998126

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCA  >  NC_000913/1997948‑1998131
                                                                                           |                                                                                            
ACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCG                                                                                            >  SRR4302535.4018790/1‑94 (MQ=60)
 CCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGG                                                                                           >  SRR4302535.298170/1‑94 (MQ=60)
  CCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGG                                                                                          >  SRR4302535.876894/1‑94 (MQ=60)
   CCCTCTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGG                                                                                         <  SRR4302535.1740480/94‑1 (MQ=60)
    CCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGT                                                                                        <  SRR4302535.4397865/94‑1 (MQ=60)
     CACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTT                                                                                       >  SRR4302535.1538887/1‑94 (MQ=60)
       CTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGA                                                                                     >  SRR4302535.2519181/1‑94 (MQ=60)
         TTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGA                                                                                   >  SRR4302535.1582910/1‑94 (MQ=60)
             TGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATC                                                                               >  SRR4302535.3749151/1‑94 (MQ=60)
              GCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCG                                                                              <  SRR4302535.1215758/94‑1 (MQ=60)
                CTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGGTGAGAAATCGTA                                                                            >  SRR4302535.2277763/1‑94 (MQ=60)
                   TCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTT                                                                         <  SRR4302535.4160149/94‑1 (MQ=60)
                     AGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTT                                                                       <  SRR4302535.1447485/94‑1 (MQ=60)
                     AGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATAGCATTTTT                                                                       >  SRR4302535.2751827/1‑94 (MQ=60)
                       ATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTC                                                                     >  SRR4302535.2580656/1‑94 (MQ=60)
                         CATCGGCTGTTTTAATGGTGCCTTGCTAACCTTTTTTACCAGCGCCTGAATTACACGAATGGTGTGGGGGGTTGAGAAATTGGATTTTTTCTTT                                                                  >  SRR4302535.2069837/1‑94 (MQ=18)
                         CATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTT                                                                   <  SRR4302535.3875849/94‑1 (MQ=60)
                          ATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTG                                                                  >  SRR4302535.3047894/1‑94 (MQ=60)
                            CGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCG                                                                <  SRR4302535.3187710/94‑1 (MQ=60)
                             GGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGC                                                               >  SRR4302535.1710732/1‑94 (MQ=60)
                             GGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGC                                                               >  SRR4302535.4232871/1‑94 (MQ=60)
                               CTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTC                                                             >  SRR4302535.569519/1‑94 (MQ=60)
                               CTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTC                                                             >  SRR4302535.2437463/1‑94 (MQ=60)
                                   TTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCA                                                         <  SRR4302535.1588461/94‑1 (MQ=60)
                                     TAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGA                                                       >  SRR4302535.1865563/1‑94 (MQ=60)
                                     TAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGA                                                       <  SRR4302535.3070436/94‑1 (MQ=60)
                                      AATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAA                                                      >  SRR4302535.386456/1‑94 (MQ=60)
                                       ATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAA                                                     <  SRR4302535.827956/94‑1 (MQ=60)
                                       ATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAA                                                     >  SRR4302535.211509/1‑94 (MQ=60)
                                       ATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAA                                                     >  SRR4302535.2781233/1‑94 (MQ=60)
                                        TGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAAT                                                    >  SRR4302535.2098980/1‑94 (MQ=60)
                                        TGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGCGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAAT                                                    <  SRR4302535.2525329/94‑1 (MQ=60)
                                            GCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTG                                                <  SRR4302535.2071238/94‑1 (MQ=60)
                                              CTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGAC                                              >  SRR4302535.3416232/1‑94 (MQ=60)
                                               TTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACC                                             >  SRR4302535.566350/1‑94 (MQ=60)
                                                TCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCT                                            <  SRR4302535.1479592/94‑1 (MQ=60)
                                                TCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCT                                            >  SRR4302535.525425/1‑94 (MQ=60)
                                                 CGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTG                                           <  SRR4302535.2204784/94‑1 (MQ=60)
                                                 CGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTG                                           <  SRR4302535.217769/94‑1 (MQ=60)
                                                 CGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTG                                           >  SRR4302535.1759900/1‑94 (MQ=60)
                                                  GTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGGAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGA                                          >  SRR4302535.751557/1‑94 (MQ=60)
                                                   TTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGAT                                         >  SRR4302535.4590463/1‑94 (MQ=60)
                                                   TTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGAT                                         <  SRR4302535.2423673/94‑1 (MQ=60)
                                                   TTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGAT                                         <  SRR4302535.3404721/94‑1 (MQ=60)
                                                    TACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATT                                        >  SRR4302535.1628556/1‑94 (MQ=60)
                                                    TACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATT                                        <  SRR4302535.3062103/94‑1 (MQ=60)
                                                       CTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAAT                                     <  SRR4302535.2976925/94‑1 (MQ=60)
                                                       CTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAAT                                     >  SRR4302535.2248575/1‑94 (MQ=60)
                                                        TTTTTTCCCCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATC                                    <  SRR4302535.734841/94‑1 (MQ=60)
                                                         TTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCA                                   <  SRR4302535.547566/94‑1 (MQ=60)
                                                         TTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCA                                   <  SRR4302535.4208947/94‑1 (MQ=60)
                                                          TTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCAC                                  >  SRR4302535.57050/1‑94 (MQ=60)
                                                           TTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACC                                 <  SRR4302535.281328/94‑1 (MQ=60)
                                                             TCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCAC                               <  SRR4302535.2510601/94‑1 (MQ=60)
                                                              CACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACA                              >  SRR4302535.1440766/1‑94 (MQ=60)
                                                              CACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACA                              >  SRR4302535.182675/1‑94 (MQ=60)
                                                                       CTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGT                     >  SRR4302535.4503203/1‑94 (MQ=60)
                                                                       CTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGT                     >  SRR4302535.3160230/1‑94 (MQ=60)
                                                                           ATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAG                 >  SRR4302535.1610839/1‑94 (MQ=60)
                                                                             ACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAG               >  SRR4302535.152402/1‑94 (MQ=60)
                                                                              CCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGT              <  SRR4302535.2886099/94‑1 (MQ=60)
                                                                                 GAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCA           <  SRR4302535.3873411/94‑1 (MQ=60)
                                                                                  AAATGTTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAG          <  SRR4302535.2961480/94‑1 (MQ=60)
                                                                                     TGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCGGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGA       >  SRR4302535.2530587/1‑94 (MQ=60)
                                                                                      GGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGAC      <  SRR4302535.360265/94‑1 (MQ=60)
                                                                                          TGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCA  >  SRR4302535.3429928/1‑94 (MQ=60)
                                                                                           |                                                                                            
ACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCA  >  NC_000913/1997948‑1998131

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