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 2,331,739 A→G intergenic (‑3922/+617) b2226 ← / ← b2229 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,331,7390AG99.1% 346.2 / ‑4.5 108intergenic (‑3922/+617)b2226/b2229hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (55/52);  total (56/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

Only 100 of 113 total aligned reads displayed.

GGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGC  >  NC_000913/2331649‑2331829
                                                                                          |                                                                                          
ggTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGCCCCGGCTCTTGTTGGCCGAGAGGAat                                                                                          >  1:1238452/1‑93 (MQ=255)
 gTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGata                                                                                         <  2:3628212/93‑1 (MQ=255)
  tATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGatat                                                                                        >  1:1311801/1‑93 (MQ=255)
     cGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATaaa                                                                                     <  1:68335/93‑1 (MQ=255)
      ggAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAAt                                                                                    >  2:2543139/1‑93 (MQ=255)
        aaacaaacACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGt                                                                                  <  1:1947405/93‑1 (MQ=255)
        aaacaaacACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGt                                                                                  <  2:1771322/93‑1 (MQ=255)
         aacaaacaCCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGtt                                                                                 <  1:4626305/93‑1 (MQ=255)
         aacaaacaCCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGtt                                                                                 <  2:3423015/93‑1 (MQ=255)
             aacaCCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTAccc                                                                             >  2:1300075/1‑93 (MQ=255)
             aacaCCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTAccc                                                                             <  2:4470761/93‑1 (MQ=255)
                cccaccGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGc                                                                          <  1:2950906/92‑1 (MQ=255)
                accaccGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGc                                                                          >  2:1586301/1‑93 (MQ=255)
                accaccGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGc                                                                          >  1:1262199/1‑93 (MQ=255)
                 ccaccGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCt                                                                         >  2:2207491/1‑93 (MQ=255)
                   accGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTgc                                                                       >  2:191702/1‑93 (MQ=255)
                     cgtcgtCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTgcgg                                                                     <  1:1619855/93‑1 (MQ=255)
                     cgtcgtCCCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTgcgg                                                                     <  1:3214953/93‑1 (MQ=255)
                      gtcgtcGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTgcggc                                                                    >  2:4049786/1‑93 (MQ=255)
                       ccgtcgCCCGATACCCACCACCCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTgcggcg                                                                   <  1:482450/92‑1 (MQ=255)
                        cgtcgCCCGATACCCACCAACCATCGCCTCGACGCGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCga                                                                  >  2:420065/1‑93 (MQ=255)
                        cgtcgCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCga                                                                  <  1:3485096/93‑1 (MQ=255)
                         gtcgCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCgag                                                                 >  2:1112066/1‑93 (MQ=255)
                           cgCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGAt                                                               <  2:1084930/93‑1 (MQ=255)
                            gCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATg                                                              >  1:170527/1‑93 (MQ=255)
                            gCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATg                                                              >  1:950973/1‑93 (MQ=255)
                              ccGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGaa                                                            <  2:4614518/93‑1 (MQ=255)
                              ccGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGaa                                                            >  1:248721/1‑93 (MQ=255)
                              ccGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGaa                                                            <  1:2706421/93‑1 (MQ=255)
                              ccGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGaa                                                            >  2:3756936/1‑93 (MQ=255)
                              ccGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGaa                                                            >  2:425103/1‑93 (MQ=255)
                               cGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGaaa                                                           <  1:103965/93‑1 (MQ=255)
                                gATACCCACCTTCCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAAt                                                          <  2:3071121/93‑1 (MQ=255)
                                gATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAAt                                                          >  2:1303847/1‑93 (MQ=255)
                                gATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAAt                                                          >  2:4697850/1‑93 (MQ=255)
                                 aTACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAAtt                                                         >  1:3873432/1‑93 (MQ=255)
                                 aTACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGAGATAAATGTTACCCGGCTGCAGCGAGATGAAAtt                                                         >  2:2765528/1‑93 (MQ=255)
                                   aCCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGc                                                       <  1:820555/93‑1 (MQ=255)
                                    cccACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCt                                                      <  1:3440235/93‑1 (MQ=255)
                                     ccaccaACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTg                                                     <  1:2202925/93‑1 (MQ=255)
                                     ccaccaACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTg                                                     >  1:2383359/1‑93 (MQ=255)
                                      caccaACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTgg                                                    >  1:4424323/1‑93 (MQ=255)
                                      caccaACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTgg                                                    >  1:927310/1‑93 (MQ=255)
                                       cccaaccaTCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGa                                                   <  1:3730670/92‑1 (MQ=255)
                                       accaaccaTCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGa                                                   <  1:2864010/93‑1 (MQ=255)
                                         caaccaTCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATg                                                 >  2:3662386/1‑93 (MQ=255)
                                          aaccaTCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATgc                                                <  1:1510380/93‑1 (MQ=255)
                                           accaTCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATgcg                                               <  1:3391595/93‑1 (MQ=255)
                                           accaTCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATgcg                                               >  1:1815781/1‑93 (MQ=255)
                                            ccaTCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATgcgc                                              <  2:2641880/93‑1 (MQ=255)
                                             caTCTCCCCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGcc                                             <  1:2286298/93‑1 (MQ=255)
                                             caTCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGcc                                             <  2:4122560/93‑1 (MQ=255)
                                             caTCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGcc                                             >  2:2915359/1‑93 (MQ=255)
                                              ctCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCt                                            <  1:994806/92‑1 (MQ=255)
                                              aTCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCt                                            <  2:2017093/93‑1 (MQ=255)
                                               tCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCtt                                           >  1:246820/1‑93 (MQ=255)
                                                cGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTc                                          <  2:2454924/93‑1 (MQ=255)
                                                 gCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTcc                                         >  1:4438954/1‑93 (MQ=255)
                                                  ccTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCa                                        >  2:1155057/1‑93 (MQ=255)
                                                      gacgaGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGttt                                    >  2:637169/1‑93 (MQ=255)
                                                       acgaGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGtttt                                   >  2:2116065/1‑93 (MQ=255)
                                                       acgaGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGtttt                                   >  1:2080023/1‑93 (MQ=255)
                                                         gaGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTAc                                 <  2:1171189/93‑1 (MQ=255)
                                                         gaGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTAc                                 <  1:4486535/93‑1 (MQ=255)
                                                          aGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTAcc                                >  2:3310903/1‑93 (MQ=255)
                                                          aGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTAcc                                >  2:4330731/1‑93 (MQ=255)
                                                          aGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTAcc                                <  1:690272/93‑1 (MQ=255)
                                                            gTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGGAATTGCTGGATGCGCCTTCCAGTTTTAcccc                              >  1:2400536/1‑93 (MQ=255)
                                                             tACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCt                             <  1:1247627/93‑1 (MQ=255)
                                                             tACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCt                             <  1:3780051/93‑1 (MQ=255)
                                                              aCAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCtt                            <  1:1892166/93‑1 (MQ=255)
                                                              aCAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCtt                            <  2:2030837/93‑1 (MQ=255)
                                                               cAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGGGCCTTCCAGTTTTACCCCttc                           >  2:3927105/1‑92 (MQ=255)
                                                               cAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTtg                           <  1:4091703/93‑1 (MQ=255)
                                                                aGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTtgc                          <  2:3210881/93‑1 (MQ=255)
                                                                 gTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTtgct                         >  2:713942/1‑93 (MQ=255)
                                                                 gTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTtgct                         <  1:1567574/93‑1 (MQ=255)
                                                                 gTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTtgct                         >  1:2804431/1‑93 (MQ=255)
                                                                   cccGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTtgctgc                       <  2:2951499/93‑1 (MQ=255)
                                                                   cccGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTtgctgc                       <  1:2874146/93‑1 (MQ=255)
                                                                    ccGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCg                      >  1:2341205/1‑93 (MQ=255)
                                                                      ggCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGc                    >  2:1461809/1‑93 (MQ=255)
                                                                        ctctTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTc                  >  2:3934884/1‑93 (MQ=255)
                                                                         tctTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCg                 >  1:2914622/1‑93 (MQ=255)
                                                                         tctTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCg                 >  2:3450061/1‑93 (MQ=255)
                                                                         tctTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCg                 >  2:150885/1‑93 (MQ=255)
                                                                         tctTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCg                 >  2:2739223/1‑93 (MQ=255)
                                                                           ttgttgGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCgga               >  2:1797510/1‑91 (MQ=255)
                                                                           ttgttgGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGaa               <  1:1631095/93‑1 (MQ=255)
                                                                           ttgttgGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGaa               >  1:980654/1‑93 (MQ=255)
                                                                           ttgttgGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGaa               >  2:1910886/1‑93 (MQ=255)
                                                                             gttgGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACg             <  1:1968222/93‑1 (MQ=255)
                                                                               tgGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGt           >  1:1012106/1‑93 (MQ=255)
                                                                                   cGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAg       >  2:3494077/1‑93 (MQ=255)
                                                                                     agagGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGcc     <  1:4488661/93‑1 (MQ=255)
                                                                                      gagGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGGCt    >  2:3041913/1‑93 (MQ=255)
                                                                                      gagGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCt    >  1:3017623/1‑93 (MQ=255)
                                                                                      gagGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGGTTAGCCt    >  2:1334066/1‑93 (MQ=255)
                                                                                       agGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGGTTAGCCTg   >  1:3148217/1‑93 (MQ=255)
                                                                                        ggAATATAAATGTTACCCGGCTGCGGCGAGATGGAATTGCTGGATGCGCCTTCCAGTTATACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGc  >  2:2367282/1‑93 (MQ=255)
                                                                                          |                                                                                          
GGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGC  >  NC_000913/2331649‑2331829

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 28 ≤ 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
CACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCA  >  NC_000913/2331646‑2331831
                                                                                             |                                                                                            
CACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGG                                                                                              <  SRR4302544.4550001/94‑1 (MQ=60)
CACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGG                                                                                              >  SRR4302544.343898/1‑94 (MQ=60)
  CGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGCCCCGGCTCTTGTTGGCCGAGAGGAAT                                                                                            >  SRR4302544.1250688/1‑94 (MQ=60)
   GGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATA                                                                                           >  SRR4302544.3185107/1‑94 (MQ=60)
    GTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATAT                                                                                          >  SRR4302544.1324700/1‑94 (MQ=60)
        CGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAAT                                                                                      >  SRR4302544.441611/1‑94 (MQ=60)
        CGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAAT                                                                                      <  SRR4302544.69010/94‑1 (MQ=60)
           AAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTT                                                                                   <  SRR4302544.1966267/94‑1 (MQ=60)
            AACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTA                                                                                  <  SRR4302544.4672180/94‑1 (MQ=60)
                  CACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGC                                                                            >  SRR4302544.1274622/1‑94 (MQ=60)
                   CCCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCT                                                                           <  SRR4302544.2979254/94‑1 (MQ=60)
                        CGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGC                                                                      <  SRR4302544.1635814/94‑1 (MQ=60)
                        CGTCGTCCCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGC                                                                      <  SRR4302544.3245621/94‑1 (MQ=60)
                          CCGTCGCCCGATACCCACCACCCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGA                                                                    <  SRR4302544.487174/94‑1 (MQ=60)
                           CGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAG                                                                   <  SRR4302544.3518310/94‑1 (MQ=60)
                              CGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATG                                                                >  SRR4302544.172208/1‑94 (MQ=60)
                              CGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATG                                                                >  SRR4302544.960322/1‑94 (MQ=60)
                                CCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAA                                                              >  SRR4302544.251330/1‑94 (MQ=60)
                                 CCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAA                                                             <  SRR4302544.2732444/94‑1 (MQ=60)
                                  CGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAAT                                                            <  SRR4302544.104921/94‑1 (MQ=60)
                                   GATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATT                                                           >  SRR4302544.3911899/1‑94 (MQ=60)
                                   CATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATT                                                           <  SRR4302544.213752/94‑1 (MQ=60)
                                      ACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCT                                                        <  SRR4302544.828581/94‑1 (MQ=60)
                                       CCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTG                                                       <  SRR4302544.3473121/94‑1 (MQ=60)
                                       CCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTG                                                       >  SRR4302544.2405927/1‑94 (MQ=60)
                                        CCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGG                                                      >  SRR4302544.4468959/1‑94 (MQ=60)
                                        CCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGG                                                      <  SRR4302544.2223676/94‑1 (MQ=60)
                                        CCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGG                                                      >  SRR4302544.936448/1‑94 (MQ=60)
                                         CACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGA                                                     <  SRR4302544.769497/94‑1 (MQ=60)
                                          CCCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGAT                                                    <  SRR4302544.3767761/94‑1 (MQ=60)
                                          ACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGAT                                                    <  SRR4302544.2891656/94‑1 (MQ=60)
                                             AACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCG                                                 <  SRR4302544.1525352/94‑1 (MQ=60)
                                             AACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCG                                                 >  SRR4302544.1833503/1‑94 (MQ=60)
                                              ACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGC                                                <  SRR4302544.3424016/94‑1 (MQ=60)
                                                CATCTCCCCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCT                                              <  SRR4302544.2307848/94‑1 (MQ=60)
                                                 CTCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTT                                             <  SRR4302544.1004646/94‑1 (MQ=60)
                                                 ATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTT                                             >  SRR4302544.249406/1‑94 (MQ=60)
                                                   CGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCC                                           >  SRR4302544.4483708/1‑94 (MQ=60)
                                                         GACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTT                                     >  SRR4302544.2100013/1‑94 (MQ=60)
                                                            GAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACC                                  <  SRR4302544.4531602/94‑1 (MQ=60)
                                                             AGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCC                                 <  SRR4302544.697151/94‑1 (MQ=60)
                                                              GGTACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGGAATTGCTGGATGCGCCTTCCAGTTTTACCCC                                >  SRR4302544.2423314/1‑94 (MQ=60)
                                                                TACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTT                              <  SRR4302544.3817491/94‑1 (MQ=60)
                                                                TACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTT                              <  SRR4302544.1259915/94‑1 (MQ=60)
                                                                 ACAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTG                             <  SRR4302544.1910441/94‑1 (MQ=60)
                                                                  CAGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGC                            <  SRR4302544.4132397/94‑1 (MQ=60)
                                                                   AGTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCT                           >  SRR4302544.2831469/1‑94 (MQ=60)
                                                                    GTCCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTG                          <  SRR4302544.1583106/94‑1 (MQ=60)
                                                                      CCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCG                        >  SRR4302544.2363384/1‑94 (MQ=60)
                                                                      CCCGGCTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCG                        <  SRR4302544.2901887/94‑1 (MQ=60)
                                                                           CTCTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCG                   >  SRR4302544.2942649/1‑94 (MQ=60)
                                                                             CTTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAA                 >  SRR4302544.990266/1‑94 (MQ=60)
                                                                              TTGTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAAC                <  SRR4302544.1647148/94‑1 (MQ=60)
                                                                                GTTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGG              <  SRR4302544.1987338/94‑1 (MQ=60)
                                                                                 TTGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGT             >  SRR4302544.1022192/1‑94 (MQ=60)
                                                                                  TGGCCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGATTTACCCCTTGCTGCGGCTCccaacgct            >  SRR4302544.2572623/1‑86 (MQ=60)
                                                                                     CCGAGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAG         <  SRR4302544.1622201/94‑1 (MQ=60)
                                                                                        AGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCT      <  SRR4302544.4533738/94‑1 (MQ=60)
                                                                                        AGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCT      >  SRR4302544.3046553/1‑94 (MQ=60)
                                                                                        AGAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCT      >  SRR4302544.2740079/1‑94 (MQ=60)
                                                                                         GAGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGGTTAGCCTG     >  SRR4302544.3178149/1‑94 (MQ=60)
                                                                                          AGGGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGC    >  SRR4302544.2480405/1‑94 (MQ=60)
                                                                                            GGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCC  >  SRR4302544.2912101/1‑94 (MQ=60)
                                                                                            GGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCA  >  SRR4302544.429191/1‑94 (MQ=60)
                                                                                            GGATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCA  <  SRR4302544.4118965/94‑1 (MQ=60)
                                                                                             |                                                                                            
CACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCA  >  NC_000913/2331646‑2331831

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: