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 647 270.1 21258240 98.4% 20918108 94.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,008,391 +C coding (474/516 nt) b4555 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,008,3911.C99.5% 718.1 / ‑4.1 211coding (474/516 nt)b4555hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (96/114);  total (97/114)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 252 total aligned reads displayed.

GCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGC‑‑‑‑G‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGT  >  NC_000913/4008299‑4008459
                                                                                                 |                                                                                 
gCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAg                                                                                    <  2:2159849/94‑1 (MQ=255)
gCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAg                                                                                    <  1:3717193/94‑1 (MQ=255)
gCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAAGACCCTGCGCTGGAAg                                                                                    >  1:1159954/1‑94 (MQ=255)
ccTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAg                                                                                    <  1:1593911/93‑1 (MQ=255)
 cTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgc                                                                                   <  2:10178825/94‑2 (MQ=255)
 cTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgc                                                                                   <  2:1641849/94‑2 (MQ=255)
 cTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgc                                                                                   <  1:7514854/94‑2 (MQ=255)
 cTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAAGCCATGCGCTGGAAgc                                                                                   >  2:8958609/1‑93 (MQ=255)
  tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcg                                                                                  >  2:7905413/1‑92 (MQ=255)
  tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcg                                                                                  <  1:7332984/94‑3 (MQ=255)
  tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAACGCCATGCGCTGGAAgcg                                                                                  >  2:2825397/1‑92 (MQ=255)
   cGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcgg                                                                                 >  1:9949682/1‑91 (MQ=255)
   cGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAAGGCCATGCGCTGGAAgcgg                                                                                 >  1:882793/1‑91 (MQ=255)
    gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGGGGGGGTTCAAAAGGCCATGCGCTGGAAgcggt                                                                                >  2:127100/1‑90 (MQ=255)
    gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcggt                                                                                <  2:895445/94‑5 (MQ=255)
    gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcggt                                                                                <  1:7037270/94‑5 (MQ=255)
         ggCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGGAGCGGTGGC‑‑‑‑G‑‑‑‑a                                                                   <  1:10385574/94‑1 (MQ=255)
         ggCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑G‑‑‑‑a                                                                   >  2:5640730/1‑94 (MQ=255)
          gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGGGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑G‑‑‑‑aa                                                                  >  2:8144401/1‑94 (MQ=255)
          gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑G‑‑‑‑aa                                                                  <  1:2817236/94‑1 (MQ=255)
          gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑G‑‑‑‑aa                                                                  <  2:3098867/94‑1 (MQ=255)
          gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑G‑‑‑‑aa                                                                  <  1:576873/94‑1 (MQ=255)
           cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑G‑‑‑‑aag                                                                 <  2:4432431/94‑2 (MQ=255)
             aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑G‑‑‑‑aagaa                                                               <  2:7118894/94‑4 (MQ=255)
             aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAAGGCATGCGCTGGAAGCGGTGGC‑‑‑‑G‑‑‑‑aagaa                                                               >  2:3979568/1‑91 (MQ=255)
              tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑G‑‑‑‑aagaaa                                                              >  1:834521/1‑90 (MQ=255)
              tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑G‑‑‑‑aagaaa                                                              >  2:5465632/1‑90 (MQ=255)
                                              gTCTGGAAGTGATGGCGGGCGGCGTTTGAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                  >  2:5460470/1‑94 (MQ=255)
                                              gTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                  <  2:6030018/94‑1 (MQ=255)
                                               tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGttt                                 >  1:790202/1‑93 (MQ=255)
                                               tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGttt                                 <  1:5051237/94‑2 (MQ=255)
                                                cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttg                                <  2:5695884/94‑3 (MQ=255)
                                                cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttg                                >  2:5958420/1‑92 (MQ=255)
                                                cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttg                                >  2:6194405/1‑92 (MQ=255)
                                                 tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg                               >  1:4852832/1‑91 (MQ=255)
                                                 tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg                               >  1:4485585/1‑91 (MQ=255)
                                                 tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg                               >  1:8323723/1‑91 (MQ=255)
                                                 tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg                               <  1:9467560/94‑4 (MQ=255)
                                                 tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg                               <  2:6631985/94‑4 (MQ=255)
                                                  ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                              <  2:8511693/94‑5 (MQ=255)
                                                        gATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                        <  1:1219882/94‑1 (MQ=255)
                                                        gATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGACTGAAGGATTGTATTGCGTTTGGTGACggg                        >  2:4202765/1‑94 (MQ=255)
                                                         aTGGCGGGCGGCGTTTCAAAAAGCCATGCGCTGGGAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACgggg                       >  2:4647869/1‑93 (MQ=255)
                                                          tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                      <  1:10271385/94‑1 (MQ=255)
                                                          tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                      >  2:10109989/1‑94 (MQ=255)
                                                          tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                      <  1:4573594/94‑1 (MQ=255)
                                                          tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                      >  1:7783819/1‑94 (MQ=255)
                                                           ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATg                     <  2:1352777/94‑1 (MQ=255)
                                                           ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATg                     >  1:1980290/1‑94 (MQ=255)
                                                             cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                   <  2:9856631/94‑1 (MQ=255)
                                                             cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                   >  2:8736980/1‑94 (MQ=255)
                                                             cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                   >  2:4511526/1‑94 (MQ=255)
                                                             cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                   >  2:2262793/1‑94 (MQ=255)
                                                             cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                   <  2:1885333/94‑1 (MQ=255)
                                                              gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaag                  >  2:3717254/1‑93 (MQ=255)
                                                              gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                  <  2:6723938/94‑1 (MQ=255)
                                                               ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacg                 >  1:10299843/1‑94 (MQ=255)
                                                               ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacg                 <  2:3462535/94‑1 (MQ=255)
                                                               ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacg                 <  2:799367/94‑1 (MQ=255)
                                                                gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga                <  1:5201171/94‑1 (MQ=255)
                                                                gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga                <  1:9298965/94‑1 (MQ=255)
                                                                gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga                >  2:6573983/1‑94 (MQ=255)
                                                                 gggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac               <  2:2430425/93‑1 (MQ=255)
                                                                 cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GCAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac               >  2:492205/1‑94 (MQ=255)
                                                                 cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac               <  1:614204/94‑1 (MQ=255)
                                                                 cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac               <  1:1914560/94‑1 (MQ=255)
                                                                 cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac               <  2:4749689/94‑1 (MQ=255)
                                                                  ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg              <  2:745645/94‑1 (MQ=255)
                                                                  ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg              >  1:3482095/1‑94 (MQ=255)
                                                                  ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg              >  1:4117323/1‑94 (MQ=255)
                                                                  ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg              >  1:4820196/1‑94 (MQ=255)
                                                                  ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacc              >  2:10361702/1‑93 (MQ=255)
                                                                   gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc             <  2:3943776/94‑1 (MQ=255)
                                                                    cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGGTGAACGACGcc            >  1:8688707/1‑94 (MQ=255)
                                                                    cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc            <  2:4805074/94‑1 (MQ=255)
                                                                    cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc            <  2:3133526/94‑1 (MQ=255)
                                                                    cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc            >  2:9058941/1‑94 (MQ=255)
                                                                    cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc            >  2:9201656/1‑94 (MQ=255)
                                                                     gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCg           >  1:7662318/1‑94 (MQ=255)
                                                                     gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCg           >  2:4874566/1‑94 (MQ=255)
                                                                     gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCg           >  2:3693511/1‑94 (MQ=255)
                                                                        ttcaaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGCACGACGCCCaa         >  1:2731701/5‑94 (MQ=255)
                                                                        ttcaaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaa         >  2:7233226/5‑94 (MQ=255)
                                                                         tcaaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa        >  2:4588664/4‑94 (MQ=255)
                                                                         tcaaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa        >  1:4648492/4‑94 (MQ=255)
                                                                          caaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt       <  1:8151447/92‑1 (MQ=255)
                                                                          caaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt       >  2:1724717/3‑94 (MQ=255)
                                                                          caaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt       <  1:2514554/92‑1 (MQ=255)
                                                                          caaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt       <  1:2292684/92‑1 (MQ=255)
                                                                          caaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt       <  2:415905/92‑1 (MQ=255)
                                                                          caaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt       <  2:5117278/92‑1 (MQ=255)
                                                                           aaaaGGCCATGCGCTGGAAGGGGTGGCGAAGG‑‑‑‑AACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAtc      >  2:865845/2‑93 (MQ=255)
                                                                           aaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg      <  2:705039/93‑1 (MQ=255)
                                                                           aaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg      >  1:5472649/2‑94 (MQ=255)
                                                                            aaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc     <  1:1012806/94‑1 (MQ=255)
                                                                              aGGCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCCACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTg   <  2:2778903/94‑1 (MQ=255)
                                                                                ggCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGt  <  2:3587287/94‑1 (MQ=255)
                                                                                ggCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGt  >  1:3013830/1‑94 (MQ=255)
                                                                                ggCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGt  <  1:6686780/94‑1 (MQ=255)
                                                                                ggCCATGCGCTGGAAGCGGTGGC‑‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGt  >  2:2802623/1‑94 (MQ=255)
                                                                                                 |                                                                                 
GCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGC‑‑‑‑G‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGT  >  NC_000913/4008299‑4008459

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 25 ≤ 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
GCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCG  >  NC_000913/4008309‑4008461
                                                                                       |                                                                             
GCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑AA                                                                    <  SRR4302545.2885434/94‑1 (MQ=60)
GCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑AA                                                                    <  SRR4302545.5624836/94‑1 (MQ=60)
GCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑AA                                                                    <  SRR4302545.589885/94‑1 (MQ=60)
                                                              tcaAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAA          >  SRR4302545.4767244/4‑94 (MQ=16)
                                                               caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT         <  SRR4302545.8579375/92‑1 (MQ=18)
                                                               caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT         <  SRR4302545.8480524/92‑1 (MQ=18)
                                                               caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT         <  SRR4302545.8435863/92‑1 (MQ=18)
                                                               caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT         <  SRR4302545.2348517/92‑1 (MQ=18)
                                                               caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT         >  SRR4302545.2354825/3‑94 (MQ=18)
                                                               caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT         <  SRR4302545.2575333/92‑1 (MQ=18)
                                                               caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT         <  SRR4302545.8387026/92‑1 (MQ=18)
                                                               AAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATG        >  SRR4302545.4910112/1‑94 (MQ=19)
                                                               AAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATG        >  SRR4302545.874445/1‑94 (MQ=19)
                                                               AAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATG        >  SRR4302545.5616317/1‑94 (MQ=19)
                                                               AAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATG        >  SRR4302545.3592584/1‑94 (MQ=19)
                                                               AAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATG        >  SRR4302545.312860/1‑94 (MQ=19)
                                                                 AAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGC       <  SRR4302545.1035173/94‑1 (MQ=21)
                                                                   AGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTG     <  SRR4302545.5627053/94‑1 (MQ=23)
                                                                   AGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTG     >  SRR4302545.287453/1‑94 (MQ=23)
                                                                   AGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTG     <  SRR4302545.10720851/94‑1 (MQ=23)
                                                                     GGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGT    >  SRR4302545.3087452/1‑94 (MQ=25)
                                                                     GGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGT    <  SRR4302545.6870088/94‑1 (MQ=25)
                                                                      GCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTC   <  SRR4302545.5377252/94‑1 (MQ=26)
                                                                      GCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTC   <  SRR4302545.7325608/94‑1 (MQ=26)
                                                                       CCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCG  <  SRR4302545.5176382/94‑1 (MQ=28)
                                                                       CCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCG  <  SRR4302545.1999026/94‑1 (MQ=28)
                                                                       CCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCG  >  SRR4302545.1430924/1‑94 (MQ=28)
                                                                       CCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCG  <  SRR4302545.10756070/94‑1 (MQ=28)
                                                                                       |                                                                             
GCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCG  >  NC_000913/4008309‑4008461

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