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.G94.5% 650.1 / 38.6 230intergenic (‑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/12);  new base G (8/200);  total (9/221)
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.C94.3% 756.5 / 40.0 230intergenic (‑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 . (1/12);  new base C (7/209);  total (9/221)
Fisher's exact test for biased strand distribution p-value = 3.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*AE015451_phaM15,555,6083.C87.8% 719.2 / 103.4 230intergenic (‑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 . (1/27);  new base C (8/194);  total (9/221)
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

Only 100 of 236 total aligned reads displayed.

TGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑G‑‑‑CC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCT  >  AE015451_phaM1/5555472‑5555729
                                                                                                                                                |||                                                                                                                          
              gTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTc                                                                                                                <  1:10261331/137‑1 (MQ=38)
                cGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCcg                                                                                                              <  1:6985356/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCgt                                                                                                           <  1:3636947/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCgt                                                                                                           <  1:248005/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCgt                                                                                                           <  2:2585074/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCgt                                                                                                           <  1:645655/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCgt                                                                                                           <  1:12776866/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCgt                                                                                                           <  1:12450745/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCgt                                                                                                           <  2:8427019/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCgt                                                                                                           <  1:11201663/137‑1 (MQ=38)
                        gCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATg                                                                                                      <  2:2655986/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TAGTGGC‑C‑‑‑CC‑GGCAATTTCCGCGTGTATGc                                                                                                     <  2:10163652/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TAGTGGC‑C‑‑‑CC‑GGCAATTTCCGCGTGTATGc                                                                                                     <  2:6259154/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                     <  2:11400432/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                     <  1:6262497/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                     <  2:14634694/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                     <  2:1558118/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                     <  1:9187689/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                     <  2:3574597/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                     <  2:4802917/137‑1 (MQ=38)
                         cAGCGATTAGGCGGGCGGACATGAGCACCTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                     <  2:14394423/137‑1 (MQ=37)
                            cGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGcaac                                                                                                  <  1:7700270/137‑1 (MQ=38)
                              aTTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGcaacaa                                                                                                <  1:14640707/137‑1 (MQ=38)
                              aTTAGTCGGGCGGACAGGAGCAACTCCTTTTGTCATGTGTCGCCGATAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGcaacaa                                                                                                <  1:4769595/137‑1 (MQ=37)
                               ttAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACaaa                                                                                               <  2:16501069/137‑1 (MQ=38)
                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCGTGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                             <  1:292897/137‑1 (MQ=38)
                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCTCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                             <  1:6820123/137‑1 (MQ=38)
                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                             <  2:12823130/137‑1 (MQ=38)
                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                             <  1:14670843/137‑1 (MQ=38)
                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                           <  2:4480582/137‑1 (MQ=38)
                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                           <  2:7873825/137‑1 (MQ=38)
                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                           <  1:10765142/137‑1 (MQ=38)
                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                           <  1:2370994/137‑1 (MQ=38)
                                    cGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgg                                                                                          <  2:5028335/137‑1 (MQ=38)
                                    cGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgg                                                                                          <  2:8578557/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                        <  2:1113326/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                        <  2:13684790/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                        <  1:1574077/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                        <  1:16108474/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                        <  1:13487623/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                        <  1:4547100/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                        <  1:16737681/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                        <  1:7811333/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                        <  1:16842173/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                        <  1:3742633/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                        <  1:9924712/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                        <  1:8539223/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                        <  1:3442406/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTGGGTCATGTGTCGCCGACAGCATAGCAACACCGCTCGGGTGTGGGGGCGGCTTTGTGTCGCGAAAGTGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTATATGCAACAACCCttttc                                                                                        <  2:10057073/137‑5 (MQ=255)
                                            cATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTa                                                                                  <  1:12112932/137‑1 (MQ=38)
                                              tGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATg                                                                                <  2:9212202/137‑1 (MQ=38)
                                                      tCCTTTTGTCTTTTGTCGCCGACAGCATAGCAACACCACTCGGGTGTGGGAGTGGCCTTGTTTCGCGAAAGGGCCACA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCACGTGTATGCAACAAACCTGGGGCGGCTATGCAGCCCtt                                                                        <  2:13163343/137‑1 (MQ=255)
                                                      tCCTTTTGTCATGTTTCGCCGACAGCATAGCACCACCGCTCATGTGTGGGAGCGGCCTTGTGTCGAGAAAGGTCCACA‑‑‑TA‑‑‑CC‑GACCCC‑GGCAATTTCCGCGTGTAGTCAACAAACCCGGGTTTCCTATGCAGCCCtt                                                                        <  2:1688712/137‑1 (MQ=255)
                                                      tCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGNGTGTATGCAACAAACCTGGGGCTGCTATGCAGcccca                                                                        <  1:16290945/137‑3 (MQ=38)
                                                      tCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCtt                                                                        <  2:7839405/137‑1 (MQ=38)
                                                      tCCTTTTGCCATGTGTCGCCGACAGCGTAACAACACCTCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACTTGGGGCTGTTATGCAGCCCtt                                                                        <  2:10314122/137‑1 (MQ=25)
                                                             ttCATGTGTCGCCGACAGCATAGCAACACCACCCAGTTGTGGGAGCGGCCTTGTTTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GT‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGTCTGCTATGCATCCCTTTCGCGac                                                                 <  2:11573290/136‑1 (MQ=255)
                                                             ttCATGTGTCGCCGACAGCATAGAAACACCTCCCAGGTGTGGGAGCGGCCTTGTGTCGAGAAAGGGCCGCA‑‑‑TA‑‑‑TC‑GGCCCC‑GACAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGac                                                                 <  2:3261764/136‑1 (MQ=37)
                                                             gTCATGTGTCGCTGACAAAATAGAAACACCGCTAAGGTGTGGGAGCGGCCTTGTGCCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GCTGCC‑CC‑GGCAATTTCCGCGTGTGTGAAACAAACCTGGGGCTGCGGTGCAGCCCTTCCGCGac                                                                 <  2:11964227/137‑1 (MQ=25)
                                                              tCATGTGTCACCGACAACATAGCAACCCCTCTCAGGTGTGGGAGCGGCCTTGTGTCGCAAAAGGGCCGCA‑‑‑TAGTGGC‑C‑‑‑CC‑GGCAATTTCCGTGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCTTTTCGCGaca                                                                <  2:2922573/137‑1 (MQ=25)
                                                               cgtgtgTCGCCGACAGCAGAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTTTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGacac                                                               <  2:16237769/135‑1 (MQ=38)
                                                               cATTTGTCGCCGACAGCATAGCAACACCGTTCAGGTTTGGGGGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTTGATGCAAAAAACCTGGGGCTGCTATGCAGCCCTTTCGCGacac                                                               <  2:9545257/137‑1 (MQ=25)
                                                               cATGTTTCGCCAAAAGCAGAGCACCACCGCCCAGGTGTGGGAGCGGCTTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGacac                                                               <  2:5578877/137‑1 (MQ=25)
                                                                 tttttCGCCGACAGCGTAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGAA‑‑‑TA‑‑‑GCTGCC‑CC‑GGAAATTTCCGCGTGTAGGAAACAAACCTGGGTCTGCTATGCAGCCCTTTCGCGACACaa                                                             <  2:4296894/133‑1 (MQ=25)
                                                                 tttgtCGCTGAAAGCGTGGCAACCCCTCTCGGGTGTGGGGGGGGCTTTGTTTTGTGAAGGGGCCGCG‑‑‑TG‑‑‑GT‑GGC‑CCTGGCAGTTTCTGTGTGGGTGAAAAAAACTTGGGGTTGTTGTGCAGCCCTTTCGCGACACaa                                                             <  2:4996810/135‑1 (MQ=255)
                                                                 tgttttgCCGACAGCATAGCAACACCGTTCAGGTGTGGGAGCGGCCTTTTGTCGCGAAAGGGCCACA‑‑‑TA‑‑‑TC‑GGCCCC‑GGAAATTTCCGCGTGCATGCAACAAACCTTGGGCTGCTGGGCAGCCCTTTCGCGACACaa                                                             <  2:16587139/131‑1 (MQ=25)
                                                                 tgtgTCGCCGACAGCATAGCAACACCTCTCAGTTGTGGGAGCGGCTTTGTGTCGCGAAAGGGCCGCA‑‑‑TAGTGGC‑C‑‑‑CC‑GGCAATTTCCGCGTGTATGCAAAAAACCTGGGGCTGCTATGCAGCCTTTTCGCGACACaa                                                             <  2:2095948/137‑1 (MQ=25)
                                                                 tgtgTCACCGACAGCGTGGCAACACCTCCCGGGTGTGGGAGTGGCCTTGTGTCGCGAAGGGGCCCCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGAAATTTCCGTGTGTAGGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACaa                                                             <  2:8648697/137‑1 (MQ=21)
                                                                      cGCCGACAGCATAGCAACCCCGCTCAGGTGTGGGAGCGGCTTGGTGTCGAGAAGGGGCCGCAGGGTG‑‑‑GC‑C‑‑‑CC‑GGAAGTTTCCGCGTGGAGGAAAAAAACCTGGGGTGGTGATGAAGCCTTTTCGCGACACAAGGccc                                                        <  2:10875886/137‑1 (MQ=255)
                                                                      cGCCGACAGCATAACAACCCCCCTCAGGTGTGGGAGCGGCCTTGTGTCGCAAAAGGTCCACA‑‑‑TA‑‑‑GC‑GTC‑CCAGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCGGCGATGCAGCCCTTTCGCGACACAAGGccc                                                        <  2:2635192/137‑1 (MQ=37)
                                                                      cGCCGACAGCAGAGCAACACTGCTCAGGTGGGGGAGGGGCCTTGTGTAGAGAAGGGTCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCATTTTCCGTGTGTATGAAACAACCCTGGGGCTGCCGTGCAGCCCTTTCGCGACACAAGGccc                                                        <  2:14909029/137‑1 (MQ=255)
                                                                           acaacaGAGCAACACCGCCAAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCACCAACCCTGGGGCGGCTATGCAGCCCTTTCGCGACACAAGGCCCACaaa                                                   <  2:3890887/137‑1 (MQ=25)
                                                                           aCAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCGTGTGTCTCGAAAGGGCCGCA‑‑‑TG‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGTGATGCAGCCCTTTCGCGACACAAGGCCCACaaa                                                   <  2:14715558/137‑1 (MQ=37)
                                                                                agaacAACACTGCTCGGGTGTGGGAGGGGCTTGGTGTCGAGAAGGGGCCGCA‑‑‑GAGGGGC‑C‑‑‑CC‑GGAAGTTTCCGCGGGGGTGAAACAAACCTGGGGCGGCGAGGAAGCCCTTTCGAGACACAAGGCCCACAAAGGAcc                                              <  2:186842/133‑1 (MQ=255)
                                                                                aTAGCAACACCGCTCAGGTGTGGGAGCGGCTTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGC‑CCTGGAAATTTCCGCGTGGATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGAAACAAGGCCCACAAAGGAcc                                              <  2:10482564/137‑1 (MQ=37)
                                                                                aTAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTTTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCTGCGTGTATGCAACAAACCTGGGGTTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGAcc                                              <  2:16339241/137‑1 (MQ=38)
                                                                                  aacAACACCGTTCGGGTGGGGGAGAGGCTTGGTGTTGTGAAGGGGCGGAA‑‑‑TA‑‑‑GT‑GGC‑CCTGGAAGTTTCTGCGGGTATGAAACAAACCTGGGGTTGCTGGGCAGCCTTTTCGAGAAACAAGGCCCACAAAGGACcgc                                            <  2:894501/135‑1 (MQ=255)
                                                                                  aGCAACACCGCTCAGGGTTGGGAGCGGCCTTGTTTCTCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGAGTGTATGCAGCAACCTTGGGTCGGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACcgc                                            <  2:14661377/137‑1 (MQ=25)
                                                                                        acCGCTCAGGTGTGGGATCGGCCTTGTGTCGCGAAAGGGACACA‑‑‑TAGAGGC‑C‑‑‑CC‑GGCAGTTTCCGTGTGTAGGCAACAAACCTGTGTTTGCGATGCAGCCCTTTCGAGAAACAAGGCCCACAAAGGACCGCGCAgga                                      <  2:15178972/137‑2 (MQ=19)
                                                                                        acCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGAA‑‑‑TA‑‑‑GC‑GGC‑CCTGGCAATTTCCGCGTGTGTGCAACAAACCTGGGGCTGCGATGCAGCCTTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGc                                      <  2:12502573/137‑1 (MQ=37)
                                                                                        acCGCTCAGGTGTGGGAGCGGCCTTGTGTCCCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GCCCCC‑GGCAATTTCCGCGTGGATACAACAAACCTGGGGCTGCTATGCAGCCCTTTCTCGACACAAGGCCCACAAAGGACCGCGCAGGc                                      <  2:12301097/137‑1 (MQ=37)
                                                                                             tcgggTGGGGGGGCGGCCTTGTGTCGAGAAAGGGCCGCG‑‑‑GG‑‑‑GG‑GGC‑CCTGGCACTTTCTGTGTGTATGAAACAAACCTGGGGTTGCTAGGCAGCCTTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAg                                 <  2:251889/134‑1 (MQ=255)
                                                                                             tCAGGTGTGGGAGCGGCCTTTTGTCGCGAAAGGTCCGCA‑‑‑TATCGGC‑C‑‑‑CC‑GTCAATTTCCGCGTGTATGCAACAAGCCTGGTGCTGCTATGCAGCCTTTTCGTGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAg                                 <  2:15723129/137‑1 (MQ=255)
                                                                                             tCAGGTGTGGGAGCGGCCTTGTGTCTCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCGACTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAg                                 <  2:7014349/137‑1 (MQ=37)
                                                                                             tCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GACAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAg                                 <  2:16317385/137‑1 (MQ=38)
                                                                                             ccAGGTGTGGGAGCGGCTTTGTTTCGAGAAAGGGCCACA‑‑‑TATCGGC‑C‑‑‑CC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCGATGCAGCCTTTTCTCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAg                                 <  2:15422532/136‑1 (MQ=255)
                                                                                                  tgtgGGAGCGGCTTTTTGTCGCGAAGGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAGTTTCCGTGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                            <  2:5438768/137‑1 (MQ=37)
                                                                                                  tgtgGGAGCGGCCTTGTTTCGCGAAAGGTCCGCA‑‑‑CA‑‑‑GC‑GCCCCC‑GGCAATTCCCGTGTGGATGCAACAAACCTGGGGCTGCTGTGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                            <  2:4832861/137‑1 (MQ=25)
                                                                                                  tgtgGGAGCGGCCTTGTTTCGCGAAAGGGCTGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                            <  2:8285315/137‑1 (MQ=38)
                                                                                                  tgtgGGAGCGGCCTTGTTTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGTGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                            <  2:10720031/137‑1 (MQ=37)
                                                                                                  tgtgGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCTGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCTCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                            <  2:4078835/137‑1 (MQ=37)
                                                                                                  tgggggAGTGGCTTGGTGTCGAGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGAAGTTTCCGTGTGTATGAAAAAAACCTGGGGCTGCTAGGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                            <  1:5869540/134‑1 (MQ=25)
                                                                                                   gtgGGAGTGGCCTTGGGTCGCGAAAGGGCCGCA‑‑‑GA‑‑‑GC‑GGCCCC‑GGCAGTTTCCGCGTGTAGGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCgg                           <  2:7614104/137‑1 (MQ=37)
                                                                                                   gtgGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGGAGACAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCgg                           <  2:3318315/137‑1 (MQ=25)
                                                                                                   gggggAGCGGCCTGGTGTCGAGAAAGGGCCGAA‑‑‑GA‑‑‑GC‑GGCCCC‑GGCAATTTCCGGGGGTGTGCAAAAAACCTGGGGTTGCGATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCgg                           <  2:933285/135‑1 (MQ=21)
                                                                                                        aGCGGCTTGGTGTCGAGAAGGGGCTGCG‑‑‑TG‑‑‑GC‑GGC‑CCTGGCAGTTTCGGTGTGTGTGCAACAAACTTGGGTCTGCTATGCAACCTTTTCGAGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTg                      <  1:4894323/137‑1 (MQ=255)
                                                                                                                          aaGGGCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAAGTTCCGCGTGTATGCAAAAAACCTGGGGCTGCTATGTAGCCCTTTCGCGAAACAAAGCCCACAAAGGAACGCGCAGGCTTCAGGCCCGGAACTGGGCCAGGCTGGCGCGCAg    >  2:7896785/1‑137 (MQ=255)
                                                                                                                          aaGGGCCGCA‑‑‑TA‑‑‑GC‑GGACCC‑GGCAAGTTCCGCGTGTCTGCAACAAACCTGGGGTTGCTATGCAGCCCCTTCGCGACACAAAGCCCAAAAAGGACCGCGCAGGATTCAGGCCCGGAACTGGCCCAGGGTGGAGCGTAg    >  2:13558190/1‑137 (MQ=255)
                                                                                                                          aaGGGCCACA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGTGTGTATGCAAAAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAg    <  2:3094175/137‑1 (MQ=38)
                                                                                                                            gggCCGCA‑‑‑TA‑‑‑GC‑GGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCt  >  2:10231799/1‑137 (MQ=255)
                                                                                                                                                |||                                                                                                                          
TGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑GC‑G‑‑‑CC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCT  >  AE015451_phaM1/5555472‑5555729

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: ‑