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,370,890 T→C V105A (GTG→GCG)  b4533 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,370,8900TC99.2% 415.7 / ‑4.2 129V105A (GTG→GCG) b4533hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (54/74);  minor base G (0/1);  total (54/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.81e-01

Only 100 of 137 total aligned reads displayed.

TCTGGGTGACGCTGGCTGCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGCGTTCATTATTGGCGGCATTGTGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAGCCTGGGTTTTGC  >  NC_000913/2370800‑2370979
                                                                                          |                                                                                         
  tGGGTGACGCTGGCTGCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATc                                                                                       <  2:3599609/93‑1 (MQ=255)
   gggTGACGCTGGCTGCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATcc                                                                                      >  1:4107557/1‑93 (MQ=255)
    ggTGACGCTGGCAGCAGTACAACTGTGACACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTAGCGGCATTGCGATCCt                                                                                     <  2:4209433/93‑1 (MQ=255)
     gTGACGCTGGCTGCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTc                                                                                    <  2:224888/93‑1 (MQ=255)
     gTGACGCTGGCTGCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTc                                                                                    >  2:1741544/1‑93 (MQ=255)
      tGACGCTGGCTGCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCg                                                                                   <  2:2015644/93‑1 (MQ=255)
        aCGCTGGCTGCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCggg                                                                                 <  2:2056179/93‑1 (MQ=255)
         cGCTGGCTGCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGa                                                                                >  2:3968958/1‑93 (MQ=255)
         cGCTGGCTGCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGa                                                                                <  1:2153897/93‑1 (MQ=255)
         cGCTGGCTGCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGa                                                                                <  1:1531395/93‑1 (MQ=255)
          gctggctgCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAg                                                                               >  2:4047390/1‑93 (MQ=255)
           ctggctgCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGt                                                                              <  1:2053752/93‑1 (MQ=255)
            tggctgCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTa                                                                             >  1:3537238/1‑93 (MQ=255)
            tggctgCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGAGGTGGTGTTCATTATTGGCGGTATTGCGCTCCTCGGGAgaa                                                                             >  1:2011526/1‑91 (MQ=255)
            tggctgCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGCGTGGGGTGGTGTTCATTATTGGCGGCCTTGCGATCCTCGGGAGtt                                                                             >  1:1205288/1‑92 (MQ=255)
             ggctgCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGATTAc                                                                            <  1:4562601/93‑1 (MQ=255)
              gctgCAGTAAAACTGTGGCACGAACCTGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACg                                                                           <  2:825446/93‑1 (MQ=255)
              gctgCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACg                                                                           <  2:3479898/93‑1 (MQ=255)
               ctgCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACgg                                                                          <  2:3158989/93‑1 (MQ=255)
               ctgCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACgg                                                                          >  1:608481/1‑93 (MQ=255)
                 gCATTAAAACTGTGCCACGACCCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGgtg                                                                        <  2:3172565/93‑1 (MQ=255)
                 gCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCCTTATTGGCGGCATTGTGTTCCTCGGGGGCACGgtg                                                                        >  2:1305102/1‑93 (MQ=255)
                 gCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGgtg                                                                        >  2:727510/1‑93 (MQ=255)
                 gCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGgtg                                                                        <  2:3172566/93‑1 (MQ=255)
                   aGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTa                                                                      <  1:778430/93‑1 (MQ=255)
                   aGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTa                                                                      >  1:261193/1‑93 (MQ=255)
                     tAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAAt                                                                    >  2:599212/1‑93 (MQ=255)
                      aaaaCTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATg                                                                   <  1:3515920/93‑1 (MQ=255)
                      aaaaCTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATg                                                                   <  2:244558/93‑1 (MQ=255)
                             tgGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGa                                                            <  1:2556612/93‑1 (MQ=255)
                              ggCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGAt                                                           >  1:2768961/1‑93 (MQ=255)
                              ggCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGAt                                                           <  2:3070692/93‑1 (MQ=255)
                                cACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGAtgt                                                         <  2:3472387/93‑1 (MQ=255)
                                 aCGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGAtgtg                                                        >  2:1033902/1‑93 (MQ=255)
                                  cGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTgg                                                       <  1:3565452/93‑1 (MQ=255)
                                    aaCCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTgggg                                                     >  1:1735481/1‑93 (MQ=255)
                                     aCCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGa                                                    <  1:3023692/93‑1 (MQ=255)
                                       cGGTATCGCCGCGTCCCTCGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGatt                                                  <  2:1252316/93‑1 (MQ=255)
                                       cGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTAGGGTGTAATGGGCCTGCTGTGGGGgtt                                                  >  1:1166762/1‑93 (MQ=255)
                                       cGGTATCCCCCGGTCAGTGGTGTCGGGTCGTTTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGatt                                                  <  2:1005621/93‑1 (MQ=255)
                                        ggTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGatta                                                 <  2:1868037/93‑1 (MQ=255)
                                        ggTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGatta                                                 >  1:4404988/1‑93 (MQ=255)
                                        ggTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGatta                                                 >  1:2250464/1‑93 (MQ=255)
                                        ggTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTAGGGAGTACGGTGTAATGGGCCTGATGTGGGGatta                                                 <  2:2751100/93‑1 (MQ=255)
                                           aTCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTc                                              <  2:3075285/93‑1 (MQ=255)
                                           aTCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTc                                              <  2:903879/93‑1 (MQ=255)
                                            tCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCa                                             <  1:477788/93‑1 (MQ=255)
                                            tCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCa                                             <  1:3333087/93‑1 (MQ=255)
                                               ccgcGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCg                                          >  1:4368575/1‑93 (MQ=255)
                                                cgcgTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGt                                         >  1:4019681/1‑93 (MQ=255)
                                                 gcgTCACTGGTGTGGGGTTGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTc                                        >  2:4151356/1‑93 (MQ=255)
                                                  cgTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCa                                       <  1:2461864/93‑1 (MQ=255)
                                                  cgTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCa                                       <  2:2391262/93‑1 (MQ=255)
                                                   gTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCat                                      <  2:1507878/93‑1 (MQ=255)
                                                   gTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCat                                      >  1:3447823/1‑93 (MQ=255)
                                                    tCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCatt                                     >  2:4476059/1‑93 (MQ=255)
                                                     cACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCatta                                    >  2:3259949/1‑93 (MQ=255)
                                                     cACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCatta                                    >  2:1704941/1‑93 (MQ=255)
                                                     cACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCatta                                    >  2:559750/1‑93 (MQ=255)
                                                      aCTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCattat                                   >  2:3328954/1‑93 (MQ=255)
                                                      aCTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCattat                                   >  2:474701/1‑93 (MQ=255)
                                                      aCTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCattat                                   <  1:3926099/93‑1 (MQ=255)
                                                      aCTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCattat                                   <  1:798855/93‑1 (MQ=255)
                                                       cTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCattatt                                  <  1:1835294/93‑1 (MQ=255)
                                                        tGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGTATTATTCAGCGTCATTATTg                                 <  2:4520705/93‑1 (MQ=255)
                                                           tgtgGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCt                              <  1:1777776/93‑1 (MQ=255)
                                                           ggtgGGGTGGTGGTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCt                              <  1:1777052/92‑1 (MQ=255)
                                                             tgGGGTGGTGTTCATTATTGGCGGCATTACGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCa                            <  1:3593991/93‑1 (MQ=255)
                                                              ggggTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAg                           >  2:4131763/1‑93 (MQ=255)
                                                                 gtggtgTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTg                        <  2:1662146/93‑1 (MQ=255)
                                                                      gtcaattattGGCGGGATTGCGATGCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGAGTTATTCAGCGTCATTATTGCCTCAGTTGCGCaa                   <  1:3134778/89‑1 (MQ=255)
                                                                      gTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCaa                   <  1:1157336/93‑1 (MQ=255)
                                                                       ttCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGGAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAAt                  >  2:2809372/1‑93 (MQ=255)
                                                                         cATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTa                >  2:189647/1‑93 (MQ=255)
                                                                         cATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTa                >  2:2514705/1‑93 (MQ=255)
                                                                          attattGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTaa               >  2:4406406/1‑93 (MQ=255)
                                                                          attattGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTaa               <  2:2685618/93‑1 (MQ=255)
                                                                          attattGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTaa               >  1:3698435/1‑93 (MQ=255)
                                                                           ttattGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAg              <  1:3142955/93‑1 (MQ=255)
                                                                           ttattGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAg              >  2:879948/1‑93 (MQ=255)
                                                                           ttattGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAg              >  1:4065513/1‑93 (MQ=255)
                                                                           ttattGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAg              >  2:4554349/1‑93 (MQ=255)
                                                                           ttattGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAg              >  2:3819900/1‑93 (MQ=255)
                                                                           ttattGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAg              >  1:3113710/1‑93 (MQ=255)
                                                                            tattGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAGc             <  2:3850684/93‑1 (MQ=255)
                                                                             attGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAGcc            <  2:2537112/93‑1 (MQ=255)
                                                                              ttGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAGCCt           <  1:3755688/93‑1 (MQ=255)
                                                                              ttGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAGCCt           <  1:3565523/93‑1 (MQ=255)
                                                                              ttGGCGGCATTGCGATACTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAGCCt           <  1:4645213/93‑1 (MQ=255)
                                                                               tGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAGCCTg          <  1:2987096/93‑1 (MQ=255)
                                                                               tGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAGCCTg          >  2:2143909/1‑93 (MQ=255)
                                                                                ggcggcATTACAATCCTCGGGAGTACAGTGTAAAGGGCCTGATGTGGGGATTATTCAGCGTCATTAATGCCTCAGTTGCGCAATTAAGCCtag         <  2:4200744/93‑3 (MQ=255)
                                                                                  cggcATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAGCCTGGGt       <  2:407138/93‑1 (MQ=255)
                                                                                   ggcATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAGCCTGGGtt      >  1:3647066/1‑93 (MQ=255)
                                                                                   ggcATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAGCCTGGGtt      <  2:3628585/93‑1 (MQ=255)
                                                                                    gcATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAGCCTGGGttt     <  2:41379/93‑1 (MQ=255)
                                                                                      aTTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAGCCTGGGTTTTg   >  2:69822/1‑93 (MQ=255)
                                                                                      aTTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAGCCTGGGTTTTg   >  2:3372742/1‑93 (MQ=255)
                                                                                      aTTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAGCCTGGGTTTTg   >  2:456397/1‑93 (MQ=255)
                                                                                       ttGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAGCCTGGGTTTTgc  >  1:1680984/1‑93 (MQ=255)
                                                                                          |                                                                                         
TCTGGGTGACGCTGGCTGCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGCGTTCATTATTGGCGGCATTGTGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAGCCTGGGTTTTGC  >  NC_000913/2370800‑2370979

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ 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
TGTCTGGGTGACGCTGGCTGCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGCGTTCATTATTGGCGGCATTGTGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAGCCTGGGTTTTGCGGCG  >  NC_000913/2370798‑2370983
                                                                                            |                                                                                             
TGTCTGGGTGACGCTGGCTGCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCG                                                                                              >  SRR4302544.3728525/1‑94 (MQ=60)
  TCTGGGTGACGCTGGCTGCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGAT                                                                                            <  SRR4302544.927347/94‑1 (MQ=60)
    TGGGTGACGCTGGCTGCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCC                                                                                          >  SRR4302544.4148396/1‑94 (MQ=60)
           CGCTGGCTGCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAG                                                                                   <  SRR4302544.2174381/94‑1 (MQ=60)
           CGCTGGCTGCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAG                                                                                   <  SRR4302544.1546598/94‑1 (MQ=60)
             CTGGCTGCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTA                                                                                 >  SRR4302544.3571079/1‑94 (MQ=60)
             CTGGCTGCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTA                                                                                 <  SRR4302544.2073517/94‑1 (MQ=60)
             CTGGCTGCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGAGGTGGTGTTCATTATTGGCGGTATTGCGCTCCTCGGGAGAA                                                                                 >  SRR4302544.2031074/1‑94 (MQ=60)
             CTGGCTGCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGCGTGGGGTGGTGTTCATTATTGGCGGCCTTGCGATCCTCGGGAGTT                                                                                 >  SRR4302544.1217226/1‑94 (MQ=60)
              TGGCTGCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTAC                                                                                >  SRR4302544.3718501/1‑94 (MQ=60)
               GGCTGCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGGGATCCTCGGGAGTACG                                                                               <  SRR4302544.2106900/94‑1 (MQ=60)
               GGCTGCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGATTACG                                                                               <  SRR4302544.4608094/94‑1 (MQ=60)
                GCTGCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGG                                                                              >  SRR4302544.614453/1‑94 (MQ=60)
                    CAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTA                                                                          >  SRR4302544.263954/1‑94 (MQ=60)
                     AGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAA                                                                         <  SRR4302544.785979/94‑1 (MQ=60)
                      GTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAAT                                                                        <  SRR4302544.3077556/94‑1 (MQ=60)
                        AAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGG                                                                      <  SRR4302544.3549483/94‑1 (MQ=60)
                            CTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCT                                                                  <  SRR4302544.288234/94‑1 (MQ=60)
                               TGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGAT                                                               <  SRR4302544.2581158/94‑1 (MQ=60)
                               TGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGAT                                                               >  SRR4302544.2795603/1‑94 (MQ=60)
                                    CGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGG                                                          <  SRR4302544.3599633/94‑1 (MQ=60)
                                     GAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGG                                                         >  SRR4302544.1752487/1‑94 (MQ=60)
                                       ACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGAT                                                       <  SRR4302544.2547285/94‑1 (MQ=60)
                                       ACCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGAT                                                       <  SRR4302544.3052665/94‑1 (MQ=60)
                                        CCGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTAGGGTGTAATGGGCCTGCTGTGGGGGTT                                                      >  SRR4302544.1178396/1‑94 (MQ=60)
                                         CGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTA                                                     >  SRR4302544.4449480/1‑94 (MQ=60)
                                         CGGTATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTA                                                     >  SRR4302544.2271612/1‑94 (MQ=60)
                                            TATCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTC                                                  <  SRR4302544.1334128/94‑1 (MQ=60)
                                              TCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAG                                                <  SRR4302544.482465/94‑1 (MQ=60)
                                              TCGCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAG                                                <  SRR4302544.3364833/94‑1 (MQ=60)
                                                GCCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCG                                              >  SRR4302544.4412879/1‑94 (MQ=60)
                                                 CCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGT                                             >  SRR4302544.4059683/1‑94 (MQ=60)
                                                 CCGCGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGT                                             <  SRR4302544.4462190/94‑1 (MQ=60)
                                                    CGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCAT                                          <  SRR4302544.2485367/94‑1 (MQ=60)
                                                    CGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCAT                                          >  SRR4302544.3554418/1‑94 (MQ=60)
                                                    CGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCAT                                          >  SRR4302544.3480755/1‑94 (MQ=60)
                                                    CGTCACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCAT                                          >  SRR4302544.3341244/1‑94 (MQ=60)
                                                        ACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATT                                      <  SRR4302544.806600/94‑1 (MQ=60)
                                                        ACTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATT                                      <  SRR4302544.3965096/94‑1 (MQ=60)
                                                         CTGGTGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTG                                     <  SRR4302544.1853119/94‑1 (MQ=60)
                                                             TGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTC                                 <  SRR4302544.1794441/94‑1 (MQ=60)
                                                             TGTGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTC                                 <  SRR4302544.1795169/94‑1 (MQ=60)
                                                             GGTGGGGTGGTGGTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTC                                 <  SRR4302544.1794442/94‑1 (MQ=60)
                                                               TGGGGTGGTGTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAG                               <  SRR4302544.3628518/94‑1 (MQ=60)
                                                               TGGGGTGGTGTTCATTATTGGCGGCATTACGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAG                               <  SRR4302544.3628519/94‑1 (MQ=60)
                                                                        gtcaATTATTGGCGGGATTGCGATGCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGAGTTATTCAGCGTCATTATTGCCTCAGTTGCGCAAT                      <  SRR4302544.3164555/90‑1 (MQ=60)
                                                                        GTTCATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAAT                      <  SRR4302544.1168860/94‑1 (MQ=60)
                                                                           CATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAA                   >  SRR4302544.3735263/1‑94 (MQ=60)
                                                                            ATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAG                  >  SRR4302544.4105992/1‑94 (MQ=60)
                                                                            ATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAG                  >  SRR4302544.4015546/1‑94 (MQ=60)
                                                                            ATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAG                  <  SRR4302544.570695/94‑1 (MQ=60)
                                                                            ATTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAG                  >  SRR4302544.3143331/1‑94 (MQ=60)
                                                                             TTATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAGC                 <  SRR4302544.3172835/94‑1 (MQ=60)
                                                                               ATTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAGCCT               <  SRR4302544.197004/94‑1 (MQ=60)
                                                                                TTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAGCCTG              <  SRR4302544.3792968/94‑1 (MQ=60)
                                                                                TTGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAGCCTG              <  SRR4302544.3599706/94‑1 (MQ=60)
                                                                                TTGGCGGCATTGCGATACTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAGCCTG              <  SRR4302544.4691207/94‑1 (MQ=60)
                                                                                 TGGCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAGCCTGG             <  SRR4302544.3015749/94‑1 (MQ=60)
                                                                                   GCGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAGCCTGGGT           <  SRR4302544.4273251/94‑1 (MQ=60)
                                                                                    CGGCATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAGCCTGGGTT          >  SRR4302544.3683458/1‑94 (MQ=60)
                                                                                        ATTGCGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAGCCTGGGTTTTGC      >  SRR4302544.1697549/1‑94 (MQ=60)
                                                                                            CGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAGCCTGGGTTTTGCGGCG  >  SRR4302544.323369/1‑94 (MQ=60)
                                                                                            |                                                                                             
TGTCTGGGTGACGCTGGCTGCAGTAAAACTGTGGCACGAACCGGTATCGCCGCGTCACTGGTGTGGGGTGGCGTTCATTATTGGCGGCATTGTGATCCTCGGGAGTACGGTGTAATGGGCCTGATGTGGGGATTATTCAGCGTCATTATTGCCTCAGTTGCGCAATTAAGCCTGGGTTTTGCGGCG  >  NC_000913/2370798‑2370983

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