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 I1 R1 1357 179.3 9500450 96.1% 9129932 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,787,410 A→G T318A (ACT→GCT)  aroH → 3‑deoxy‑D‑arabinoheptulosonate‑7‑phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,787,4100AG99.0% 334.9 / NA 104T318A (ACT→GCT) aroH3‑deoxy‑D‑arabinoheptulosonate‑7‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (22/81);  minor base T (0/1);  total (22/82)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

ATCCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGC  >  NC_000913/1787320‑1787499
                                                                                          |                                                                                         
aTCCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCcgctcgct                                                                                         <  1:3389383/93‑1 (MQ=255)
aTCCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCcgctcgct                                                                                         <  2:4116988/93‑1 (MQ=255)
 tCCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCtt                                                                                        <  2:507364/93‑1 (MQ=255)
  ccGCGATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTa                                                                                       >  1:3653201/1‑93 (MQ=255)
  ccGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTa                                                                                       <  2:2089340/93‑1 (MQ=255)
  ccGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTa                                                                                       <  2:3991465/93‑1 (MQ=255)
   cGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTAc                                                                                      >  2:1376269/1‑93 (MQ=255)
     cAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACgg                                                                                    <  2:662806/93‑1 (MQ=255)
      aaTGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCTCGAAGGAAGGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACggc                                                                                   >  1:3320664/1‑92 (MQ=255)
      aaTGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGt                                                                                   <  2:3526424/93‑1 (MQ=255)
      aaTGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGt                                                                                   >  1:3081996/1‑93 (MQ=255)
       aTGGCTCTACGGCGATTGCTCGAATTATGGCGGGAAGGTTCTTGCGCTAAGGACTGCAACAAATCGTCGTCAGACTGCCGCTCGTTTACAGTc                                                                                  >  2:3219326/1‑93 (MQ=255)
            tctACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATcc                                                                             >  2:4117718/1‑93 (MQ=255)
            tctACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATcc                                                                             <  1:1750268/93‑1 (MQ=255)
              tACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCAt                                                                           <  1:1758771/93‑1 (MQ=255)
               aCGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCAtt                                                                          <  1:3086708/93‑1 (MQ=255)
               aCGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCAtt                                                                          >  1:4610278/1‑93 (MQ=255)
                cGGCGATTGCTGGAATTATGGCGGAAGGTTTCCTGCGCTAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTa                                                                         <  1:2339710/93‑1 (MQ=255)
                cGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTa                                                                         <  2:2256702/93‑1 (MQ=255)
                cGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTa                                                                         <  2:334914/93‑1 (MQ=255)
                cGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTa                                                                         <  2:2256703/93‑1 (MQ=255)
                cGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTa                                                                         <  1:3559963/93‑1 (MQ=255)
                  gCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTAcc                                                                       <  1:577901/93‑1 (MQ=255)
                      ttGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGAcc                                                                   <  1:1453063/93‑1 (MQ=255)
                      ttGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGAcc                                                                   >  2:2078281/1‑93 (MQ=255)
                       tGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGAccc                                                                  <  1:3423570/93‑1 (MQ=255)
                        gcgtgAATTAGTGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCg                                                                 <  1:2953219/89‑1 (MQ=255)
                         gtGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCgt                                                                <  1:1521606/92‑1 (MQ=255)
                           ggAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCgtgt                                                              <  1:1604670/93‑1 (MQ=255)
                             aaTTATGACGAAAAGTTCCCTGCGCGAAGCAACGCAAAAAACCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTct                                                            <  2:1805317/93‑1 (MQ=255)
                              aTTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctg                                                           <  2:3954361/93‑1 (MQ=255)
                              aTTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctg                                                           <  2:2588791/93‑1 (MQ=255)
                              aTTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctg                                                           <  2:2194972/93‑1 (MQ=255)
                              aTTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctg                                                           >  1:1204834/1‑93 (MQ=255)
                              aTTATGGCGGAAAGTTTCCTGCGCGAAAGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctg                                                           <  2:3954360/93‑1 (MQ=255)
                               ttttGGGGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgg                                                          <  2:278132/90‑1 (MQ=255)
                               ttATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgg                                                          <  1:3942606/93‑1 (MQ=255)
                               ttATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgg                                                          <  1:2360504/93‑1 (MQ=255)
                                tATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctggg                                                         <  1:492488/93‑1 (MQ=255)
                                tATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctggg                                                         <  1:893493/93‑1 (MQ=255)
                                 aTGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgggc                                                        <  2:2682165/93‑1 (MQ=255)
                                 aTGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgggc                                                        <  2:577827/93‑1 (MQ=255)
                                 aTGGCGGAAAGTTGCCTGCGCAAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgggc                                                        <  2:1042342/93‑1 (MQ=255)
                                  tGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgggct                                                       <  1:4738074/93‑1 (MQ=255)
                                  tGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgggct                                                       <  1:3420434/93‑1 (MQ=255)
                                    gCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgggctgg                                                     <  2:4522396/93‑1 (MQ=255)
                                     cGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgggctggg                                                    <  2:1952001/93‑1 (MQ=255)
                                     cGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgggctggg                                                    >  2:2490730/1‑93 (MQ=255)
                                      ggAAAGTTTCCTGCGCGAAGGAGCGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGgga                                                   <  2:930631/93‑1 (MQ=255)
                                      ggAAACTTTCCTGCGCGAAGGAGCGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGgga                                                   <  2:930630/93‑1 (MQ=255)
                                       gAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGggag                                                  >  2:4732046/1‑93 (MQ=255)
                                       gAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGggag                                                  <  2:710136/93‑1 (MQ=255)
                                       gAAAGTTTCCTGCGCGAAGGAACGCAAAAAACCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGggag                                                  <  1:256801/93‑1 (MQ=255)
                                         aaGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGggagga                                                <  2:1848623/93‑1 (MQ=255)
                                          aGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGAt                                               <  1:1510467/93‑1 (MQ=255)
                                           gTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATa                                              >  2:3194728/1‑93 (MQ=255)
                                            tttCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATAc                                             <  2:671686/93‑1 (MQ=255)
                                            tttCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATAc                                             <  1:42942/93‑1 (MQ=255)
                                            tttCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATAc                                             <  2:935945/93‑1 (MQ=255)
                                            tttCCTGCGCGAAGGAACGCAAAAAATCGGCGGCAGTCAGCCGCTCGCTTACGGTCAATCCTTTACCGACCCGTGTCTGGGCTGGGAGGATAc                                             <  2:3743367/93‑1 (MQ=255)
                                             ttCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATAcc                                            >  2:1261505/1‑93 (MQ=255)
                                             ttCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATAcc                                            <  1:1380154/93‑1 (MQ=255)
                                              tCCTGCGCGAAAGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCg                                           <  1:1236629/93‑1 (MQ=255)
                                               ccTGCGCGAAGGAACGCAAAAAATCGTAGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGa                                          <  1:3896737/93‑1 (MQ=255)
                                                cTGCGCGAAGGACCGCAAAAATTCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGGCCCGTGTCTGGGCTGGGAGGATCCCGaa                                         <  2:3798406/93‑1 (MQ=255)
                                                cTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGaa                                         >  1:2898199/1‑93 (MQ=255)
                                                cTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGaa                                         <  1:2903090/93‑1 (MQ=255)
                                                 cgcgcGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAAc                                        <  2:954428/92‑1 (MQ=255)
                                                    gcgAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGcc                                     <  1:377765/93‑1 (MQ=255)
                                                    gcgAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGcc                                     <  1:1716500/93‑1 (MQ=255)
                                                     cgAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCt                                    >  2:4582033/1‑93 (MQ=255)
                                                     cgAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCt                                    >  2:4643432/1‑93 (MQ=255)
                                                      gAAGTAACGCAAAAAATCGTCGGCAGTCAGCCGCTCTCTTACTGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTg                                   <  2:4070992/93‑1 (MQ=255)
                                                      gAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTCCCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTg                                   <  2:2375088/93‑1 (MQ=255)
                                                      gAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTg                                   <  2:3680953/93‑1 (MQ=255)
                                                      gAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTg                                   <  1:3071653/93‑1 (MQ=255)
                                                      gAAGGAACGCAAAAAATCCGCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTg                                   <  1:72043/93‑1 (MQ=255)
                                                        aGGAACGCACAAAACCGTCTACCCCCAGCCGCTCTCTTACGGTCAATCCATTAACGACCCGTGTCTGGGCTGGGAGTATACAGAACGCCTGGt                                 <  2:1627885/93‑1 (MQ=255)
                                                        aGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGt                                 >  1:808446/1‑93 (MQ=255)
                                                        aGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGt                                 <  1:1361786/93‑1 (MQ=255)
                                                         ggAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTc                                <  2:469672/93‑1 (MQ=255)
                                                            aCGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGaa                             >  1:2282446/1‑93 (MQ=255)
                                                              gCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGaaaa                           >  2:3430007/1‑93 (MQ=255)
                                                               cAAAAAACCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGaaaaa                          <  2:1285096/93‑1 (MQ=255)
                                                                 aaaaaTCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACt                        >  2:3045235/1‑93 (MQ=255)
                                                                  aaaaTCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTc                       <  1:900583/93‑1 (MQ=255)
                                                                     aTCGTCTGCAGGCACCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGcc                    <  2:3437833/93‑1 (MQ=255)
                                                                     aTCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGcc                    <  1:435719/93‑1 (MQ=255)
                                                                       cgtcgGCAGTCAGCCTCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCctc                  <  1:1080359/93‑1 (MQ=255)
                                                                       cgtcgGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCctc                  <  2:516759/93‑1 (MQ=255)
                                                                       cgtcgGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCctc                  <  2:1559953/93‑1 (MQ=255)
                                                                        gtcgGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCctct                 <  2:2381508/93‑1 (MQ=255)
                                                                          cgGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGc               <  1:2635498/93‑1 (MQ=255)
                                                                           ggCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCg              <  1:147824/93‑1 (MQ=255)
                                                                           ggCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCg              <  2:4621879/93‑1 (MQ=255)
                                                                           ggCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCg              <  2:3423510/93‑1 (MQ=255)
                                                                                 cAGCCTCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGAt        <  1:1806147/93‑1 (MQ=255)
                                                                                   gCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATAc      <  1:840206/93‑1 (MQ=255)
                                                                                     cgctcgctTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATAccc    <  1:438207/93‑1 (MQ=255)
                                                                                       cTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGc  <  2:1126995/93‑1 (MQ=255)
                                                                                          |                                                                                         
ATCCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGC  >  NC_000913/1787320‑1787499

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AATCCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTT  >  NC_000913/1787319‑1787501
                                                                                           |                                                                                           
AATCCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCT                                                                                           >  SRR4302544.793196/1‑94 (MQ=60)
 ATCCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTT                                                                                          >  SRR4302544.3180986/1‑94 (MQ=60)
 ATCCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTT                                                                                          <  SRR4302544.3421773/94‑1 (MQ=60)
  TCCGCGATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTA                                                                                         >  SRR4302544.3689652/1‑94 (MQ=60)
      CAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCTCGAAGGAAGGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGC                                                                                     >  SRR4302544.3352299/1‑94 (MQ=60)
      CAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGT                                                                                     >  SRR4302544.3111355/1‑94 (MQ=60)
            CTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCC                                                                               >  SRR4302544.1557401/1‑94 (MQ=60)
             TCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCA                                                                              <  SRR4302544.1767398/94‑1 (MQ=60)
               TACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATT                                                                            >  SRR4302544.4656066/1‑94 (MQ=60)
               TACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATT                                                                            <  SRR4302544.1775983/94‑1 (MQ=60)
                ACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTA                                                                           <  SRR4302544.3116116/94‑1 (MQ=60)
                 CGGCGATTGCTGGAATTATGGCGGAAGGTTTCCTGCGCTAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTAC                                                                          <  SRR4302544.2361875/94‑1 (MQ=60)
                 CGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTAC                                                                          <  SRR4302544.3594092/94‑1 (MQ=60)
                   GCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCG                                                                        <  SRR4302544.583548/94‑1 (MQ=60)
                       TTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCC                                                                    <  SRR4302544.1467499/94‑1 (MQ=60)
                        TGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCG                                                                   <  SRR4302544.3456301/94‑1 (MQ=60)
                        TGCTGATATTATGGCCAAAAGTTTCCAGCGCGAGGAAACACAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCG                                                                   <  SRR4302544.3524962/94‑1 (MQ=60)
                         gcgtgaattagtGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGT                                                                  <  SRR4302544.2981586/82‑1 (MQ=60)
                          GTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTG                                                                 <  SRR4302544.1536705/94‑1 (MQ=60)
                            GGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTC                                                               <  SRR4302544.1620495/94‑1 (MQ=60)
                              AATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTG                                                             >  SRR4302544.1216771/1‑94 (MQ=60)
                                TTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGG                                                           <  SRR4302544.2382853/94‑1 (MQ=60)
                                TTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGG                                                           <  SRR4302544.3981776/94‑1 (MQ=60)
                                 TATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGC                                                          <  SRR4302544.902335/94‑1 (MQ=60)
                                 TATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGC                                                          <  SRR4302544.497331/94‑1 (MQ=60)
                                 TATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGC                                                          >  SRR4302544.4201532/1‑94 (MQ=60)
                                   TGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTG                                                        <  SRR4302544.4784775/94‑1 (MQ=60)
                                   TGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTG                                                        <  SRR4302544.3453140/94‑1 (MQ=60)
                                        GAAAGTTTCCTGCGCGAAGGAACGCAAAAAACCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGG                                                   <  SRR4302544.259525/94‑1 (MQ=60)
                                           AGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATA                                                <  SRR4302544.1525440/94‑1 (MQ=60)
                                             TTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACC                                              <  SRR4302544.43363/94‑1 (MQ=60)
                                              TTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCG                                             <  SRR4302544.1393934/94‑1 (MQ=60)
                                               TCCTGCGCGAAAGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGA                                            <  SRR4302544.1248854/94‑1 (MQ=60)
                                                CCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAA                                           >  SRR4302544.2926114/1‑94 (MQ=60)
                                                CCTGCGCGAAGGAACGCAAAAAATCGTAGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAA                                           <  SRR4302544.3935439/94‑1 (MQ=60)
                                                 CTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAAC                                          <  SRR4302544.2931045/94‑1 (MQ=60)
                                                     GCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCT                                      <  SRR4302544.381617/94‑1 (MQ=60)
                                                     GCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCT                                      <  SRR4302544.1733375/94‑1 (MQ=60)
                                                       GAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGG                                    <  SRR4302544.3100926/94‑1 (MQ=60)
                                                       GAAGGAACGCAAAAAATCCGCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGG                                    <  SRR4302544.72740/94‑1 (MQ=60)
                                                        AAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGT                                   >  SRR4302544.816333/1‑94 (MQ=60)
                                                         AGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTC                                  <  SRR4302544.1375308/94‑1 (MQ=60)
                                                            AACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAA                               >  SRR4302544.2303951/1‑94 (MQ=60)
                                                                   AAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCG                        <  SRR4302544.909494/94‑1 (MQ=60)
                                                                      ATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCT                     <  SRR4302544.440108/94‑1 (MQ=60)
                                                                        CGTCGGCAGTCAGCCTCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCT                   <  SRR4302544.1091098/94‑1 (MQ=60)
                                                                           CGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCG                <  SRR4302544.2660754/94‑1 (MQ=60)
                                                                            GGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGG               <  SRR4302544.149236/94‑1 (MQ=60)
                                                                                 TCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGAT          <  SRR4302544.24745/94‑1 (MQ=60)
                                                                                  CAGCCTCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATA         <  SRR4302544.1823799/94‑1 (MQ=60)
                                                                                    GCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACC       <  SRR4302544.848500/94‑1 (MQ=60)
                                                                                      CGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCG     <  SRR4302544.442610/94‑1 (MQ=60)
                                                                                         TCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTT  >  SRR4302544.281054/1‑94 (MQ=60)
                                                                                           |                                                                                           
AATCCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTT  >  NC_000913/1787319‑1787501

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