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,013,946 G→A Q26Q (CAG→CAA fliH → FliH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,013,9460GA97.6% 426.3 / ‑3.0 124Q26Q (CAG→CAAfliHFliH
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (60/61);  total (62/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 128 total aligned reads displayed.

GCGAGATGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAGGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACTGCAAAT  >  NC_000913/2013857‑2014035
                                                                                         |                                                                                         
  gagaTGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCaga                                                                                      <  1:1067130/93‑1 (MQ=255)
  gagaTGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCaga                                                                                      >  2:429726/1‑93 (MQ=255)
  gagaTGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCaga                                                                                      >  1:1957852/1‑93 (MQ=255)
   agaTGGTAATTGGCAGCGGCGAGGATACCTATGTCGGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCagag                                                                                     >  1:837363/1‑93 (MQ=255)
    gaTGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGt                                                                                    >  1:4266829/1‑93 (MQ=255)
     aTGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGtt                                                                                   <  2:1037846/93‑1 (MQ=255)
      tGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGttt                                                                                  <  1:2432689/93‑1 (MQ=255)
       ggTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTtg                                                                                 <  2:1390095/93‑1 (MQ=255)
       ggTAATTGGCAGCGGCGAAGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTtg                                                                                 >  1:2452452/1‑93 (MQ=255)
        gTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTtgt                                                                                <  2:225556/93‑1 (MQ=255)
        gTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTtgt                                                                                <  1:4518551/93‑1 (MQ=255)
         tAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTtgtg                                                                               >  1:2614370/1‑93 (MQ=255)
           aTTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGcc                                                                             <  1:2176072/93‑1 (MQ=255)
           aTTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGcc                                                                             <  1:4563310/93‑1 (MQ=255)
            ttGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGccc                                                                            >  2:3461712/1‑93 (MQ=255)
              ggCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCAt                                                                          <  1:1565251/93‑1 (MQ=255)
                 cgcggcgAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGt                                                                       <  1:1788607/92‑1 (MQ=255)
                 aGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGt                                                                       <  2:2764747/93‑1 (MQ=255)
                 aGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGt                                                                       <  2:3201097/93‑1 (MQ=255)
                  gcggcgAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTc                                                                      >  1:2256185/1‑93 (MQ=255)
                    ggcgAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCgg                                                                    >  2:3341964/1‑92 (MQ=255)
                    ggcgAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGa                                                                    <  1:3121124/93‑1 (MQ=255)
                     gcgAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAg                                                                   <  2:1170833/93‑1 (MQ=255)
                      cgAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGc                                                                  <  1:1299538/93‑1 (MQ=255)
                       gAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGcc                                                                 >  2:1758925/1‑93 (MQ=255)
                        aGGATGCCTATGTCTGATACTCTGCCGTGGAACACCTGGACGCCGGACGTTCTCGCGCCACCACAGGCAGAGTTTGTGCCCATAGTCGGGCCg                                                                <  2:3609250/93‑1 (MQ=255)
                        aGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCg                                                                >  2:2964959/1‑93 (MQ=255)
                        aGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCg                                                                >  1:1061865/1‑93 (MQ=255)
                        aGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCg                                                                >  1:2287969/1‑93 (MQ=255)
                          gATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGga                                                              >  2:1757651/1‑93 (MQ=255)
                             aCCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGgaaga                                                           >  2:201859/1‑93 (MQ=255)
                             aCCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGgaaga                                                           >  1:2383825/1‑93 (MQ=255)
                               cTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGaaa                                                         >  1:1088720/1‑93 (MQ=255)
                               cTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGaaa                                                         <  2:1673638/93‑1 (MQ=255)
                               cTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGaaa                                                         >  1:52857/1‑93 (MQ=255)
                                tATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAAc                                                        >  2:866804/1‑93 (MQ=255)
                                tATGTCTGATAATCTGCAGTGGAAAACCTGGACGCCGGACGATCTCGCGCCTCGACAGGCAGAGATTGAGCCCATAGTAGAGCCGGAAGCAAc                                                        >  2:2765103/1‑93 (MQ=255)
                                aaTGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAAc                                                        <  1:2138140/92‑1 (MQ=255)
                                  tGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACca                                                      >  1:2846575/1‑93 (MQ=255)
                                  tGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACca                                                      >  1:2242118/1‑93 (MQ=255)
                                  tGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATATCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACca                                                      >  2:1281635/1‑93 (MQ=255)
                                   gTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACcat                                                     >  1:488647/1‑93 (MQ=255)
                                   gTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACcat                                                     >  1:488648/1‑93 (MQ=255)
                                   gTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACcat                                                     <  2:2639153/93‑1 (MQ=255)
                                    tCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACcatc                                                    <  2:3402641/93‑1 (MQ=255)
                                    tCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACcatc                                                    >  2:4183240/1‑93 (MQ=255)
                                    tCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACcatc                                                    <  2:568748/93‑1 (MQ=255)
                                      tGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAGGCAGAGTTTGTGCCCATAGACGAGCCGGAAGAAACcatcat                                                  >  1:3995917/1‑93 (MQ=255)
                                      tGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACcatcat                                                  >  2:1822329/1‑93 (MQ=255)
                                       gATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCTCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCAtt                                                 >  1:3623875/1‑93 (MQ=255)
                                       gATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCAtt                                                 >  1:259800/1‑93 (MQ=255)
                                          aaTCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGaa                                              >  2:3131698/1‑93 (MQ=255)
                                          aaTCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGaa                                              >  2:2512039/1‑93 (MQ=255)
                                            tCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAaga                                            >  2:4620046/1‑93 (MQ=255)
                                            tCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAaga                                            <  1:3609702/93‑1 (MQ=255)
                                            tCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAaga                                            <  1:2307691/93‑1 (MQ=255)
                                             cTGCCGTGGAAAACCTGGACGCCGGACGAACTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAagag                                           <  2:1703111/93‑1 (MQ=255)
                                              tGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAgg                                          >  1:947834/1‑93 (MQ=255)
                                              tGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAgg                                          >  1:2835583/1‑93 (MQ=255)
                                              cgCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAgg                                          <  2:3534172/92‑1 (MQ=255)
                                               gCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGc                                         <  1:866606/93‑1 (MQ=255)
                                                  gTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGa                                      >  1:1769486/1‑93 (MQ=255)
                                                  gTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGa                                      <  1:1644149/93‑1 (MQ=255)
                                                   tGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGaa                                     >  1:3775076/1‑93 (MQ=255)
                                                     gAAAACCTGGACGCCGGACGATCTCGCGCCACCACAATCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAAcc                                   >  2:3479981/1‑93 (MQ=255)
                                                     gAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAAcc                                   >  2:1389823/1‑93 (MQ=255)
                                                     gAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAAcc                                   >  1:3464317/1‑93 (MQ=255)
                                                       aaaCCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCa                                 <  1:3779049/93‑1 (MQ=255)
                                                       aaaCCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCa                                 >  2:676911/1‑93 (MQ=255)
                                                       aaaCCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCa                                 <  2:838795/93‑1 (MQ=255)
                                                        aaCCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAg                                <  2:3121880/93‑1 (MQ=255)
                                                         aCCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGc                               <  1:3386147/93‑1 (MQ=255)
                                                           cTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCt                             <  2:2154870/93‑1 (MQ=255)
                                                             ggACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTg                           <  2:54090/93‑1 (MQ=255)
                                                             ggACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTg                           >  1:3772748/1‑93 (MQ=255)
                                                               aCGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGag                         <  2:4580393/93‑1 (MQ=255)
                                                                 gcgggACGACCTCGCGCCACCACAAGCAGAGTTTGTGCCCAAAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGagca                       <  2:1246444/90‑1 (MQ=255)
                                                                 gCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGagca                       <  1:464091/93‑1 (MQ=255)
                                                                 gCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGagca                       >  1:1565359/1‑93 (MQ=255)
                                                                    ggACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAgcagca                    >  1:4320476/1‑93 (MQ=255)
                                                                      aCGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAAc                  <  2:412157/93‑1 (MQ=255)
                                                                       cGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACt                 <  2:3523393/93‑1 (MQ=255)
                                                                        gATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTg                <  2:2112828/93‑1 (MQ=255)
                                                                        gATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTg                <  2:2085047/93‑1 (MQ=255)
                                                                        gATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTg                <  2:263176/93‑1 (MQ=255)
                                                                          tctcGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGgc              >  1:40504/1‑93 (MQ=255)
                                                                              gcgcCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCaa          <  2:2233890/93‑1 (MQ=255)
                                                                               cgcCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAAc         <  1:1625757/93‑1 (MQ=255)
                                                                               cgcCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAAc         <  2:2059595/93‑1 (MQ=255)
                                                                                gcCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACt        <  2:834889/93‑1 (MQ=255)
                                                                                 ccaccaCAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGCACCCAGCCTTGAACAGCAACTGGCGCAACtc       >  2:1201987/1‑92 (MQ=255)
                                                                                 ccaccaCAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACTg       >  2:2663913/1‑93 (MQ=255)
                                                                                 ccaccaCAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACTg       <  2:4720715/93‑1 (MQ=255)
                                                                                  caccacAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACTGc      >  2:2719256/1‑93 (MQ=255)
                                                                                   accacAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACTGCa     >  1:4399248/1‑93 (MQ=255)
                                                                                   accacAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACTGCa     >  2:1024501/1‑93 (MQ=255)
                                                                                   accacAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACTGCa     >  1:875939/1‑93 (MQ=255)
                                                                                     cacaGGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACTGCaaa   >  2:3995363/1‑93 (MQ=255)
                                                                                     cacaAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACTGCaaa   <  1:3362656/93‑1 (MQ=255)
                                                                                      acaAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACTGCAAAt  <  2:4050631/93‑1 (MQ=255)
                                                                                         |                                                                                         
GCGAGATGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAGGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACTGCAAAT  >  NC_000913/2013857‑2014035

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CTGGCGAGATGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAGGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACTGCAAATGCA  >  NC_000913/2013854‑2014038
                                                                                            |                                                                                            
CTGGCGAGATGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAG                                                                                             >  SRR4302544.3557054/1‑94 (MQ=60)
    CGAGATGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGA                                                                                         >  SRR4302544.1976844/1‑94 (MQ=60)
     GAGATGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAG                                                                                        <  SRR4302544.1077742/94‑1 (MQ=60)
     GAGATGGTAATTGGCAGCGGCGAGGATACCTATGTCGGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAG                                                                                        >  SRR4302544.845616/1‑94 (MQ=60)
      AGATGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGT                                                                                       >  SRR4302544.4310378/1‑94 (MQ=60)
         TGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTG                                                                                    <  SRR4302544.2455867/94‑1 (MQ=60)
         TGGTAATTGGCAGCGGCGAAGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTG                                                                                    >  SRR4302544.2475846/1‑94 (MQ=60)
           GTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTG                                                                                  <  SRR4302544.4563802/94‑1 (MQ=60)
           GTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTG                                                                                  >  SRR4302544.2639470/1‑94 (MQ=60)
              ATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCC                                                                               <  SRR4302544.2196688/94‑1 (MQ=60)
              ATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCC                                                                               <  SRR4302544.4608808/94‑1 (MQ=60)
                 GGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATA                                                                            <  SRR4302544.1580772/94‑1 (MQ=60)
                   CAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGT                                                                          >  SRR4302544.527201/1‑94 (MQ=60)
                   CAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGT                                                                          <  SRR4302544.1723469/94‑1 (MQ=60)
                    CGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTC                                                                         <  SRR4302544.1806090/94‑1 (MQ=60)
                    AGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCCCCACAAGCAGAGTTTGTGCCCATAGTC                                                                         >  SRR4302544.3647371/1‑94 (MQ=60)
                    AGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTC                                                                         >  SRR4302544.2277415/1‑94 (MQ=60)
                       GGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAG                                                                      <  SRR4302544.3150792/94‑1 (MQ=60)
                         CGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCC                                                                    <  SRR4302544.1312322/94‑1 (MQ=60)
                          GAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCG                                                                   >  SRR4302544.1072421/1‑94 (MQ=60)
                          GAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCG                                                                   >  SRR4302544.2309541/1‑94 (MQ=60)
                               TACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGA                                                              >  SRR4302544.2406398/1‑94 (MQ=60)
                                 CCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGAGGATATCGCGCCACCACGAGAAGCGTTTGTGCGCATAGTCGAGCgggagacaa                                                            >  SRR4302544.53396/1‑85 (MQ=60)
                                 CCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAA                                                            >  SRR4302544.1099546/1‑94 (MQ=60)
                                 CCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAA                                                            >  SRR4302544.53395/1‑94 (MQ=60)
                                   AATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACC                                                          <  SRR4302544.2158566/94‑1 (MQ=60)
                                    ATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCA                                                         >  SRR4302544.2263210/1‑94 (MQ=60)
                                    ATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCA                                                         >  SRR4302544.2874067/1‑94 (MQ=60)
                                     TGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCAT                                                        >  SRR4302544.493438/1‑94 (MQ=60)
                                     TGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCAT                                                        >  SRR4302544.493439/1‑94 (MQ=60)
                                        CTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAGGCAGAGTTTGTGCCCATAGACGAGCCGGAAGAAACCATCAT                                                     >  SRR4302544.4035691/1‑94 (MQ=60)
                                        CTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCAT                                                     <  SRR4302544.4675707/94‑1 (MQ=60)
                                         TGATAATCTGCCGTGGAAACCCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATT                                                    <  SRR4302544.2282104/94‑1 (MQ=60)
                                         TGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCTCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATT                                                    >  SRR4302544.3659984/1‑94 (MQ=60)
                                         TGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATT                                                    >  SRR4302544.262551/1‑94 (MQ=60)
                                             AATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAG                                                <  SRR4302544.3270633/94‑1 (MQ=60)
                                               TCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAG                                              <  SRR4302544.3644922/94‑1 (MQ=60)
                                               TCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAG                                              <  SRR4302544.2329504/94‑1 (MQ=60)
                                                CTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGG                                             >  SRR4302544.957161/1‑94 (MQ=60)
                                                CTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGG                                             >  SRR4302544.2862975/1‑94 (MQ=60)
                                                  GCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCT                                           <  SRR4302544.875181/94‑1 (MQ=60)
                                                  GCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCT                                           <  SRR4302544.2067545/94‑1 (MQ=60)
                                                   CCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTG                                          <  SRR4302544.1579973/94‑1 (MQ=60)
                                                    CGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGA                                         <  SRR4302544.4020935/94‑1 (MQ=60)
                                                    CGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGA                                         >  SRR4302544.1786804/1‑94 (MQ=60)
                                                     GTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAA                                        >  SRR4302544.3812490/1‑94 (MQ=60)
                                                     GTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAA                                        <  SRR4302544.1660359/94‑1 (MQ=60)
                                                       GGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACC                                      >  SRR4302544.3497365/1‑94 (MQ=60)
                                                        GAAAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCC                                     <  SRR4302544.242801/94‑1 (MQ=60)
                                                          AAACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAG                                   <  SRR4302544.3816483/94‑1 (MQ=60)
                                                            ACCTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCC                                 <  SRR4302544.3418502/94‑1 (MQ=60)
                                                              CTGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTT                               <  SRR4302544.2782106/94‑1 (MQ=60)
                                                               TGGACGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTG                              >  SRR4302544.3810141/1‑94 (MQ=60)
                                                                   CGCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCA                          >  SRR4302544.1580880/1‑94 (MQ=60)
                                                                    GCCGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAG                         <  SRR4302544.468646/94‑1 (MQ=60)
                                                                      CGGACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCA                       >  SRR4302544.4364484/1‑94 (MQ=60)
                                                                        GACGATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAAC                     <  SRR4302544.900995/94‑1 (MQ=60)
                                                                           GATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGG                  >  SRR4302544.4348149/1‑94 (MQ=60)
                                                                            ATCTCGCGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGC                 >  SRR4302544.40907/1‑94 (MQ=60)
                                                                                  CGCCACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACT           <  SRR4302544.1641776/94‑1 (MQ=60)
                                                                                     CACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACTGCA        >  SRR4302544.884628/1‑94 (MQ=60)
                                                                                     CACCACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACTGCA        >  SRR4302544.4443713/1‑94 (MQ=60)
                                                                                        CACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACTGCAAAT     <  SRR4302544.3394747/94‑1 (MQ=60)
                                                                                         ACAAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACTGCAAATG    >  SRR4302544.847191/1‑94 (MQ=60)
                                                                                           AAGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACTGCAAATGCA  <  SRR4302544.3005517/94‑1 (MQ=60)
                                                                                            |                                                                                            
CTGGCGAGATGGTAATTGGCAGCGGCGAGGATACCTATGTCTGATAATCTGCCGTGGAAAACCTGGACGCCGGACGATCTCGCGCCACCACAGGCAGAGTTTGTGCCCATAGTCGAGCCGGAAGAAACCATCATTGAAGAGGCTGAACCCAGCCTTGAGCAGCAACTGGCGCAACTGCAAATGCA  >  NC_000913/2013854‑2014038

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