Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451_phaM1 5,681,415 +CGGG intergenic (‑133/+71) PP_4986 CDS (Channel protein_ hemolysin III family) ← / ← PP_4987 CDS (putative Chemotaxis protein) –/–

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM15,681,4151.C98.4% 474.5 / ‑3.1 127intergenic (‑133/+71)PP_4986 CDS (Channel protein_ hemolysin III family)/PP_4987 CDS (putative Chemotaxis protein)–/–
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base C (83/42);  minor base A (2/0);  total (85/42)
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.64e-01
*AE015451_phaM15,681,4152.G99.2% 464.5 / NA 127intergenic (‑133/+71)PP_4986 CDS (Channel protein_ hemolysin III family)/PP_4987 CDS (putative Chemotaxis protein)–/–
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base G (84/42);  minor base A (1/0);  total (85/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*AE015451_phaM15,681,4153.G99.2% 480.2 / ‑6.3 127intergenic (‑133/+71)PP_4986 CDS (Channel protein_ hemolysin III family)/PP_4987 CDS (putative Chemotaxis protein)–/–
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base G (84/42);  minor base A (1/0);  total (85/42)
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.99e-01
*AE015451_phaM15,681,4154.G100.0% 492.0 / NA 127intergenic (‑133/+71)PP_4986 CDS (Channel protein_ hemolysin III family)/PP_4987 CDS (putative Chemotaxis protein)–/–
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (85/42);  total (85/42)

Only 100 of 148 total aligned reads displayed.

ACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGC‑‑‑‑GCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTT  >  AE015451_phaM1/5681279‑5681529
                                                                                                                                         ||||                                                                                                                  
aCATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACcgc                                                                                                                        <  2:15973974/137‑1 (MQ=255)
aCATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACcgc                                                                                                                        <  2:15218409/137‑1 (MQ=255)
aCATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACcgc                                                                                                                        <  1:8950050/137‑1 (MQ=255)
aCATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACcgc                                                                                                                        <  2:2753926/137‑1 (MQ=255)
  aTGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACcgccg                                                                                                                      <  2:18938586/137‑3 (MQ=255)
   tGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACcgccgg                                                                                                                     <  2:11232808/137‑4 (MQ=255)
     cAAACACATACCTCCTGAGGTGGGCTGGGTTCCTTTTGTGAGGCTGATGGCTTCGTTCTTTGTGGGAGCGGTCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACcgccgggg                                                                                                                   <  1:12253562/137‑6 (MQ=38)
     cAAACACATACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACcgccgggg                                                                                                                   <  2:1044508/137‑6 (MQ=255)
     cAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACcgccgggg                                                                                                                   <  2:9421807/137‑6 (MQ=255)
     cAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACcgccgggg                                                                                                                   <  1:71633/137‑6 (MQ=255)
     cAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACcgccgggg                                                                                                                   <  1:6634362/137‑6 (MQ=255)
     cAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACcgccgggg                                                                                                                   <  1:18756196/137‑6 (MQ=255)
     cAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACcgccgggg                                                                                                                   <  2:81114/137‑6 (MQ=255)
     cAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACcgccgggg                                                                                                                   <  2:6974063/137‑6 (MQ=255)
     cAAACACAGACCCCCTTAGGTGGGCTGGGTTCCGGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACcgccgggg                                                                                                                   <  1:6634363/137‑6 (MQ=38)
                     taGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  1:15433753/136‑1 (MQ=37)
                     gAGGTGGGCTGGGTTCCTTTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTTTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  1:7265336/137‑1 (MQ=37)
                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTNGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  1:15752313/137‑1 (MQ=37)
                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  2:11710178/137‑1 (MQ=37)
                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  1:1386644/137‑1 (MQ=37)
                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  2:1732484/137‑1 (MQ=37)
                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  1:4541854/137‑1 (MQ=37)
                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  2:17232242/137‑1 (MQ=37)
                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  2:1691980/137‑1 (MQ=37)
                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  1:5069490/137‑1 (MQ=37)
                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  1:5189038/137‑1 (MQ=37)
                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  2:1691979/137‑1 (MQ=37)
                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  1:12856803/137‑1 (MQ=37)
                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  2:9403935/137‑1 (MQ=37)
                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  2:7608282/137‑1 (MQ=37)
                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  2:11380613/137‑1 (MQ=37)
                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  2:1948341/137‑1 (MQ=37)
                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  1:15432583/137‑1 (MQ=37)
                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  2:3104068/137‑1 (MQ=37)
                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  1:18498450/137‑1 (MQ=37)
                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  1:9474214/137‑1 (MQ=37)
                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  2:5104131/137‑1 (MQ=37)
                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  1:16247892/137‑1 (MQ=37)
                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  2:609816/137‑1 (MQ=37)
                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  1:15981029/137‑1 (MQ=37)
                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  1:15926123/137‑1 (MQ=37)
                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  2:10220897/137‑1 (MQ=37)
                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  2:10670686/137‑1 (MQ=37)
                     gAGGTGGGCTGGGTTCCTGTTGTGATCCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  1:19306309/137‑1 (MQ=37)
                      aGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTAt                                                                                                  <  1:15820146/137‑1 (MQ=37)
                                                                                        aaaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAg                                >  2:8674868/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCCTGAACCGCAGGGGCTGCTGCGCAGTCATATCGTGAAAACAGGCCGATCCCCCTCCAGACCACAAGGTCCTGTGCCTCCCATCCTAACAAGGGTCAgg                               >  2:19268849/1‑137 (MQ=255)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGTTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:9137918/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCCGACCACCCGGCCCTGCGCCTCCCATCCTAACCAAGGTCAgg                               >  1:10531913/1‑137 (MQ=37)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  2:15343773/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  2:16713467/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  2:14266771/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  2:14082551/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  2:14025463/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  2:13922731/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  2:13180562/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  2:18454818/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  2:4183295/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:10365722/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:7790275/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:4232160/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:4812744/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:5044632/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:16063580/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:5440188/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:3669175/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:14486436/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:7532166/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:8176260/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:8229096/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:8359285/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:19096257/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:19087366/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:19066603/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCACTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:8307219/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGAAACACGGCCGCTCCCACTCCAGACCACCAGGCCCTGTGCCTCACATCCTAACCAAGGTCAgg                               >  1:8737752/1‑137 (MQ=255)
                                                                                         aaGGGGCGCGAAGCGGCCCCAGGATTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCATCACATCCTAACCAAGGTCAgg                               >  1:1850378/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCCCCAAGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  2:19500400/1‑137 (MQ=38)
                                                                                         aaGGGGCGCGAAGCGGCACAAGGCTTTCATTGCCACTGCATGAACAGCAGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  2:3723602/1‑137 (MQ=38)
                                                                                              gcgcGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTACCCAAGGTCAGTTTGCa                          >  1:11426125/1‑137 (MQ=38)
                                                                                              gcgcGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCa                          >  1:9024421/1‑137 (MQ=38)
                                                                                                    aGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCcg                    >  1:18058707/1‑137 (MQ=38)
                                                                                                      cGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCcgcg                  >  1:5818960/1‑137 (MQ=38)
                                                                                                          cccAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCAc              >  2:4787711/1‑137 (MQ=38)
                                                                                                          cccAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCAc              >  2:11678486/1‑137 (MQ=38)
                                                                                                           ccAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACg             >  1:103360/1‑137 (MQ=38)
                                                                                                           ccAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACg             >  1:3851930/1‑137 (MQ=38)
                                                                                                            cAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGc            >  2:4765817/1‑137 (MQ=38)
                                                                                                                 tttCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAg       >  1:13639227/1‑137 (MQ=38)
                                                                                                                      ttGCCACTTCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTCACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt  >  1:5439214/1‑137 (MQ=38)
                                                                                                                      ttGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCCCGCACCAGCCttt  >  2:18050704/1‑137 (MQ=38)
                                                                                                                      ttGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt  >  1:18332628/1‑137 (MQ=38)
                                                                                                                      ttGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt  >  2:176228/1‑137 (MQ=38)
                                                                                                                      ttGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt  >  2:6450614/1‑137 (MQ=38)
                                                                                                                      ttGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt  >  1:13286761/1‑137 (MQ=38)
                                                                                                                      ttGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt  >  2:127571/1‑137 (MQ=38)
                                                                                                                      ttGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt  >  1:10876289/1‑137 (MQ=38)
                                                                                                                      ttGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt  >  2:10523067/1‑137 (MQ=38)
                                                                                                                      ttGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGAAACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt  >  2:6949754/1‑137 (MQ=38)
                                                                                                                      ttGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCACTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt  >  1:6213315/1‑137 (MQ=38)
                                                                                                                                         ||||                                                                                                                  
ACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGC‑‑‑‑GCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTT  >  AE015451_phaM1/5681279‑5681529

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