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.G90.3% 740.2 / 93.6 266intergenic (‑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/25);  new base G (7/225);  total (7/259)
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,6082.C90.5% 840.9 / 86.1 266intergenic (‑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/25);  new base C (5/234);  total (7/259)
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.C80.0% 774.4 / 208.0 266intergenic (‑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/53);  new base C (6/206);  total (7/259)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 273 total aligned reads displayed.

TGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑G‑‑‑CC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAG  >  AE015451_phaM1/5555472‑5555727
                                                                                                                                                |||                                                                                                                        
tgCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑g‑‑‑c‑g                                                                                                                             <  2:14603725/137‑1 (MQ=255)
tgCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAACATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑g‑‑‑c‑g                                                                                                                             <  2:4982207/137‑1 (MQ=255)
    ggAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑g‑‑‑c‑ggccc                                                                                                                         <  2:6273480/137‑5 (MQ=255)
    ggAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑g‑‑‑c‑ggccc                                                                                                                         <  2:11980759/137‑5 (MQ=255)
      aaaTTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑g‑‑‑c‑ggccccg                                                                                                                       <  1:11597054/137‑7 (MQ=255)
      aaaTTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑g‑‑‑c‑ggccccg                                                                                                                       <  1:9433781/137‑7 (MQ=255)
              gTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTc                                                                                                              <  2:16649389/137‑1 (MQ=38)
                cGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCcg                                                                                                            <  1:12013006/137‑1 (MQ=38)
                 ggCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCcgc                                                                                                           <  2:14065069/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATNTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCgt                                                                                                         <  1:15263227/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCgt                                                                                                         <  1:2836237/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCgt                                                                                                         <  1:10253029/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCgt                                                                                                         <  1:7914000/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCgt                                                                                                         <  2:12999478/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCgt                                                                                                         <  2:13764433/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCgt                                                                                                         <  2:15055900/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCgt                                                                                                         <  2:18399221/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCgt                                                                                                         <  2:3114110/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCgt                                                                                                         <  2:4579757/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCgt                                                                                                         <  2:5459637/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCgt                                                                                                         <  1:2129009/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCAGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCgt                                                                                                         <  1:3579851/137‑1 (MQ=38)
                     tGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCgtgt                                                                                                       <  2:16822654/137‑1 (MQ=38)
                     tGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCgtgt                                                                                                       <  2:9221960/137‑1 (MQ=38)
                     tGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCgtgt                                                                                                       <  1:7933542/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCNTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                   <  2:665927/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                   <  1:14993769/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                   <  2:2504568/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                   <  2:7106750/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                   <  2:16005922/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                   <  1:13835926/137‑1 (MQ=38)
                          aGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGca                                                                                                  <  1:5179186/137‑1 (MQ=38)
                            cGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGcaac                                                                                                <  2:16962246/137‑1 (MQ=38)
                              aTTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGcaacaa                                                                                              <  1:8942799/137‑1 (MQ=38)
                              aTTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGcaacaa                                                                                              <  1:9419740/137‑1 (MQ=38)
                                tAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAAc                                                                                            <  2:15283664/137‑1 (MQ=38)
                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                           <  1:12082595/137‑1 (MQ=38)
                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                           <  2:10294511/137‑1 (MQ=38)
                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                           <  2:18223768/137‑1 (MQ=38)
                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                         <  2:16699189/137‑1 (MQ=38)
                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                         <  2:15512896/137‑1 (MQ=38)
                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                         <  1:3332208/137‑1 (MQ=38)
                                    cGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgg                                                                                        <  1:1055324/137‑1 (MQ=38)
                                    cGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgg                                                                                        <  1:10928196/137‑1 (MQ=38)
                                      tgCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATGGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                      <  2:9416703/136‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCTCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                      <  1:3815729/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                      <  1:17022219/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                      <  1:10039462/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                      <  1:1605578/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                      <  1:12478710/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                      <  1:5950758/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                      <  2:13969733/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                      <  1:15104290/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                      <  1:633771/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                      <  1:8039380/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                      <  1:14575738/137‑1 (MQ=38)
                                            caggAGCAATTCCTTTTGTTATGGGTCGCCGACAGCATAGTAACACCACTCATGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATGTCCGCGTGTATGCAACAAACCTGGGGCTGCTa                                                                                <  2:1259931/134‑1 (MQ=25)
                                            cATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTa                                                                                <  2:8519142/137‑1 (MQ=38)
                                                        cTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCATC‑‑‑G‑‑‑C‑GGGCCC‑GGTGATTTTCGCGTGTGTGTAAAAAACCTGGGGGTACTCTGTAGCACtttt                                                                    >  2:15551272/1‑136 (MQ=19)
                                                             ttCATGTGTCGCCGACAGCATAAAAACACCGCTAGGGTGGGGAAGCGGCCTTGTTTCGCGAAAGGGCCGCATAGTGG‑‑‑C‑C‑‑‑CC‑GGCAATTTCCGTGTGTGTGCAACAAACCTGGGGCTGCTGTGCAGCCCTTTCGCGac                                                               <  2:14864456/136‑1 (MQ=25)
                                                             tcaaTGGGTCGCTGACAGCGGAGCACCACCCCTCAGGTGTGGGAGCGTCCTTGTGCCGCGAAGGGGCCGCAGA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAAAAAACCTGGGGCTGCTGTGCAGCCCTTTCGCGac                                                               <  2:11587056/134‑1 (MQ=19)
                                                               cATGTTTCGCCGACAGAATAGCAACACCATTCAGGTGTGGGAGCGGCCTTGTGTCGCAAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCACGTGTATGCAACAAACCTGGGGCGGCGATGCAGCCCTTTCGCGacac                                                             <  2:4868744/137‑1 (MQ=25)
                                                                 tgttTCGCCGACAGCACAGCAACACCTCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCGTA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACaa                                                           <  2:10781321/137‑1 (MQ=37)
                                                                 tgtgTCTCCGACGGCATAGAAACACCGCTCAGGTGTGGGAGCGGCCTTGTTTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGTGTGGATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACaa                                                           <  2:4376702/137‑1 (MQ=25)
                                                                 tgtgTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAGGGGCAGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACaa                                                           <  2:7729934/137‑1 (MQ=38)
                                                                 ggtgtCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAGA‑‑‑G‑‑‑C‑GGC‑CCTGGCAATTTCCGCGTGGATGCAAAAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACaa                                                           <  2:2233303/136‑1 (MQ=25)
                                                                  ttgtCGCCGACAGCATAGCAACACCGCTCGGGTGTGGGAGTGACCTTGTGTCGCGAAAGGTCCGCATG‑‑‑G‑‑‑C‑GGCCCC‑GACAGTTTCCGCGTGTATGCAACAAACTTGGTGTTGCTGTGCAGCCCTTTCGCGACACAAg                                                          <  2:4264238/136‑1 (MQ=255)
                                                                      cGCCGACAGCAAAACAACCCCCTTCAGGTGGGGGAGCGGCTTTGTGTCGGGAAAGGGCCGAAGA‑‑‑G‑‑‑G‑GGC‑CCAGGCAGTTTCCGCGGGTATGAAACAAACTTGGGGTTGCTATGCAGCCCTTTCGCGACACAAGGccc                                                      <  2:14280234/137‑1 (MQ=255)
                                                                                agaacAACACCGCTCAGGTGGGGGAGTGGCTTGGTGTAGAGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGAAATTTCCGTGTGTGGGAAACAAACCTGGGGCTGCTGTGCAGCCCTTTCGCGACACAAGGCCCACAAAGGAcc                                            <  2:13827636/133‑1 (MQ=19)
                                                                                acaGCACCACCTCTCAGGTGTGGGAGCGGCCTTGTTTCGCGAAAGGGCCGAATG‑‑‑GTGAC‑C‑‑‑CC‑GGAAGTTTCCGCGGGGATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGAcc                                            <  2:10919311/135‑1 (MQ=25)
                                                                                aTAGCAACCCCTCTCAGTTGTGGGAGTGGCTTTGTTTCGCGAAAGGGCTGAATG‑‑‑G‑‑‑T‑GGC‑CCTGGCACTTTCCGCTTGCATGCAACAACCCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGAcc                                            <  2:11816472/137‑1 (MQ=255)
                                                                                aTAGCAACACCTCTCAGGTGTGGAAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTTGCGACACAAGGCCCACAAAGGAcc                                            <  2:18190858/137‑1 (MQ=37)
                                                                                aTAGAAACACCGCTAGGGTGGGGGGGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑GTGTC‑C‑‑‑CC‑GGAAATTTCTGTGTGGAGGAAAAAAACTTGGGGTGGTCATGCATCCTTTTCGCAACACAAGGCCCACAAAGGAcc                                            <  2:10576481/137‑1 (MQ=255)
                                                                                  agaaaCACCTCTCAGGTGTGGGAGCGGCTTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGTGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACcgc                                          <  2:11454884/134‑1 (MQ=25)
                                                                                  aacAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATATCGG‑‑‑C‑C‑‑‑CC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACcgc                                          <  2:4436433/135‑1 (MQ=38)
                                                                                  aGCAACACCTCTCAGTTGGGGGAGCGGCCTTTTGTCGCGAAAGGGCCGCAGA‑‑‑G‑‑‑T‑GGCCCC‑GGCACTTTCCGCGTGTATGCAACAACCCTGGGGCTGTGATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACcgc                                          <  2:9844749/137‑1 (MQ=25)
                                                                                  aGCAACACCGTTCGGGTGTGGGGGTGGCTTTGTGTTGCGAAGGGGCCGCATG‑‑‑G‑‑‑T‑GGCCCC‑GGCAGTTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACcgc                                          <  2:3470507/137‑1 (MQ=255)
                                                                                  aGCAACACCGTCCAGGTGGGGGAGCGGCCTTGTTTCGCAAAGGGGCTTCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCACTTTCCGCGGGCAGAAAACAAACCTGGGGCTGCTATGCGGCCCTTTCGCGACACAAGGCCCACAAAGGACcgc                                          <  2:2933258/137‑1 (MQ=255)
                                                                                  aGCAACACCGCTCAGGTTTGGGAGCGGCCTTGTGTCGCAAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCGGCTATGCAACCCTTTCGCGACACAAGGCCCACAAAGGACcgc                                          <  2:7467381/137‑1 (MQ=37)
                                                                                  aGCAACACCGCTCAGGGGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAGTTTCCGCGTGGATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACcgc                                          <  2:17517639/137‑1 (MQ=37)
                                                                                       cacCGCTCAGGTGTGGGGGGGGCCTTGTGTCGCGAAGGGGCAGAATA‑‑‑G‑‑‑C‑GGCCCC‑GGCACTTTCCGTGTGTAGGCAACAAACCGGGGGTTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAgg                                     <  1:15762797/137‑1 (MQ=25)
                                                                                        acCGCTCAGGTGGGGGAGCGGCCTTGTGTCGCGAAAGGTCCGCGTA‑‑‑G‑‑‑C‑GGCCCC‑GGCAGTTTCCGTGTGTATGCAACAAACCTGGGGCTGCTATGCATCCTTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGc                                    <  2:6466056/137‑1 (MQ=25)
                                                                                             tCAGGTGTGGGAGCGGCCTTGTGTTGCGAAAGGGCCGCATA‑‑‑T‑‑‑C‑GTCCCC‑GGCAATTTCCGCGTGTATGCAAAAAACCTGGGGTTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAg                               <  2:16434041/137‑1 (MQ=37)
                                                                                             tCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAGTTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAg                               <  2:8380859/137‑1 (MQ=37)
                                                                                                  tgtgGGAGTGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACTTGGGGCTGCTATGCAGCCCTTTCTCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                          <  1:13486052/137‑1 (MQ=37)
                                                                                                  tgtgGGAGCGGCTTTGTGTCGCGAAAGGGCCGCATAGTGG‑‑‑C‑C‑‑‑CC‑GGCAGTTTCCGCGTGTAGGAAACAAACCTGGGGCGGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                          <  2:1888200/137‑1 (MQ=37)
                                                                                                  tgtgGGAGCGGCTTTGTGCCAGAAGAGGGCCGCAGAGAGG‑‑‑C‑C‑‑‑CC‑GGAAGTTTCACCGGGTATGCAAAAAACTTGGGGTTGCTATGAGGCCTTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                          <  2:12299493/137‑1 (MQ=255)
                                                                                                  tgggggATCGGCCTTGTGTCGCGAAAGGGCCGCGGA‑‑‑G‑‑‑C‑GGC‑CCAGGCAATTTCCGAGTGTGTGCAACAAACCTGGGGCTGCTGGGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                          <  2:12861525/134‑1 (MQ=25)
                                                                                                  ggggggAGCGGCCTTGGGTAGAGAAGGGCCCGCAGA‑‑‑G‑‑‑C‑GGC‑CCAGTCAGTTTCCGTGTGTATGCAACAAACTTNTGGTTGCTATGCAGCCCTTTCCCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                          <  1:17403120/134‑1 (MQ=255)
                                                                                                   ttgGGAGCGGCTCTTTGTCGCGAAAGGGCTGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTCCCGCGGGTATGCAACAAACTTGGGGCTGCTATGCACCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCgg                         <  2:15919102/136‑1 (MQ=255)
                                                                                                   gtgGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGC‑CCTGGCAATTTCCGTGTGTAGGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCgg                         <  2:10612428/137‑1 (MQ=37)
                                                                                                        aGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGAAACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTg                    <  2:2072253/137‑1 (MQ=38)
                                                                                                             ctttGTTTCGCGAAAGGGCCGCATA‑‑‑G‑‑‑CTGCC‑CC‑GGCAATTTCCGCGTGTATGCAACAAACCTTGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCa               <  2:18126712/135‑1 (MQ=255)
                                                                                                             ccTTTTGTCGCGAAGGGGCCGAATATT‑GT‑‑C‑C‑‑‑CT‑GGCATTTTCCGCTTGTATGCAACAACCTTGGGGCTTCGGTACATCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCa               <  1:18311082/137‑1 (MQ=255)
                                                                                                                  tgtagagAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGC‑CCTGGCATTTTCCGTGTGTAGGCAACAAACTTGGGTCTGCTATGCATCCCTTTCGCAACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTg          <  1:9978748/131‑1 (MQ=255)
                                                                                                                  tgtCGCGAGGGGGCCGCATA‑‑‑G‑‑‑C‑GGC‑CCTGGCACTTTCCGTGTGTGTAAAACAACCTTGGGTCTGCTATGCATCCTTTTCTCAACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTg          <  1:7531626/137‑1 (MQ=255)
                                                                                                                  tgtCGCGAAAGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGCATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTg          <  2:6824878/137‑1 (MQ=38)
                                                                                                                  tgtCGCGAAAGGGCCGCAGA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTAGGAAACAAACTTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTg          <  1:18761518/137‑1 (MQ=37)
                                                                                                                          aaGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCCAAAAGGAACGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAg  >  2:8020629/1‑137 (MQ=255)
                                                                                                                          aaGGGCCGCATA‑‑‑G‑‑‑C‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGATGCTCTGCAGCCCTTTTTCGAAACAAAGCCCACAAAGGACCGCGCAGGATTCCGGCCCGGAACTGGCCCAGGCTGGTGCGCAg  >  2:10960442/1‑137 (MQ=255)
                                                                                                                                                |||                                                                                                                        
TGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATA‑‑‑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 < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑