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.G92.2% 760.9 / 78.5 278intergenic (‑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/21);  new base G (11/236);  total (11/267)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-01
*AE015451_phaM15,555,6082.C92.5% 891.3 / 71.8 278intergenic (‑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/21);  new base C (11/246);  total (11/267)
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,555,6083.C85.6% 852.5 / 154.7 278intergenic (‑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/40);  new base C (11/227);  total (11/267)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 294 total aligned reads displayed.

TGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑G‑‑‑CC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAG  >  AE015451_phaM1/5555472‑5555727
                                                                                                                                                 |||                                                                                                                        
tgCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑g‑c‑g                                                                                                                             <  1:9150785/137‑1 (MQ=255)
tgCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑g‑c‑g                                                                                                                             <  1:10117573/137‑1 (MQ=255)
tgCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑g‑c‑g                                                                                                                             <  2:2837638/137‑1 (MQ=255)
    ggAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑g‑c‑ggccc                                                                                                                         <  2:11174582/137‑5 (MQ=255)
                cGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCcg                                                                                                            <  1:3270068/137‑1 (MQ=38)
                cGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCcg                                                                                                            <  2:942978/137‑1 (MQ=38)
                cGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCcg                                                                                                            <  1:1270916/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCgt                                                                                                         <  1:10230755/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCgt                                                                                                         <  1:8106876/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCgt                                                                                                         <  2:17632760/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCgt                                                                                                         <  2:18012241/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCgt                                                                                                         <  1:14537057/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCgt                                                                                                         <  1:14124466/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCgt                                                                                                         <  2:8552340/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCgt                                                                                                         <  2:8913862/137‑1 (MQ=38)
                        gCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATg                                                                                                    <  2:15734436/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                   <  2:12310166/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                   <  1:11201551/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                   <  1:5163168/137‑1 (MQ=38)
                          aGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGca                                                                                                  <  2:422050/137‑1 (MQ=38)
                              aTTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGcaacaa                                                                                              <  2:15606153/137‑1 (MQ=38)
                              aTTAGTCGGGCGGACATGAGCAACCCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGcaacaa                                                                                              <  2:16947989/137‑1 (MQ=38)
                               ttAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                   >  1:1662503/1‑131 (MQ=38)
                               ttAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                   <  2:1662503/131‑1 (MQ=38)
                               ttAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACaaa                                                                                             >  2:9712553/1‑137 (MQ=38)
                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                           <  1:8913921/137‑1 (MQ=38)
                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                           <  2:7128336/137‑1 (MQ=38)
                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                           <  2:231460/137‑1 (MQ=38)
                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                           <  1:2868859/137‑1 (MQ=38)
                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                           <  1:1722163/137‑1 (MQ=38)
                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTGGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                           <  1:7790921/137‑1 (MQ=38)
                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                         <  2:9289168/137‑1 (MQ=38)
                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                         <  2:14984636/137‑1 (MQ=38)
                                    cGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgg                                                                                        <  1:13610137/137‑1 (MQ=38)
                                    cGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgg                                                                                        <  2:7712137/137‑1 (MQ=38)
                                     gggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTggg                                                                                       <  2:5584445/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                      <  2:16204980/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                      <  2:4612775/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                      <  1:2349252/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                      <  1:9712553/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                      <  2:16102956/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                      <  1:17864160/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                      <  2:8351/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                      <  1:13122712/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                      <  1:11712647/137‑1 (MQ=38)
                                         ggACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGgctg                                                                                   <  1:572903/137‑1 (MQ=38)
                                            cATGAGCAACTTCTTTTGTCATGTGTCGCCGACAGCGTAGCAACACCGCTCGGGTGTGGGAGCGGCTTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTa                                                                                <  2:12606733/137‑1 (MQ=37)
                                                 gCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAg                                                                           <  2:18141697/137‑1 (MQ=38)
                                                   aaCTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑T‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGcc                                                                         <  2:18764647/137‑1 (MQ=38)
                                                      tCCTTTTGTCGTGTTTCTCCGACAGCATAGCAACACCGCTCGGGTGTGGGAGCGGCCTTGTGTCGAGAAGGGGCCGCA‑‑‑TATCGG‑C‑C‑‑‑CC‑GGCAATTCCCGCGGGTATGCAACAACCCTGGGGCTGCTATGCAGCCCtt                                                                      <  2:12099080/137‑1 (MQ=25)
                                                      tCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAGGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAGTTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCtt                                                                      <  2:1397225/137‑1 (MQ=38)
                                                      tCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCtt                                                                      <  2:13268172/137‑1 (MQ=38)
                                                             ttCATGTGTCGCCGACGACATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑GA‑‑‑G‑C‑GGCCCC‑GGCACTTTCCGCGTGTGTGCAACAAACCTGGGGCTGCTATTCACCCCTTTCGCGac                                                               <  2:15515124/136‑1 (MQ=25)
                                                               catttttCGCCGAAAGCATAGCAACACCGTTCAGGTGTGGGAGCGGTCTTGTGTCGCGAAAGGGCCGCA‑‑‑TG‑‑‑G‑G‑GGC‑CCTGGCAATTTCCGCGTGTGTGCAACAAACTTGGGGCTGCTATGCAGCCCTTTCGCGacac                                                             <  2:14344875/131‑1 (MQ=25)
                                                               cATGTTTCGCCGACAGCATATCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAGGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGacac                                                             <  2:17290555/137‑1 (MQ=37)
                                                               cATGTTTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGTGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCGATGCAGCCCTTTCGCGacac                                                             <  2:2680562/137‑1 (MQ=37)
                                                               cATGTGTCGCTAACAGCATAGCAACACTGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑T‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGTTGCTATGCAGCCCTTTCGCGacac                                                             <  2:10232581/137‑1 (MQ=37)
                                                               cATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGTGGCCTTGTGTCGTGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCTGCGTGTATACAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGacac                                                             <  2:13734205/137‑1 (MQ=38)
                                                               cATGTGTCGCCGACAGCACAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGTTATGCAGCCCTTTCGCGacac                                                             <  2:6179132/137‑1 (MQ=37)
                                                                 tgtgTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTTGCGAAGGGTCCGCA‑‑‑TA‑‑‑G‑T‑GGCCCC‑GGCAATTTCCGAGGGTATGCAACAAACTTGGGGCTGCGATGCAGCCCTTTCGCGACACaa                                                           <  2:4122140/137‑1 (MQ=37)
                                                                 tgtgTCGCCGACAGCACAGCAACCCCGTTCGGGTGTGGAAGCGGCTTTGTTTCGCGAAGGGTCCGCA‑‑‑TATT‑GTC‑C‑‑‑TC‑GGCAGTTTCCGCGGGTGTACAGTAAACCTGGGGTTGCTGTGCAGCCCTTTCGCGACACaa                                                           <  2:7656122/137‑1 (MQ=255)
                                                                      cGCCAACAGCATAGCAACACCCCTCAGGTGTGGGGGCGGCCTTGTGCCGGAAAGGGGCCGCAGAGTG‑‑‑G‑C‑C‑‑‑CC‑GGCAGTTTCCGTGTGTGTGAAAAAAACTTGGGGCTGCTATGCAGCCTTTTCGCGACACAAGGccc                                                      <  2:6206024/137‑1 (MQ=255)
                                                                           aaagcaaagcaaCACCGCTCAGGTGTGGAAGGGACTTGGGGTAGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCGACGGGGAGGAAAAAAACCGGGGGCTGCGATGAAGCCCTTTCGCGAAACAAGGCCCACaaa                                                 <  2:15950442/135‑1 (MQ=255)
                                                                                atggCAACACCGCTAAGGTGGGGGAGCGGCCTTGTTTCGCGAAAGGGCCGCAGGGTG‑‑‑G‑C‑C‑‑‑CC‑GGCAATTTCCGCGTGTATGCAACAAACTTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGAcc                                            <  2:2075385/134‑1 (MQ=37)
                                                                                agagaaaCACCTCTCGGGTGGGGGAGGGGCTTTGTGTCGCGAAGGGTCCGCG‑‑‑GG‑‑‑G‑C‑GGC‑CCTGGAAGTTCCCGTGTGTATGAAAAAAACTTGGGGCTGCTATGCAGCCTTTTCGCGACACAAGGCCCACAAAGGAcc                                            <  1:13363630/132‑1 (MQ=255)
                                                                                agagaaaCACCGCTCAGGTGTGGGAGCGGCCTTGTTTCGCGAAAGTGCCGCA‑‑‑GA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGAAAAAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGAcc                                            <  2:7902684/132‑1 (MQ=25)
                                                                                aTAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCTCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAGTTTCCGTGTGTGTGCAACAAACCTGGGGCGGCTAGGCAGCCCTTTCGCGACACAAGGCCCACAAAGGAcc                                            <  2:5784301/137‑1 (MQ=25)
                                                                                aTAGCAACACCGCTAAGGTGTGGGAGCGGCCTTGTTTCGCGAAGGGGCTGCA‑‑‑GA‑‑‑G‑C‑GGC‑CCTGGCAATTTCCGTGGGGGGGAAACAAACTTGGGGCGGCTATGCGGCCTTTTCGGGACACAAGGCCCACAAAGGAcc                                            <  2:6371833/137‑1 (MQ=255)
                                                                                aTAGCAACACCCCTCAGGTGGGGGAGTGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TATCGG‑C‑C‑‑‑CC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCGGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGAcc                                            <  2:3627116/137‑1 (MQ=37)
                                                                                aTAGCAACACCCCCCAGGTGTGGGAGCGACCTTGTGCCGCAAAAGGACCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGAAATTCCCGCGGGTATGCAACAAACCTGGGGCGGCTATGCAGCCCTTTCACGACACAAGGCCCACAAAGGAcc                                            <  2:12972521/137‑1 (MQ=255)
                                                                                aTAACAACACCGCTCAGGTGTGGGAGCGGCCTTGGGTCGCGAAAGGGCCGCG‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGAcc                                            <  2:10964588/137‑1 (MQ=37)
                                                                                  agaaaCACCGCTCGGGTGTGGGAGCGGCTTTGTGTCGAGAAGGGGCCGCA‑‑‑GA‑‑‑G‑C‑GGCCCC‑GGCAGTTTCCGTGGGGATGAAAAAAACTTGGGGCTGCGATGAAGCCCTTTCGCGACACAAGGCCCACAAAGGACcgc                                          <  2:10467218/134‑1 (MQ=16)
                                                                                  aacAACACCGCTAAGGGGGGGGAGCGGCCTTGTGTCGAGAAGGGGCCGCAGAGTG‑‑‑G‑C‑C‑‑‑CC‑GGCAGTTTCTGTGTGGAGGAAAAAAACCTGGGGTTGCGGGGAAGCCTTTTCGCGACACAAGGCCCACAAAGGACcgc                                          <  2:9049202/135‑1 (MQ=255)
                                                                                  aaaaaCACCGCTAAGGTGGGGGAGCGGCTTTGTGTTGCGAAAGGGCAGCA‑‑‑GA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGTGTGTGGGCAACAAACTTGGGGCTGTTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACcgc                                          <  2:3362179/134‑1 (MQ=21)
                                                                                  aGCAACACCGCTCAGGTGGGGGATCGGCCTTGTGTAGCGAAAGGGCCGAA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAAAAAACCTGGGGCTGCTAGGCGGCCTTTTCGCGACACAAGGCCCACAAAGGACcgc                                          <  2:10051444/137‑1 (MQ=255)
                                                                                     aaCACCGTTCGGGGGTGGGAGTGGCTTTGTGTCGGGAAAGGGCTGCA‑‑‑GA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGGGGGCAACAAATTGGGGGTTGCGATGCAGCCCTTTCGTGACACAAGGCCCACAAAGGACCGCGCa                                       <  1:15258028/137‑1 (MQ=255)
                                                                                       cccctcTCGGGTTTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑T‑GGCCCC‑GGCAATTTCCGTGTGTGTGCAACAAACCTGGGGTTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAgg                                     <  2:690413/132‑1 (MQ=25)
                                                                                        acCGCTCGGGTGTGGGAGTGGCTTTGTGTTGCGAAAGGGCCGCA‑‑‑TG‑‑‑G‑C‑GGC‑CCTGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGc                                    <  2:655177/137‑1 (MQ=37)
                                                                                          cGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTggg                                                                                       >  2:3107742/1‑84 (MQ=18)
                                                                                             tcgggTGGGGGAGCGGCTTTGTTTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACTTGGGGCTGCGGTGCAGCCCTTTCGCGAAACAAGGCCCACAAAGGACCGCGCAGGCTTCAg                               <  2:1981446/134‑1 (MQ=255)
                                                                                             taaGGTGTGGGAGCGGCCTTTTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGGAGGCAAAAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAg                               <  2:6756902/135‑1 (MQ=21)
                                                                                             ccAGGTGTGGGAGCGGCCTTGTTCCGCGAAAGGGCCGAA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTAGGAAACAAACCGGGGGTTGCTAGGCACCCCTTTCGCGAAACAAGGCCCACAAAGGACCGCGCAGGCTTCAg                               <  2:3911670/136‑1 (MQ=255)
                                                                                             ccAGGTGGGGGAGCGGCTTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTCCCGCGTGTGTGAAAAAAACCTGGGGTGGCTGTGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAg                               <  1:2017607/136‑1 (MQ=25)
                                                                                                  tgtgGGAGTGGCTTTGTTTTGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAGTTTCCGTGGGTATGCAACAAACCGGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                          <  1:18453142/137‑1 (MQ=25)
                                                                                                  tgtgGGAGCGGCCTTGTGTCGCGAAGGGGCCGCA‑‑‑TAGTGG‑C‑C‑‑‑CC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCTTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                          <  2:17314359/137‑1 (MQ=38)
                                                                                                  tgtgGGAGCGGCCTTGTGTCGCAAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GTCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGAAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                          <  1:9034022/137‑1 (MQ=38)
                                                                                                  tgtgGGAGCGGCCTTGTGTCGCAAAAGGACCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GACAATTTCCGCGTGTATGCAACAACCCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                          <  2:10692372/137‑1 (MQ=38)
                                                                                                  tgggggAGCGGTCTTGTGTCGAGAAGGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAGTTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                          <  2:15138254/134‑1 (MQ=37)
                                                                                                   ttgGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGTCTGTTATGCATCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCgg                         <  2:17478410/136‑1 (MQ=37)
                                                                                                   gtgGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑T‑C‑GGCCCC‑GGCAATTTCCGCGTGTATACAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCgg                         <  2:5470737/137‑1 (MQ=38)
                                                                                                       gggCGGCCGTGTGTCGCGAAAGGGCCACA‑‑‑TA‑‑‑G‑CTGCC‑CC‑GACAATTTCCGCGTGTGTGAAACAAACCGGGTGCTGTTATGCAGCCTTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACt                     <  1:15943428/135‑1 (MQ=255)
                                                                                                       gAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑GA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGTGTGTATGAAACAAACTTGGGGTTGCTGTGCAGCCTTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACt                     <  2:6288510/137‑1 (MQ=25)
                                                                                                                  tttCGCGAAGGGGCCGCA‑‑‑GG‑‑‑G‑C‑GTCCCC‑GGCACTTTCCGTGTGTATGAAACAAACTTGGGGTTGCTATGCAGCCTTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTg          <  1:4349788/135‑1 (MQ=255)
                                                                                                                  tgtAGCGAAAGGGCCACA‑‑‑TA‑‑‑G‑C‑GGC‑CCTGGCAATTTCCGCGTGGATGCAACAAACTTGGGGCTGCTAGGAAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTg          <  2:11380129/137‑1 (MQ=255)
                                                                                                                    tCGCGAAGGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATACAACAAACCTGGGGTTGCTGGGCAGCCTTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGgc        <  2:14562415/137‑1 (MQ=255)
                                                                                                                        gAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAGTTTCCGCGTGTATGCAACAAACCTGGGTCTGCTATGCAACCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGgcgcgc    <  2:7450597/137‑1 (MQ=255)
                                                                                                                         aaaGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAACCCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCa   <  2:4471974/137‑1 (MQ=255)
                                                                                                                          aaaggCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCCGGGGCTGCTCTGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAAGCTTGCGCGCAg  >  2:18359889/4‑137 (MQ=255)
                                                                                                                          aaGGGCCGCA‑‑‑TA‑‑‑G‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAAACTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCACGCAGGCTTCCGGCCCGGAACTGGCCCCGGCTGGCGCGCAg  >  2:10837388/1‑137 (MQ=255)
                                                                                                                          aaGGGCCACA‑‑‑TA‑‑‑G‑C‑GGCCCC‑AGCAATTTCCGCGTGTATGCCACAAACCTGGGGCTGCTCTGCAGCCCTTTCGCGACACACGGCCCACAAAGGACCGCGCAGGCTTCCGGCCCTGAACTGGCCCAGGCTGGCGCGCAg  >  2:1476296/1‑137 (MQ=255)
                                                                                                                                                 |||                                                                                                                        
TGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑C‑G‑‑‑CC‑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: ‑