Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F1 I2 R1 1624 250.6 12620048 98.0% 12367647 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,008,393 C→G A159G (GCT→GGT)  b4555 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,008,3930CG100.0% 263.9 / NA 83A159G (GCT→GGT) b4555hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (34/49);  total (34/49)

GCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AA‑‑‑‑‑CTGGCTACAGCTGAAGGATTGATTGCG‑‑TT‑‑GTGACGGGATGAACGACGCCGAAATGCTG  >  NC_000913/4008309‑4008458
                                                                                       |                                                                                 
gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑a                                                                         >  1:882407/1‑93 (MQ=255)
gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑a                                                                         <  2:722109/93‑1 (MQ=255)
gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑a                                                                         >  2:2003960/1‑93 (MQ=255)
gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAAGGGCATGCGCTGGAAGCGGTGGCG‑‑‑‑a                                                                         >  2:4689360/1‑93 (MQ=255)
 cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aa                                                                        <  2:4357872/93‑1 (MQ=255)
 cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aa                                                                        <  2:2103880/93‑1 (MQ=255)
   aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                                      <  2:219070/93‑3 (MQ=255)
   aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                                      <  2:5143571/93‑3 (MQ=255)
   aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                                      >  1:4242722/1‑91 (MQ=255)
   aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                                      <  2:3925526/93‑3 (MQ=255)
   aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                                      >  1:3268770/1‑91 (MQ=255)
   aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                                      <  2:4328059/93‑3 (MQ=255)
   aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCACTGGCTGGAAGGGGTGGCG‑‑‑‑aaaa                                                                      >  1:1410130/1‑91 (MQ=255)
    tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaa                                                                     >  1:2474435/1‑90 (MQ=255)
    tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAAGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaa                                                                     >  2:2923090/1‑90 (MQ=255)
     cgcgTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaaa                                                                    <  1:3617388/93‑5 (MQ=255)
     cgcgTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaaa                                                                    <  2:2759649/93‑5 (MQ=255)
                                     tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tt                                 <  1:3067422/93‑1 (MQ=255)
                                     tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tt                                 >  2:5367497/1‑93 (MQ=255)
                                      cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑ttt                                >  2:2690426/1‑92 (MQ=255)
                                      cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑ttt                                <  1:3423220/93‑2 (MQ=255)
                                      cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑ttt                                >  2:2607369/1‑92 (MQ=255)
                                       tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttg                               >  2:1156237/1‑91 (MQ=255)
                                       tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttg                               <  1:1051706/93‑3 (MQ=255)
                                       tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttg                               >  1:5951006/1‑91 (MQ=255)
                                       tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttg                               >  1:5995372/1‑91 (MQ=255)
                                       tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttg                               >  2:4708075/1‑91 (MQ=255)
                                        ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttgg                              >  2:630739/1‑90 (MQ=255)
                                                tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGa                      <  1:100493/93‑1 (MQ=255)
                                                tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGa                      <  2:349093/93‑1 (MQ=255)
                                                tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGa                      <  2:5970655/93‑1 (MQ=255)
                                                tGGCGGGCGGCGTTTCAAAAAGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGa                      >  2:6045798/1‑93 (MQ=255)
                                                 ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGAt                     <  2:2148939/93‑1 (MQ=255)
                                                 ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑AAGAAACCGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGAt                     <  2:1686766/93‑1 (MQ=255)
                                                  gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATg                    <  2:2236559/93‑1 (MQ=255)
                                                  gcgggcggCGTGTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATg                    <  2:790947/93‑1 (MQ=255)
                                                   cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGa                   <  1:5478423/93‑1 (MQ=255)
                                                    gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGaa                  <  2:3445771/93‑1 (MQ=255)
                                                    gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGaa                  <  1:2635131/93‑1 (MQ=255)
                                                     ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAac                 <  2:2772791/93‑1 (MQ=255)
                                                     ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAac                 >  1:3993000/1‑93 (MQ=255)
                                                     ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAac                 <  1:4957056/93‑1 (MQ=255)
                                                     ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAac                 <  2:331868/93‑1 (MQ=255)
                                                      gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacg                >  1:4814006/1‑93 (MQ=255)
                                                      gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacg                <  1:2957114/93‑1 (MQ=255)
                                                       cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacga               <  2:2166732/93‑1 (MQ=255)
                                                       cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacga               <  1:2546223/93‑1 (MQ=255)
                                                        ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTT‑‑GTGACGGGATGAacgac              >  2:5773184/1‑93 (MQ=255)
                                                        ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgac              >  2:3119511/1‑93 (MQ=255)
                                                        ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgac              >  2:5281397/1‑93 (MQ=255)
                                                        ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgac              <  2:4498313/93‑1 (MQ=255)
                                                        ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgac              <  2:4279944/93‑1 (MQ=255)
                                                         gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgacg             >  2:4300379/1‑93 (MQ=255)
                                                          cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGc            >  2:5938791/1‑93 (MQ=255)
                                                           gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGcc           <  2:4808827/93‑1 (MQ=255)
                                                           gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGcc           <  1:4418546/93‑1 (MQ=255)
                                                             ttcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGa         <  2:2580522/89‑1 (MQ=255)
                                                             ttcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGa         >  2:5932837/5‑93 (MQ=255)
                                                             ttcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGa         >  1:6270373/5‑93 (MQ=255)
                                                              tcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGaa        >  2:476989/4‑93 (MQ=255)
                                                               caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGaaa       <  1:3434171/91‑1 (MQ=255)
                                                               caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGaaa       <  1:4553544/91‑1 (MQ=255)
                                                               caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGaaa       >  1:4825533/3‑93 (MQ=255)
                                                               caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGaaa       >  1:4825534/3‑93 (MQ=255)
                                                               caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGaaa       <  1:5766444/91‑1 (MQ=255)
                                                               caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGaaa       <  2:723500/91‑1 (MQ=255)
                                                               caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGaaa       <  1:995345/91‑1 (MQ=255)
                                                               caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGaaa       <  2:808001/91‑1 (MQ=255)
                                                                aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAAt      >  2:3678743/2‑93 (MQ=255)
                                                                aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAAt      >  2:4130459/2‑93 (MQ=255)
                                                                aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAAt      >  1:5825289/2‑93 (MQ=255)
                                                                aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAAt      >  2:1526869/2‑93 (MQ=255)
                                                                aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAAt      <  2:2946514/92‑1 (MQ=255)
                                                                 aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGCAATg     >  2:2085646/1‑93 (MQ=255)
                                                                 aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATg     >  1:290367/1‑93 (MQ=255)
                                                                  aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGc    <  1:6192030/93‑1 (MQ=255)
                                                                  aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGc    >  2:3848605/1‑93 (MQ=255)
                                                                  aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGc    <  1:2935716/93‑1 (MQ=255)
                                                                   aGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCt   >  2:6045723/1‑93 (MQ=255)
                                                                   aGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCt   <  2:193059/93‑1 (MQ=255)
                                                                   aGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCt   <  1:1197056/93‑1 (MQ=255)
                                                                   aGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCt   <  1:3735094/93‑1 (MQ=255)
                                                                    ggCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTg  <  2:2269624/93‑1 (MQ=255)
                                                                    ggCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTg  <  1:5385678/93‑1 (MQ=255)
                                                                    ggCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTg  <  1:241840/93‑1 (MQ=255)
                                                                    ggCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTg  <  1:4092752/93‑1 (MQ=255)
                                                                                       |                                                                                 
GCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AA‑‑‑‑‑CTGGCTACAGCTGAAGGATTGATTGCG‑‑TT‑‑GTGACGGGATGAACGACGCCGAAATGCTG  >  NC_000913/4008309‑4008458

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTC  >  NC_000913/4008308‑4008460
                                                                                       |                                                                            
GGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑A                                                                    >  SRR4302533.887590/1‑94 (MQ=60)
GGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑A                                                                    <  SRR4302533.5337889/94‑1 (MQ=60)
                                                               tcaAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAA         >  SRR4302533.4854775/4‑94 (MQ=16)
                                                               tcaAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAA         >  SRR4302533.4854776/4‑94 (MQ=16)
                                                                caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT        >  SRR4302533.5860484/3‑94 (MQ=18)
                                                                caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT        <  SRR4302533.3455385/92‑1 (MQ=18)
                                                                caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT        <  SRR4302533.5801257/92‑1 (MQ=18)
                                                                caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT        <  SRR4302533.1001159/92‑1 (MQ=18)
                                                                caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT        <  SRR4302533.4581160/92‑1 (MQ=18)
                                                                AAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATG       >  SRR4302533.292135/1‑94 (MQ=19)
                                                                   AAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCT     <  SRR4302533.6229189/94‑1 (MQ=21)
                                                                   AAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCT     <  SRR4302533.2953665/94‑1 (MQ=21)
                                                                    AGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTG    <  SRR4302533.1204246/94‑1 (MQ=23)
                                                                    AGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTG    >  SRR4302533.2757655/1‑94 (MQ=23)
                                                                    AGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTG    <  SRR4302533.3758008/94‑1 (MQ=23)
                                                                    AGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTG    >  SRR4302533.49854/1‑94 (MQ=23)
                                                                     GGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGT   <  SRR4302533.5418036/94‑1 (MQ=25)
                                                                     GGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGT   <  SRR4302533.4117665/94‑1 (MQ=25)
                                                                     GGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGT   <  SRR4302533.243338/94‑1 (MQ=25)
                                                                      GCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTC  >  SRR4302533.39867/1‑94 (MQ=26)
                                                                      GCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTC  <  SRR4302533.5623384/94‑1 (MQ=26)
                                                                      GCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTC  <  SRR4302533.3654174/94‑1 (MQ=26)
                                                                      GCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTC  <  SRR4302533.20499/94‑1 (MQ=26)
                                                                      GCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTC  <  SRR4302533.1174613/94‑1 (MQ=26)
                                                                                       |                                                                            
GGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTC  >  NC_000913/4008308‑4008460

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