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 4,477,041 Δ1 bp intergenic (‑151/‑269) b4256 ← / → yjgN hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,477,0410G.99.2% 562.6 / ‑3.4 124intergenic (‑151/‑269)b4256/yjgNhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (44/79);  minor base A (1/0);  total (45/79)
Fisher's exact test for biased strand distribution p-value = 3.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 127 total aligned reads displayed.

CGTCAAGATTCGGATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCGAGGATTATCCTGCAGCGTCAGGCGACGATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAA  >  NC_000913/4476953‑4477126
                                                                                        |                                                                                       
  tCAAGATTCGGATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTa                                                                                  >  2:4179155/1‑93 (MQ=255)
  tCAAGATTCGGATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTa                                                                                  >  1:4442894/1‑93 (MQ=255)
  tCAAGATTCGGATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTa                                                                                  >  2:1755103/1‑93 (MQ=255)
   cAAGATTCGGATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCTGCGGCTACCTGACGAATTACGCGGGCGATAGCAGGATTAt                                                                                 <  1:485688/93‑1 (MQ=255)
    aaGATTCGGATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATc                                                                                >  2:3487712/1‑93 (MQ=255)
     aGATTCGGATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATcc                                                                               >  2:402046/1‑93 (MQ=255)
     aGATTCGGATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATcc                                                                               <  1:4648659/93‑1 (MQ=255)
      gATTCGGATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCt                                                                              >  1:1431510/1‑93 (MQ=255)
        ttCGGATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGc                                                                            >  2:2888401/1‑93 (MQ=255)
           ggATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGc                                                                         >  2:2268902/1‑93 (MQ=255)
           ggATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGc                                                                         >  2:340816/1‑93 (MQ=255)
           cggtCGGCCGCGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGc                                                                         <  1:391565/90‑1 (MQ=255)
            gATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGCCGACTGACGCTAGCGATAGCAGGATTATCCTGCAGCg                                                                        >  1:1354890/1‑93 (MQ=255)
            gATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCg                                                                        >  1:1354889/1‑93 (MQ=255)
            gATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCg                                                                        <  1:529221/93‑1 (MQ=255)
               cGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCa                                                                     <  2:2782767/93‑1 (MQ=255)
                ggCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAg                                                                    <  1:3356668/93‑1 (MQ=255)
                 gCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATCCCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAgg                                                                   <  1:2552221/93‑1 (MQ=255)
                 gCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAgg                                                                   <  2:324560/93‑1 (MQ=255)
                 gCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAgg                                                                   <  1:4707250/93‑1 (MQ=255)
                  cGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGc                                                                  >  1:2071603/1‑93 (MQ=255)
                  cGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGc                                                                  <  2:3459996/93‑1 (MQ=255)
                  cGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGc                                                                  <  2:4070325/93‑1 (MQ=255)
                   gACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGcg                                                                 <  1:4268582/93‑1 (MQ=255)
                   gACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGcg                                                                 >  2:3484156/1‑93 (MQ=255)
                    aCGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGcga                                                                <  2:1051063/93‑1 (MQ=255)
                     cGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGcgac                                                               >  1:2208822/1‑93 (MQ=255)
                         gtAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGcgacgcat                                                          >  1:4182178/1‑90 (MQ=255)
                               tttATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAgg                                                    <  1:1910419/93‑1 (MQ=255)
                                ttATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGt                                                   <  2:2136153/93‑1 (MQ=255)
                                ttATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGt                                                   >  1:1529983/1‑93 (MQ=255)
                                  aTCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGa                                                 <  1:1867763/93‑1 (MQ=255)
                                   tCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGAt                                                <  1:4164086/93‑1 (MQ=255)
                                   tCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGAt                                                <  2:2071407/93‑1 (MQ=255)
                                   tCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGAt                                                <  2:1824470/93‑1 (MQ=255)
                                   tCAGCGGGAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGAt                                                >  2:4652619/1‑93 (MQ=255)
                                    cAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGAtt                                               <  2:3128903/93‑1 (MQ=255)
                                     aGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTg                                              <  2:3535618/93‑1 (MQ=255)
                                      gCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTgc                                             >  2:4666671/1‑93 (MQ=255)
                                        ggTAAGACCGTATTCGGCTGATACCTGACGAATCACGGGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTgcgg                                           <  2:3297325/93‑1 (MQ=255)
                                        ggTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTgcgg                                           <  2:191703/93‑1 (MQ=255)
                                         gTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTgcggc                                          >  1:3156267/1‑93 (MQ=255)
                                           aaGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGgcgc                                        >  2:1146126/1‑93 (MQ=255)
                                           aaGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGgcgc                                        <  2:2023581/93‑1 (MQ=255)
                                           aaGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGgcgc                                        <  1:3851070/93‑1 (MQ=255)
                                             gACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCaa                                      >  1:1616557/1‑93 (MQ=255)
                                              aCCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAAt                                     <  1:2910103/93‑1 (MQ=255)
                                              aCCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAAt                                     >  1:2617580/1‑93 (MQ=255)
                                              aCCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAAt                                     <  2:3561578/93‑1 (MQ=255)
                                               ccGTAAACGGCGGATACCTGACGAATGAAGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCTCacgt                                    >  2:3780247/1‑90 (MQ=255)
                                                 gTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGGAATGtt                                  >  2:676323/1‑93 (MQ=255)
                                                  tATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTa                                 <  2:163593/93‑1 (MQ=255)
                                                  tATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTa                                 >  2:4184879/1‑93 (MQ=255)
                                                   aTTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTAt                                <  2:2696679/93‑1 (MQ=255)
                                                    ttCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTAtt                               >  1:2908395/1‑93 (MQ=255)
                                                     tCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTc                              <  1:1143431/93‑1 (MQ=255)
                                                     tCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTc                              <  1:4500598/93‑1 (MQ=255)
                                                      cggcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCa                             <  1:1973876/93‑1 (MQ=255)
                                                      cggcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGCCGCAATGTTATTCa                             >  1:479304/1‑93 (MQ=255)
                                                       ggcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCAt                            >  1:2858586/1‑93 (MQ=255)
                                                       ggcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCAt                            <  2:2552922/93‑1 (MQ=255)
                                                        gcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATa                           <  1:4170091/93‑1 (MQ=255)
                                                        gcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATa                           <  2:2001457/93‑1 (MQ=255)
                                                         cggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAg                          <  2:2858586/93‑1 (MQ=255)
                                                          ggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTTTTCATAGt                         >  1:1402860/1‑93 (MQ=255)
                                                          ggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGt                         <  1:2402857/93‑1 (MQ=255)
                                                          ggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGt                         >  2:213860/1‑93 (MQ=255)
                                                          ggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGt                         >  2:3415916/1‑93 (MQ=255)
                                                           gATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGtt                        <  1:185470/93‑1 (MQ=255)
                                                           gATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGtt                        >  2:4057760/1‑93 (MQ=255)
                                                           gATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGtt                        >  1:3733446/1‑93 (MQ=255)
                                                           gATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGtt                        >  2:4057761/1‑93 (MQ=255)
                                                              aCCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTaa                     >  1:1238512/1‑93 (MQ=255)
                                                              aCCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTaa                     >  1:3401881/1‑93 (MQ=255)
                                                               ccTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAAc                    >  2:4685728/1‑93 (MQ=255)
                                                               ccTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAAc                    <  1:1164653/93‑1 (MQ=255)
                                                                cTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAAcc                   <  2:3914288/93‑1 (MQ=255)
                                                                cTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAAcc                   >  2:91799/1‑93 (MQ=255)
                                                                 tGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAAccc                  <  1:1046616/93‑1 (MQ=255)
                                                                 tGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAAccc                  <  1:3844330/93‑1 (MQ=255)
                                                                 tGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAAccc                  <  2:4283208/93‑1 (MQ=255)
                                                                  gACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCa                 <  2:1591309/93‑1 (MQ=255)
                                                                  gACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCa                 <  1:1989457/93‑1 (MQ=255)
                                                                    cGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaat               <  2:38772/93‑1 (MQ=255)
                                                                    cGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaat               <  2:677799/93‑1 (MQ=255)
                                                                     gAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaata              <  2:3303625/93‑1 (MQ=255)
                                                                     gAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaata              <  1:1619901/93‑1 (MQ=255)
                                                                      aaTGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataa             <  1:1627348/93‑1 (MQ=255)
                                                                      aaTGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataa             <  2:3737841/93‑1 (MQ=255)
                                                                          aCGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataaaaat         >  1:2392051/1‑93 (MQ=255)
                                                                          aCGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataaaaat         >  1:4641761/1‑93 (MQ=255)
                                                                           cgcgGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataaaaata        <  2:2987267/93‑1 (MQ=255)
                                                                           cgcgGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataaaaata        <  2:2961499/93‑1 (MQ=255)
                                                                             cgGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataaaaataaa      <  1:4518208/93‑1 (MQ=255)
                                                                               ggCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGt    <  2:3103718/93‑1 (MQ=255)
                                                                               ggCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGt    <  2:4693929/93‑1 (MQ=255)
                                                                                gCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTa   <  2:3626138/93‑1 (MQ=255)
                                                                                gCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTa   <  2:2550317/93‑1 (MQ=255)
                                                                                gCGATAGCAAGATTATCCTGCAGCGTCAGGCGACGGATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTa   <  2:4108960/93‑1 (MQ=255)
                                                                                 cGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTaa  <  1:3288213/93‑1 (MQ=255)
                                                                                        |                                                                                       
CGTCAAGATTCGGATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCGAGGATTATCCTGCAGCGTCAGGCGACGATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAA  >  NC_000913/4476953‑4477126

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CAGCTCGTCAAGATTCGGATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCGAGGATTATCCTGCAGCGTCAGGCGACGATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTA  >  NC_000913/4476948‑4477129
                                                                                             |                                                                                         
CAGCTCGTCAAGATTCGGATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCA                                                                                           >  SRR4302544.1267613/1‑94 (MQ=60)
    TCGTCAAGATTCGGATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGAT                                                                                      <  SRR4302544.3554462/94‑1 (MQ=60)
     CGTCAAGATTCGGATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATT                                                                                     <  SRR4302544.3146948/94‑1 (MQ=60)
      GTCAAGATTCGGATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTA                                                                                    >  SRR4302544.4487678/1‑94 (MQ=60)
        TAAGATGCGGATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATC                                                                                  <  SRR4302544.490445/94‑1 (MQ=60)
        CAAGATTCGGATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCTGCGGCTACCTGACGAATTACGCGGGCGATAGCAGGATTATC                                                                                  <  SRR4302544.490446/94‑1 (MQ=60)
          AGATTCGGATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCT                                                                                <  SRR4302544.4694676/94‑1 (MQ=60)
          AGATTCGGATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCT                                                                                >  SRR4302544.1445752/1‑94 (MQ=60)
               CGGATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGC                                                                           <  SRR4302544.2703035/94‑1 (MQ=60)
                cggtcggccgCGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCG                                                                          <  SRR4302544.395537/84‑1 (MQ=60)
                GGATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGCCGACTGACGCTAGCGATAGCAGGATTATCCTGCAGCG                                                                          >  SRR4302544.1368328/1‑94 (MQ=60)
                GGATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCG                                                                          >  SRR4302544.1368327/1‑94 (MQ=60)
                 GATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGT                                                                         <  SRR4302544.534474/94‑1 (MQ=60)
                     GGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGG                                                                     <  SRR4302544.175758/94‑1 (MQ=60)
                     GGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGG                                                                     <  SRR4302544.3388699/94‑1 (MQ=60)
                      GCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATCCCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGC                                                                    <  SRR4302544.2576718/94‑1 (MQ=60)
                      GCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGC                                                                    <  SRR4302544.4753761/94‑1 (MQ=60)
                      GCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGC                                                                    >  SRR4302544.2091508/1‑94 (MQ=60)
                        GACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGA                                                                  <  SRR4302544.4312145/94‑1 (MQ=60)
                         ACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGAC                                                                 >  SRR4302544.2229610/1‑94 (MQ=60)
                             TGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGcat                                                             >  SRR4302544.4223595/1‑91 (MQ=60)
                               TAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTA                                                           >  SRR4302544.1431195/1‑94 (MQ=60)
                                    TTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGT                                                      >  SRR4302544.1545176/1‑94 (MQ=60)
                                    TTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGT                                                      <  SRR4302544.4328465/94‑1 (MQ=60)
                                    TTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGT                                                      <  SRR4302544.1928880/94‑1 (MQ=60)
                                       ATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGAT                                                   <  SRR4302544.1885818/94‑1 (MQ=60)
                                        TCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATT                                                  <  SRR4302544.4205334/94‑1 (MQ=60)
                                        TCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATT                                                  <  SRR4302544.3557029/94‑1 (MQ=60)
                                             GGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGC                                             >  SRR4302544.3186292/1‑94 (MQ=60)
                                                AAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCA                                          <  SRR4302544.3889297/94‑1 (MQ=60)
                                                 AGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAA                                         >  SRR4302544.1632493/1‑94 (MQ=60)
                                                  GACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAAT                                        >  SRR4302544.2642701/1‑94 (MQ=60)
                                                   ACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGGGGCGCAATG                                       >  SRR4302544.4645718/1‑94 (MQ=60)
                                                   ACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATG                                       <  SRR4302544.2938103/94‑1 (MQ=60)
                                                   ACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATG                                       <  SRR4302544.1315492/94‑1 (MQ=60)
                                                      GTATTCTGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTA                                    <  SRR4302544.1521271/94‑1 (MQ=60)
                                                        ATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATT                                  >  SRR4302544.2936380/1‑94 (MQ=60)
                                                          TCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCA                                <  SRR4302544.4545760/94‑1 (MQ=60)
                                                          TCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCA                                <  SRR4302544.1154798/94‑1 (MQ=60)
                                                          TCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGCCGCAATGTTATTCA                                >  SRR4302544.484001/1‑94 (MQ=60)
                                                           CGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCAT                               <  SRR4302544.1993058/94‑1 (MQ=60)
                                                           CGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCAT                               >  SRR4302544.2886181/1‑94 (MQ=60)
                                                             GCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAG                             <  SRR4302544.4211405/94‑1 (MQ=60)
                                                              CGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTTTTCATAGT                            >  SRR4302544.1416834/1‑94 (MQ=60)
                                                               GGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTT                           >  SRR4302544.3770563/1‑94 (MQ=60)
                                                               GGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTT                           <  SRR4302544.4042634/94‑1 (MQ=60)
                                                               GGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTT                           <  SRR4302544.2425657/94‑1 (MQ=60)
                                                                GATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTT                          <  SRR4302544.187313/94‑1 (MQ=60)
                                                                GATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTT                          <  SRR4302544.1463285/94‑1 (MQ=60)
                                                                  TACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAA                        >  SRR4302544.3434406/1‑94 (MQ=60)
                                                                  TACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAA                        >  SRR4302544.1250748/1‑94 (MQ=60)
                                                                    CCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACC                      <  SRR4302544.1176259/94‑1 (MQ=60)
                                                                      TGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCA                    <  SRR4302544.1056994/94‑1 (MQ=60)
                                                                      TGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCA                    <  SRR4302544.3882472/94‑1 (MQ=60)
                                                                       GACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAA                   <  SRR4302544.2008770/94‑1 (MQ=60)
                                                                          GAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAA                <  SRR4302544.1635860/94‑1 (MQ=60)
                                                                           AATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAA               <  SRR4302544.1643375/94‑1 (MQ=60)
                                                                              GACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAAT            >  SRR4302544.4687735/1‑94 (MQ=60)
                                                                              GACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAAT            >  SRR4302544.2414714/1‑94 (MQ=60)
                                                                                  CGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAG        <  SRR4302544.4563454/94‑1 (MQ=60)
                                                                                   GGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGT       >  SRR4302544.2298114/1‑94 (MQ=60)
                                                                                      CGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAAT    <  SRR4302544.3319542/94‑1 (MQ=60)
                                                                                        ATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTA  >  SRR4302544.1352899/1‑94 (MQ=60)
                                                                                             |                                                                                         
CAGCTCGTCAAGATTCGGATCGGCGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCGAGGATTATCCTGCAGCGTCAGGCGACGATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTA  >  NC_000913/4476948‑4477129

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