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
A12 F6 I1 R1 78 1766.9 106897164 98.9% 105721295 100.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451_phaM1 5,312,212 2 bp→CG noncoding (469‑470/1522 nt) PP_16SG rRNA (16S ribosomal RNA) ←

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM15,312,2120TC83.7% 16.7 / 23.4 86noncoding (470/1522 nt)PP_16SG rRNA (16S ribosomal RNA)
Reads supporting (aligned to +/- strand):  ref base T (5/9);  new base C (72/0);  total (77/9)
Fisher's exact test for biased strand distribution p-value = 4.35e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*AE015451_phaM15,312,2130CG84.7% 17.0 / 33.3 85noncoding (469/1522 nt)PP_16SG rRNA (16S ribosomal RNA)
Reads supporting (aligned to +/- strand):  ref base C (4/9);  new base G (72/0);  total (76/9)
Fisher's exact test for biased strand distribution p-value = 1.74e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCACCTACGCGCGCTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCTGTATTACCGCGGCTGCTGGCACAGAGTTAGCCGGTGCTTATTCTGTCGGTAACGTCAAAACAGCAAGGTATTAACTTACTGCCCTTCCTCCCAACTTAAAGTGCTTTACAATCCGAAGACCTTCTTCACACACGCGGCATGGC  >  AE015451_phaM1/5312114‑5312309
                                                                                                  ||                                                                                                
ccccccTACGCGCGCTTTACGCCCAGTAATTCCGATTACCGCTTGCACCCTCTGTCTTACCGCGGCTGCTGGCACAGAGTTAGCCGGTGCTTATTCTGTc                                                                                                  >  7:29243142/5‑100 (MQ=1)
   acttACGCGCGCTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCTGTATTACCGCGGCTGCTGGCACAGAGTTAGCCGGTGCTTATTCTGTCGGt                                                                                               >  8:42862188/4‑100 (MQ=1)
    tctgcgcgcgCTTTACGCCCAGTAATTCCGATTAGCGCTTGCACCCTCTGTATTACCGCGGCTGCTGGCACAGAGTTAGCCGGTGCTTATTCTGTCGGTa                                                                                              >  8:48094176/5‑100 (MQ=0)
     ctgcgcgcgCTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCTGTATTACCGCGGCTGCGGGCACAGAGTTAGCCGGTGCTTATTCTGTCGGTaa                                                                                             <  7:36344425/97‑1 (MQ=0)
      tgcgtgcgcTTTGCGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:353687/6‑100 (MQ=1)
      tgcgtgcgcTTTACGNCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:284170/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCCCCCTCCGTCTTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:103786/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCTTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:289024/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:340535/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:25645/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:33722/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:326833/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:308910/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:308474/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:28466/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:282537/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:276453/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:276452/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:101881/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:263824/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:260951/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:258524/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:276074/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:348450/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:350790/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:35595/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:373053/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:5288/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:52904/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:5704/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:57091/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:7075/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:70865/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:71741/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:77615/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:78205/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:83784/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:88210/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:96572/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:147561/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:168262/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:165434/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:165197/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:163652/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:162477/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:158167/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:157315/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:149471/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:149155/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:243904/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:142717/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:135704/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:130211/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:12148/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:119267/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:110767/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:109002/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:104667/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:176080/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:253410/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:241886/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:236900/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:236634/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:235626/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:228283/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:225435/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:206470/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:205620/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:202672/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:202244/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:199322/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:194705/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:18504/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:206005/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:173690/6‑100 (MQ=1)
      tgcgtgcgcTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTCTTCTGCGGGTAAc                                                                                            >  4:18020/6‑100 (MQ=1)
              ctctACGCCCAGTAATTCCGATTAACGCTTGCACCCTCTGTATTACCGCGGCTGCTGGCACAGAGTTAGCCGGTGCTTATTCTGTCGGTAACGTCAAAAc                                                                                    <  8:28388522/97‑1 (MQ=1)
              ctctACGCCCAGTAATTCCGATTAACGCTTGCACCCTCTGTATTACCGCGGCTGCTGGCACAGAGTTAGCCGGTGCTTATCCTGTCGGTAACGTCAAAAc                                                                                    <  8:30473138/97‑1 (MQ=0)
                ttccGCCCAGTACTTCCGATTAACGCTTGCACCCTCTGTATTACCGCGGCTGCTGGCACAGAGTTAGCCGGTGCTTATTCTGTCGGTAACGTCAAAACAg                                                                                  >  7:42979768/4‑100 (MQ=0)
                             gtcagATTAACGCTTGCACCCGCTGTATTACCGCGGCTGCTGGCACAGAGTTAGCCGGTGCTTATTCTGTCGGTAACGTCAAAACAGCAAGGTATTAACt                                                                     <  8:6647216/96‑1 (MQ=0)
                                            gcgcccTCTGTATTACCGCGGCTGCTGGCACAGAGTTAGCCGGTGCTTATTCTGTCGGTAACGTCAAAACAGCAAGGTATTAACTTACGGCCCTtcctcc                                                      <  7:8518117/97‑1 (MQ=0)
                                                 ctttGTATTACCGCGGCTGCTGGTACAGAGTTAGCCGGTGCTTACTCTGTCGGTAACGTCAAAACATCAAGGTATTAACTTACTGCCCTTCCTCCTAACt                                                 <  8:48808800/97‑1 (MQ=1)
                                                     ttagtACCGCGGCTGCTGGCACAGAGGTAGCCGGTGCTGAGTATGTCGGTAACGTCAAAACAGCAAGGTAGTAACGTACTGCCCTGCCTCCCAACTGaaa                                             <  7:45189901/96‑1 (MQ=1)
                                                                   gcgggCACAGAGTTAGCCGGTGCTGATTCTGTCGGGAACGTCAAAACAGCAAGGTATTAACTTACTGCACTTCCTCCCAACTTAAAGGGCGTTACAATcc                               <  8:2540142/97‑1 (MQ=1)
                                                                               ttggCCTGTGCTTATTCTGTCGGTAACGTCAAAACAGCAAGGTATTAACTTACTGCCCTTCCTCCCAACTTAAAGTGCTTTACAATCCGAAGACcttctt                   <  8:42215849/97‑1 (MQ=0)
                                                                                         ctcaTTCTGTCGGTAACGTCAAAACAGCAAGGTATTAACTTACTGCCCTTCCTCCCAACTTAAAGTGCTTTACAATCCGAAGACCTTCTTCACACAcgcg         >  7:43433416/4‑100 (MQ=1)
                                                                                                tggcGGTAACGTCAAAACAGCAAGGTATTAACTTACTGCCCTTCCTCCCAACTTAAAGTGCTTTACAATCCGAAGACCTTCTTCACACACGCGGCATGGc  <  8:31034541/97‑1 (MQ=1)
                                                                                                  ||                                                                                                
ACCACCTACGCGCGCTTTACGCCCAGTAATTCCGATTAACGCTTGCACCCTCTGTATTACCGCGGCTGCTGGCACAGAGTTAGCCGGTGCTTATTCTGTCGGTAACGTCAAAACAGCAAGGTATTAACTTACTGCCCTTCCTCCCAACTTAAAGTGCTTTACAATCCGAAGACCTTCTTCACACACGCGGCATGGC  >  AE015451_phaM1/5312114‑5312309

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑

GATK/CNVnator alignment

N/A