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 2,297,335 A→G G88G (GGT→GGC napH ← ferredoxin‑type protein: electron transfer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,297,3350AG99.1% 347.3 / ‑4.6 113G88G (GGT→GGCnapHferredoxin‑type protein: electron transfer
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (45/67);  total (46/67)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.40e-01

Only 100 of 129 total aligned reads displayed.

CGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGGTGTCGAATAA  >  NC_000913/2297245‑2297424
                                                                                          |                                                                                         
cGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATcaccgcgcc                                                                                         <  2:2392175/93‑1 (MQ=255)
 gCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATcaccgcgccg                                                                                        >  2:2472103/1‑93 (MQ=255)
  cTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCACGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCgg                                                                                       <  1:1311138/93‑1 (MQ=255)
   tAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGt                                                                                      >  1:4558915/1‑93 (MQ=255)
   tAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGACGGt                                                                                      >  2:3831931/1‑93 (MQ=255)
    aaGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTc                                                                                     >  2:2669389/1‑93 (MQ=255)
     aGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCa                                                                                    <  2:1315212/93‑1 (MQ=255)
      gTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCaa                                                                                   >  1:4277078/1‑93 (MQ=255)
       tCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAAc                                                                                  <  2:817279/93‑1 (MQ=255)
         ggTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGc                                                                                <  2:2396841/93‑1 (MQ=255)
         ggTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGc                                                                                <  2:1736608/93‑1 (MQ=255)
          gTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGcc                                                                               >  2:3438139/1‑93 (MQ=255)
          gTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGcc                                                                               >  1:3704991/1‑93 (MQ=255)
          gTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGcc                                                                               <  2:650724/93‑1 (MQ=255)
            aaTCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCAc                                                                             >  2:3239289/1‑93 (MQ=255)
             aTCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCAcc                                                                            >  1:2196224/1‑93 (MQ=255)
             aTCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCAcc                                                                            <  1:4420738/93‑1 (MQ=255)
                 ggTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTgg                                                                        <  2:1664656/93‑1 (MQ=255)
                  gTTCGGCGGGCAGACCCAGCTGCAAAATAAACGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTggc                                                                       >  2:1129777/1‑93 (MQ=255)
                    tCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTggcag                                                                     >  2:2909417/1‑93 (MQ=255)
                      aGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTggcaggc                                                                   <  1:426646/93‑1 (MQ=255)
                      aGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTggcaggc                                                                   <  1:2818095/93‑1 (MQ=255)
                       gCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTggcaggca                                                                  >  2:348149/1‑93 (MQ=255)
                         gggCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGa                                                                >  1:2251650/1‑93 (MQ=255)
                          ggCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGAt                                                               <  2:14782/93‑1 (MQ=255)
                              gACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGAcc                                                           >  2:1903599/1‑93 (MQ=255)
                                 ccAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACt                                                        >  2:3421134/1‑93 (MQ=255)
                                  cAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTg                                                       <  1:3764133/93‑1 (MQ=255)
                                   aGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTgg                                                      <  2:3750295/93‑1 (MQ=255)
                                    gCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGc                                                     >  2:4597343/1‑93 (MQ=255)
                                        cAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGa                                                 <  2:1367129/93‑1 (MQ=255)
                                        cAAAATAATCGCTTTCCCACCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGa                                                 <  1:2613490/93‑1 (MQ=255)
                                          aaaTAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACt                                               >  2:4402750/1‑93 (MQ=255)
                                           aataatCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACtt                                              <  2:3305980/93‑1 (MQ=255)
                                            ataatCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACttt                                             <  1:3730452/93‑1 (MQ=255)
                                            ataatCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACttt                                             <  2:1311438/93‑1 (MQ=255)
                                            ataatCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACttt                                             <  2:1654048/93‑1 (MQ=255)
                                            ataatCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACttt                                             <  2:3071792/93‑1 (MQ=255)
                                            ataatCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACttt                                             <  1:542354/93‑1 (MQ=255)
                                            ataatCGCTTCCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACttt                                             <  2:1774035/93‑1 (MQ=255)
                                             taatCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTg                                            >  1:4418927/1‑93 (MQ=255)
                                              aatCCCTGTCCCGCCAGGGCGTGGAGCACGGTCATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGc                                           <  2:3252970/93‑1 (MQ=255)
                                               atCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCa                                          <  2:2894244/93‑1 (MQ=255)
                                               atCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCa                                          >  1:3548287/1‑93 (MQ=255)
                                                tCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAg                                         >  1:3741784/1‑93 (MQ=255)
                                                 cGCTTTCCAGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGc                                        >  2:4575264/1‑93 (MQ=255)
                                                  ccTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCg                                       <  2:2153873/92‑1 (MQ=255)
                                                   cTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGt                                      <  1:3057679/93‑1 (MQ=255)
                                                     ttCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCGCCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGtca                                    <  1:1419781/93‑1 (MQ=255)
                                                     tccccGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGtca                                    <  1:3173865/91‑1 (MQ=255)
                                                      tCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGGCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGtcat                                   <  2:3632407/93‑1 (MQ=255)
                                                       cccGCCAGGGCGTAGAGCACGGTGATAATCGCCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCTtcatc                                  <  2:1335991/93‑1 (MQ=255)
                                                        ccGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGtcatca                                 <  2:2880735/93‑1 (MQ=255)
                                                        ccGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGtcatca                                 <  2:732598/93‑1 (MQ=255)
                                                         cGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAg                                <  2:2416241/93‑1 (MQ=255)
                                                          gCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGc                               >  2:3719454/1‑93 (MQ=255)
                                                           ccAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCg                              >  2:1310642/1‑93 (MQ=255)
                                                           ccAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCg                              >  1:4192188/1‑93 (MQ=255)
                                                           ccAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCg                              >  1:1548261/1‑93 (MQ=255)
                                                           ccAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCg                              >  1:2911270/1‑93 (MQ=255)
                                                            cAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCgg                             <  2:437430/93‑1 (MQ=255)
                                                             aGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGa                            <  1:1749208/93‑1 (MQ=255)
                                                              gggCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGAt                           <  2:1734309/93‑1 (MQ=255)
                                                              gggCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGAt                           <  2:2025031/93‑1 (MQ=255)
                                                              gggCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGAt                           <  1:1512897/93‑1 (MQ=255)
                                                                gCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCg                         <  2:1623389/93‑1 (MQ=255)
                                                                 cGTAGAGCCCGGTGATAATCCCCGCGGCGGTCAACGCCCCCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCgg                        <  2:1798233/93‑1 (MQ=255)
                                                                 cGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCgg                        <  1:247715/93‑1 (MQ=255)
                                                                   tAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTc                      >  2:2350255/1‑93 (MQ=255)
                                                                   tAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTc                      >  2:3888419/1‑93 (MQ=255)
                                                                    agagCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCa                     <  1:1806945/93‑1 (MQ=255)
                                                                    agagCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCa                     >  1:834751/1‑93 (MQ=255)
                                                                    agagCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCa                     <  2:756841/93‑1 (MQ=255)
                                                                     gagCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAg                    >  1:2583099/1‑93 (MQ=255)
                                                                         aCGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGa                >  1:2318753/1‑93 (MQ=255)
                                                                         aCGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGa                >  1:2586491/1‑93 (MQ=255)
                                                                          cGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGaa               >  1:1541390/1‑93 (MQ=255)
                                                                          cGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGaa               <  2:1622285/93‑1 (MQ=255)
                                                                          cGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGaa               <  2:3466590/93‑1 (MQ=255)
                                                                           ggTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATGAGCGGCTCGGTCAGCGGAAc              <  2:3830846/93‑1 (MQ=255)
                                                                           ggTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAAc              <  1:3455784/93‑1 (MQ=255)
                                                                            gTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACg             <  1:3506399/93‑1 (MQ=255)
                                                                             tGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACgg            >  2:4590075/1‑93 (MQ=255)
                                                                             tGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACgg            <  1:1970201/93‑1 (MQ=255)
                                                                             tGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACgg            <  2:4049078/93‑1 (MQ=255)
                                                                             tGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACgg            <  1:1200961/93‑1 (MQ=255)
                                                                             tGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACgg            >  2:659657/1‑93 (MQ=255)
                                                                              gATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGgt           >  1:327222/1‑93 (MQ=255)
                                                                               atcaTCCCCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGgtg          <  1:109402/90‑1 (MQ=255)
                                                                               aTAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGgtg          >  2:4571755/1‑93 (MQ=255)
                                                                                tAATCACCGCGCCGGTCAGCGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGgtgt         <  2:1781592/93‑1 (MQ=255)
                                                                                tAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGgtgt         >  2:149818/1‑93 (MQ=255)
                                                                                tAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGgtgt         <  2:554586/93‑1 (MQ=255)
                                                                                  aTCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGGTGTCg       <  1:2771779/93‑1 (MQ=255)
                                                                                    caccgcgccgGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGGTGTCGaa     <  1:1963860/93‑1 (MQ=255)
                                                                                     accgcgccgGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGGTGTCGAAt    >  2:4385317/1‑93 (MQ=255)
                                                                                     accgcgccgGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGGTGTCGAAt    >  1:3471328/1‑93 (MQ=255)
                                                                                     accgcgccgGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGGTGTCGAAt    <  2:2010140/93‑1 (MQ=255)
                                                                                      ccgcgccgGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGGTGTCGAATa   <  2:2865842/93‑1 (MQ=255)
                                                                                       cgcgccgGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGGTGTCGAATaa  <  1:671233/93‑1 (MQ=255)
                                                                                          |                                                                                         
CGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGGTGTCGAATAA  >  NC_000913/2297245‑2297424

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 27 ≤ 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
CTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGGTGTCGAATAACAG  >  NC_000913/2297247‑2297427
                                                                                        |                                                                                            
CTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGT                                                                                         >  SRR4302535.4590040/1‑94 (MQ=60)
CTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCACGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGT                                                                                         <  SRR4302535.1320509/94‑1 (MQ=60)
 TAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTC                                                                                        >  SRR4302535.3492831/1‑94 (MQ=60)
   AGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAA                                                                                      >  SRR4302535.4306351/1‑94 (MQ=60)
       GGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCC                                                                                  >  SRR4302535.3730650/1‑94 (MQ=60)
          AATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACC                                                                               <  SRR4302535.3722433/94‑1 (MQ=60)
          AATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACC                                                                               >  SRR4302535.2211876/1‑94 (MQ=60)
          AATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACC                                                                               >  SRR4302535.1655566/1‑94 (MQ=60)
           ATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCG                                                                              <  SRR4302535.4451006/94‑1 (MQ=60)
                GTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCA                                                                         <  SRR4302535.2436036/94‑1 (MQ=60)
                    AGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCA                                                                     <  SRR4302535.2838113/94‑1 (MQ=60)
                    AGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCA                                                                     <  SRR4302535.429617/94‑1 (MQ=60)
                      CGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGA                                                                   >  SRR4302535.2267684/1‑94 (MQ=60)
                              CCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACT                                                           <  SRR4302535.1495227/94‑1 (MQ=60)
                                CAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGG                                                         <  SRR4302535.3790136/94‑1 (MQ=60)
                                 AGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGC                                                        >  SRR4302535.2750751/1‑94 (MQ=60)
                                    TGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGacag                                                     >  SRR4302535.440237/1‑90 (MQ=60)
                                      CAAAATAATCGCTTTCCCACCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGAC                                                   <  SRR4302535.2632235/94‑1 (MQ=60)
                                          ATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTG                                               <  SRR4302535.3756285/94‑1 (MQ=60)
                                          ATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTG                                               <  SRR4302535.546143/94‑1 (MQ=60)
                                          ATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTG                                               >  SRR4302535.4449183/1‑94 (MQ=60)
                                            AATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCA                                             >  SRR4302535.3572920/1‑94 (MQ=60)
                                             ATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAG                                            >  SRR4302535.3767675/1‑94 (MQ=60)
                                                 CTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTC                                        <  SRR4302535.3079071/94‑1 (MQ=60)
                                                  TTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCA                                       <  SRR4302535.488911/94‑1 (MQ=60)
                                                   TTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCGCCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCAT                                      <  SRR4302535.1429919/94‑1 (MQ=60)
                                                   TCCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCAT                                      <  SRR4302535.3195832/94‑1 (MQ=60)
                                                        GCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCG                                 >  SRR4302535.4220876/1‑94 (MQ=60)
                                                        GCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCG                                 >  SRR4302535.1559414/1‑94 (MQ=60)
                                                        GCCAGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCG                                 >  SRR4302535.2931833/1‑94 (MQ=60)
                                                           AGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGAT                              >  SRR4302535.2170425/1‑94 (MQ=60)
                                                           AGGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGAT                              <  SRR4302535.1761860/94‑1 (MQ=60)
                                                            GGGCGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATC                             <  SRR4302535.1523796/94‑1 (MQ=60)
                                                               CGTAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGT                          <  SRR4302535.249414/94‑1 (MQ=60)
                                                                 TAGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCA                        >  SRR4302535.840604/1‑94 (MQ=60)
                                                                  AGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAG                       <  SRR4302535.3881773/94‑1 (MQ=60)
                                                                  AGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAG                       >  SRR4302535.2601658/1‑94 (MQ=60)
                                                                  AGAGCACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAG                       <  SRR4302535.1819999/94‑1 (MQ=60)
                                                                      CACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGA                   >  SRR4302535.2605066/1‑94 (MQ=60)
                                                                      CACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGA                   >  SRR4302535.2335351/1‑94 (MQ=60)
                                                                       ACGGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAA                  >  SRR4302535.1552504/1‑94 (MQ=60)
                                                                         GGTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACG                <  SRR4302535.3479804/94‑1 (MQ=60)
                                                                          GTGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGG               <  SRR4302535.3530768/94‑1 (MQ=60)
                                                                           TGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGGT              >  SRR4302535.329490/1‑94 (MQ=60)
                                                                           TGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGGT              <  SRR4302535.1209565/94‑1 (MQ=60)
                                                                           TGATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGGT              <  SRR4302535.1984394/94‑1 (MQ=60)
                                                                            GATAATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGGTG             <  SRR4302535.3032817/94‑1 (MQ=60)
                                                                             ATCATCCCCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGGTGT            <  SRR4302535.110177/94‑1 (MQ=60)
                                                                                ATCACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGGTGTCGA         <  SRR4302535.2791516/94‑1 (MQ=60)
                                                                                  CACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGGTGTCGAAT       <  SRR4302535.1978019/94‑1 (MQ=60)
                                                                                  CACCGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGGTGTCGAAT       >  SRR4302535.3495450/1‑94 (MQ=60)
                                                                                     CGCGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGGTGTCGAATAAC    <  SRR4302535.675919/94‑1 (MQ=60)
                                                                                       CGCCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGGTGTCGAATAACAG  <  SRR4302535.1982298/94‑1 (MQ=60)
                                                                                        |                                                                                            
CTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGGTGTCGAATAACAG  >  NC_000913/2297247‑2297427

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