Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451_phaM1 5,555,608 +GCC intergenic (‑84/+87) PP_4887 CDS (conserved protein of unknown function) ← / ← PP_4888 CDS (putative methyl‑accepting chemotaxis transducer) –/–

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM15,555,6081.G93.0% 678.0 / 60.2 251intergenic (‑84/+87)PP_4887 CDS (conserved protein of unknown function)/PP_4888 CDS (putative methyl‑accepting chemotaxis transducer)–/–
Reads supporting (aligned to +/- strand):  ref base . (0/17);  new base G (17/208);  total (17/234)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01
*AE015451_phaM15,555,6082.C93.1% 817.3 / 56.0 251intergenic (‑84/+87)PP_4887 CDS (conserved protein of unknown function)/PP_4888 CDS (putative methyl‑accepting chemotaxis transducer)–/–
Reads supporting (aligned to +/- strand):  ref base . (0/17);  new base C (13/217);  total (17/234)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*AE015451_phaM15,555,6083.C84.4% 762.2 / 152.0 251intergenic (‑84/+87)PP_4887 CDS (conserved protein of unknown function)/PP_4888 CDS (putative methyl‑accepting chemotaxis transducer)–/–
Reads supporting (aligned to +/- strand):  ref base . (0/39);  new base C (16/195);  total (17/234)
Fisher's exact test for biased strand distribution p-value = 1.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 262 total aligned reads displayed.

TGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑G‑‑‑C‑C‑‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAG  >  AE015451_phaM1/5555472‑5555727
                                                                                                                                                 |||                                                                                                                          
tgCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑g‑c‑g                                                                                                                               <  2:8539652/137‑1 (MQ=255)
tgCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑g‑c‑g                                                                                                                               <  1:119641/137‑1 (MQ=255)
tgCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑g‑c‑g                                                                                                                               <  2:14064616/137‑1 (MQ=255)
tgCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑g‑c‑g                                                                                                                               <  2:7865792/137‑1 (MQ=255)
tgCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑g‑c‑g                                                                                                                               <  2:11263520/137‑1 (MQ=255)
   aGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTTGCGAAAGGGACGCA‑‑‑TA‑‑‑g‑c‑ggcc                                                                                                                            >  1:5384628/1‑134 (MQ=255)
    ggAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑g‑c‑ggccc                                                                                                                           <  1:15267844/137‑5 (MQ=255)
              gTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTc                                                                                                              <  1:13370413/137‑1 (MQ=38)
              gTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTc                                                                                                              <  2:5883212/137‑1 (MQ=38)
                   ccTGGGCAGCGATTATTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCgt                                                                                                         <  2:15139502/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCcgcgc                                                                                                         <  2:11098211/137‑2 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCgt                                                                                                         <  2:2560470/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCgt                                                                                                         <  1:3804976/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCgt                                                                                                         <  1:5080192/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCgt                                                                                                         <  1:6963185/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCgt                                                                                                         <  2:12672933/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCACA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAACTTCCGCgt                                                                                                         <  2:15139500/137‑1 (MQ=38)
                    cTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCgtg                                                                                                        <  1:16621241/137‑1 (MQ=38)
                     tGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCgtgt                                                                                                       <  1:13584926/137‑1 (MQ=38)
                       ggCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTAt                                                                                                     <  2:4810458/137‑1 (MQ=38)
                        gCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATg                                                                                                    <  2:862095/137‑1 (MQ=38)
                        gCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATg                                                                                                    <  2:85049/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGc                                                                                                   <  1:16134758/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGc                                                                                                   <  2:14138771/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGc                                                                                                   <  1:14028565/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGc                                                                                                   <  2:1311842/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGc                                                                                                   <  1:14664042/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGc                                                                                                   <  1:7675802/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGc                                                                                                   <  1:3189748/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGc                                                                                                   <  1:5835101/137‑1 (MQ=38)
                          aGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGca                                                                                                  <  1:1884286/137‑1 (MQ=38)
                              aTTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGcaacaa                                                                                              <  2:14012099/137‑1 (MQ=38)
                              aTTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGcaacaa                                                                                              <  1:9766875/137‑1 (MQ=38)
                                tAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAAc                                                                                            <  2:10419814/137‑1 (MQ=38)
                                 aGTCGGGCGTACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑CTGCC‑C‑C‑‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                           <  2:7192077/137‑1 (MQ=38)
                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                           <  2:4996622/137‑1 (MQ=38)
                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                           <  2:7548825/137‑1 (MQ=38)
                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                           <  1:15670226/137‑1 (MQ=38)
                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                           <  1:2655370/137‑1 (MQ=38)
                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                           <  1:6861663/137‑1 (MQ=38)
                                  gTCGTGCGGACACGAGCAACTCCTTTTGTCCCTTGTCGCCTACATCATAACAACACCGCTCAGGCGTGGGAGCTGCCTCGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGC‑CTC‑‑GGCAATTTCCGCGTGTATGCAACAAACCt                                                                                          <  2:12187673/137‑1 (MQ=21)
                                   tCGGGCGGACATGAGCAACTCCTTTTTTCATNTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                         <  2:182093/137‑1 (MQ=38)
                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                         <  2:14483588/137‑1 (MQ=38)
                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                         <  2:1580764/137‑1 (MQ=38)
                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                         <  1:15008285/137‑1 (MQ=38)
                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                         <  2:4366755/137‑1 (MQ=38)
                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                         <  2:5712257/137‑1 (MQ=38)
                                    cGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAACCTgg                                                                                        <  1:16969992/137‑1 (MQ=38)
                                     gggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGGGGGAGCGGCCTTGTGTCGCGAAATGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAACCTggg                                                                                       <  2:17643966/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                      <  1:10824034/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                      <  2:593878/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                      <  1:14116001/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                      <  2:10962559/137‑1 (MQ=38)
                                            cATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTa                                                                                <  2:2205363/137‑1 (MQ=38)
                                             aTGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTAt                                                                               <  2:3543169/137‑1 (MQ=38)
                                             aTGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTAt                                                                               <  2:17136387/137‑1 (MQ=38)
                                                      tCCTTTTTTCATGTGTCGCCGAAAGCAGAGCAACACCACTAAGGTGGGGGAAAGGCCTTGTGTCGCGAAAGGACCACA‑‑‑AA‑‑‑G‑C‑GGCCC‑C‑‑GACAATTTCCGCGTGCATGAAACAAACCTGGGGCTGCTGTGCATCCCtt                                                                      <  2:5995171/137‑1 (MQ=255)
                                                      tCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAAAAAACCTGGGGCTGCTATGCAGCCCtt                                                                      <  2:5469542/137‑1 (MQ=38)
                                                          tttGTCGTGTGTCGCCGACAGCATAGAAACACCGTTAGGGTGGGGGAGTGGCTTTGTGTCGCGAAAGGGCAGAA‑‑‑TA‑‑‑G‑T‑GGC‑C‑CG‑GGCAGTTTCCGGGTGGAGGAAACAAACTTGGGGTGGCGATGCAGCCCTTTcgc                                                                  <  2:6808442/137‑1 (MQ=255)
                                                             ttCATTTGTCTCCGACAGCATAACACCACCGCCCAGGTGTGGAAGCGGCCTGGTGTCGAGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGC‑C‑CA‑GGCACTTTCCGCGTGTATGCAACAAACTTGGGGCTGCTATGCATCCCTTTAGAGac                                                               <  2:2629972/136‑1 (MQ=255)
                                                              tCATGTTTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGTGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGaca                                                              <  2:6739209/137‑1 (MQ=37)
                                                               cATGTGTCGCCGACAGCATAGAAACACCGCTAAGGTGGGGGAGCGGCCTTGTTTCGAGAAAGGGCCGAA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGGAGGAAACAAACCTGGGGCTGCTGTGCAGCCCTTTCGCGacac                                                             <  2:2883630/137‑1 (MQ=25)
                                                               cATGTGTCACCGACAGCATACCAACACCTCTCAGGTGTGGGAGCGGCCTTGTTTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGTGTGTATGCAACAAACTTGGGTCTGCTATGCAGCCCTTTCGCGacac                                                             <  1:16154212/137‑1 (MQ=25)
                                                                 tttttCGCCGACAGCATAGCAACACCGTTCAGGTGTGGGAGTGGCTTGGTGTCGAGAAAGGGCCGAA‑‑‑TAGTGG‑C‑C‑‑‑C‑C‑‑GGCAGTTTCCGCGGGGAGGCAAAAAACTTGGGGCTGCTATGCAGCCCTTTCGCGACACaa                                                           <  2:1534130/133‑1 (MQ=19)
                                                                 tgtgTCGCTGAAAGCAGAGAAACACTACTCAGGTGGGGGAGTGGCCTTGTGCAGCGAAAGGGCCGCA‑‑‑TA‑‑‑T‑C‑GGC‑C‑CT‑GGCAGTTTCTGAGTGGGTGCAACAAACCTGGGGCTGCCATGCAGCCCTTTCGCGACACaa                                                           <  2:14909809/137‑1 (MQ=255)
                                                                 tgtgTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCTTTGTGCCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GTG‑G‑‑‑C‑CCTGGAAATTCCCGCGGGTATGAAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACaa                                                           <  2:455704/137‑1 (MQ=255)
                                                                 tgtgTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACaa                                                           <  2:10874071/137‑1 (MQ=38)
                                                                 tgtgTCGCCGAAAGCATAAAAACACCGCTCGGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑GA‑‑‑G‑C‑GGCCC‑C‑‑GGCAGTTTCCGTGTGTAGGCAACAAACTTGGGGCTGTTGTGCAGCCTTTTCGCGACACaa                                                           <  2:9014185/137‑1 (MQ=21)
                                                                           aCAGCAGAGAAACACCGCTCAGGTGTGGGGGTGGCTTTGTGTCGCGAAGGGGCCGCG‑‑‑TA‑‑‑G‑T‑GGC‑C‑CT‑GGCAATTTCTGTGTGTATGCAACAACCCTGGGGCTGCTATGCGGCCTTTTCGCGACACAAGGCCCACaaa                                                 <  2:637122/137‑1 (MQ=255)
                                                                                aTAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGAGAAGGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCACGTGTATGCAACAACCCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGAcc                                            <  2:5319930/137‑1 (MQ=37)
                                                                                aTAACAACACCGCTCAGGGGTGGGAGTGGCTTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GACAATTTCCGTGTGTATGCAAAAAACTTGGGGTTGCTATGCAGCCTTTTCGAGACACAAGGCCCACAAAGGAcc                                            <  2:14551882/137‑1 (MQ=21)
                                                                                  agaaaCACCGCTAAGGTGGGGGAGTGGCCTTGTGTCGAGAAAGGGCCGAA‑‑‑GA‑‑‑G‑C‑GGCCC‑C‑‑GGAAGTTTCCGTGTGGATGCAACAAACTGGGGGTGGCTATGAAGCCTTTTCGCGACACAAGGCCCACAAAGGACcgc                                          <  2:10065939/134‑1 (MQ=255)
                                                                                  aGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑GA‑‑‑G‑C‑GGCCC‑C‑‑GGAAATTTCCGTGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACcgc                                          <  2:8733457/137‑1 (MQ=25)
                                                                                       ccccGCTCAGTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑GA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTAGGCAACAACCCTGGGGCTGCTATGCAGCCCTTTCTCGACACAAGGCCCACAAAGGACCGCGCAgg                                     <  2:16275260/135‑1 (MQ=25)
                                                                                        acCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGc                                    <  2:8716800/137‑1 (MQ=38)
                                                                                        acCGCCCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCACCAAACCTTGGTCTGCTATGCATCCCTTCCGCGACACAAGGCCCACAAAGGACCGCGCAGGc                                    <  2:17163649/137‑1 (MQ=25)
                                                                                        acCCCTCAGTTCTGGAAGCGGCTTTGCGTCCCAAAAGGGCCACG‑‑‑TA‑‑‑C‑C‑GCCCC‑C‑‑GCCAATTTCCCCTTGCCTTCAACAAACCTTCGACTGCTCTGCAGCCCCTTTCCAACACAAGGCCCACAAAGGACCGCGCAGGc                                    <  2:13203908/137‑1 (MQ=255)
                                                                                             taaGGTGGGGGAGCGACCTGTTGTCGCGAAAGGACAGCA‑‑‑TAGTGG‑C‑C‑‑‑C‑C‑‑GGCAATTCCTGCGTGCATACAACAAACCTGGGGTTGCGAGAAAACCCTTTCGAAACACAAGGCCCACAAAGGACCGCGCAGGCTTCAg                               <  1:12745727/135‑1 (MQ=255)
                                                                                             tCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAACCCTGGGGCTGCTATACAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAg                               <  2:4074819/137‑1 (MQ=37)
                                                                                               aGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                      <  1:10499644/80‑1 (MQ=18)
                                                                                                  tgtgGGAGTGGCCTGGTGTTGTGAAAGGGTCACA‑‑‑TAGTGG‑C‑C‑‑‑C‑C‑‑GGCACTTTCTGCGTGTATGAAACAACCTTGGGGCTGCGATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                          <  2:16272616/137‑1 (MQ=255)
                                                                                                  tgtgGGAGCGGCTTTTTGTCGCGAAGGGGCCGCA‑‑‑GG‑‑‑G‑C‑GGC‑C‑CT‑GGAAGTTTCTGTGGGGGTGCAACAACCTTGGGGCGGTGGTGAGACCCTTTTGCAACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                          <  1:4108519/137‑1 (MQ=255)
                                                                                                  tgtgGGAGCGGCTTTGTTTCGCGAAAGGGCCGAAGGGTG‑‑‑G‑C‑C‑‑‑C‑T‑‑GGCACTTTCCGCGGGGATACAACAAACCTGGGGCTGCTATGCATCCTTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                          <  2:4624486/137‑1 (MQ=255)
                                                                                                  tgtgGGAGCGGCCTTGTGTCGCGAAGGTGCCGCA‑‑‑GA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGGAGGCAACAAACTGGGGGTTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                          <  2:3860716/137‑1 (MQ=25)
                                                                                                  tgtgGGAGCGGCCTTGTGTCGCGAAAGGTCCGCA‑‑‑GA‑‑‑G‑C‑GGC‑C‑CT‑GGAAATTTCCGCGTGTATGAAAAAAACCTGGGGCTGTGGGGAAGCCTTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                          <  1:17297017/137‑1 (MQ=21)
                                                                                                  tgtgGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTGTGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                          <  2:12583760/137‑1 (MQ=37)
                                                                                                  tgtgGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                          <  2:16582406/137‑1 (MQ=38)
                                                                                                  tgtgGGAGCGGCCTTGTGTCGCAAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                          <  2:3320748/137‑1 (MQ=38)
                                                                                                  tgtgGGAGCGGCCTTGTGTCGCAAAAGGGCCGAA‑‑‑GA‑‑‑G‑C‑GGC‑C‑CT‑GGCAATTTCCGCGGGTGTGCAACAACCCTTGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                          <  2:7348680/137‑1 (MQ=25)
                                                                                                  tgggggAGCGGCCTTGTGTCGCAAAAGGTCCACA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATACAACAACCCTGGGGCTGCTATTCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                          <  2:1606226/134‑1 (MQ=255)
                                                                                                  tgggggAGAGGCTTTGTGGAGAGAAAGGGCCGCA‑‑‑GA‑‑‑G‑G‑GGC‑C‑CA‑GGAAATTTCTGTGGGGAGGAAAAAAACTGGGGGTTGCGATGCAGCCCTTTCGAGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                          <  1:11313491/134‑1 (MQ=255)
                                                                                                         gCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGAAACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTgg                   >  1:8351429/1‑137 (MQ=38)
                                                                                                             ccTTGTGTCGCGAAAGGGCCGCA‑‑‑GA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAACCTTGGGCTGCGATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCa               <  2:16662704/137‑1 (MQ=38)
                                                                                                               ttGTGTCGAGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGAAGTTTCCGTGTGGAGGAAAAAAACCTGGGGCTGCTATGCAGCCTTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAgg             <  2:11576111/137‑1 (MQ=255)
                                                                                                                  tttcacgAAAGGCCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTGTGAAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTg          <  2:17426357/132‑1 (MQ=255)
                                                                                                                          aaGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCCGGCCCGGAACTGGCCCAAGCTGGCGCGCAg  >  2:1369003/1‑137 (MQ=255)
                                                                                                                          aaGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGAACGCGTAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAg  >  2:2514492/1‑137 (MQ=255)
                                                                                                                          aaGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGACC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAACCCGGGGTTGCTATGCAACCCTTTCGTGACACAAGGCCCACCAAGGACCGCGCAGGCTTCCAGCCCGGAACTTGCCCAAGCTGGCGCGCAg  >  2:7772652/1‑137 (MQ=255)
                                                                                                                          aaGGGCCGCA‑‑‑CA‑‑‑G‑C‑GGCCC‑C‑‑GGAAATTTCCGCGTGTATGCCAAAAACCTGGGGGTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAg  >  2:4204023/1‑137 (MQ=255)
                                                                                                                          aaGGGACGCA‑‑‑TA‑‑‑G‑C‑GGCCC‑C‑‑GGCAATTTCCGCGTGTATGCAACAAAACTGGGGCTGCTATGCAGCCCTTTCGCGAAACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCaa  >  2:16408117/1‑136 (MQ=255)
                                                                                                                                                 |||                                                                                                                          
TGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑G‑‑‑C‑C‑‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAG  >  AE015451_phaM1/5555472‑5555727

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