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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM15,555,6081.G93.4% 702.1 / 52.5 232intergenic (‑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/15);  new base G (15/196);  total (16/216)
Fisher's exact test for biased strand distribution p-value = 6.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
*AE015451_phaM15,555,6082.C93.5% 772.0 / 48.8 232intergenic (‑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/15);  new base C (14/201);  total (16/216)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*AE015451_phaM15,555,6083.C90.1% 759.4 / 83.9 232intergenic (‑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/23);  new base C (16/193);  total (16/216)
Fisher's exact test for biased strand distribution p-value = 3.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 243 total aligned reads displayed.

TGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CG‑‑‑CC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCT  >  AE015451_phaM1/5555472‑5555729
                                                                                                                                                  |||                                                                                                                          
tgCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑g‑‑‑cg                                                                                                                               <  1:17444247/137‑1 (MQ=255)
tgCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑g‑‑‑cg                                                                                                                               <  2:10144163/137‑1 (MQ=255)
      aaaTTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑g‑‑‑cggccccg                                                                                                                         <  2:14798998/137‑7 (MQ=255)
              gTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACTGCTCAGTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTc                                                                                                                <  2:18510120/137‑1 (MQ=37)
                  gCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCcgcg                                                                                                            <  2:5650493/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCgt                                                                                                           <  2:9810082/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCgt                                                                                                           <  1:6151847/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCgt                                                                                                           <  1:7175255/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCgt                                                                                                           <  2:14934065/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCgt                                                                                                           <  2:16046520/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCgt                                                                                                           <  2:2558486/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCgt                                                                                                           <  2:4463037/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCgt                                                                                                           <  2:7016591/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCgt                                                                                                           <  2:8341350/137‑1 (MQ=38)
                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCgt                                                                                                           <  2:8988308/137‑1 (MQ=38)
                        gCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATg                                                                                                      <  1:12619504/137‑1 (MQ=38)
                        gCAGCGATGAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATg                                                                                                      <  1:2095295/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                     <  2:2222295/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                     <  2:1220617/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                     <  2:9341651/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                     <  2:1220538/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                     <  1:1597670/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                     <  1:13277455/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                     <  1:17128335/137‑1 (MQ=38)
                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                     <  1:14959507/137‑1 (MQ=38)
                          aGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGca                                                                                                    <  2:4001396/137‑1 (MQ=38)
                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                             <  1:8274158/137‑1 (MQ=38)
                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                             <  1:12098913/137‑1 (MQ=38)
                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                             <  2:10992475/137‑1 (MQ=38)
                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                             <  1:5132624/137‑1 (MQ=38)
                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                             <  2:13240748/137‑1 (MQ=38)
                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                             <  1:2213074/137‑1 (MQ=38)
                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                             <  2:18205557/137‑1 (MQ=38)
                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                             <  1:13914113/137‑1 (MQ=38)
                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                           <  2:17802781/137‑1 (MQ=38)
                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                           <  1:15449774/137‑1 (MQ=38)
                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                           <  1:14108723/137‑1 (MQ=38)
                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                           <  1:17350428/137‑1 (MQ=38)
                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                           <  1:770092/137‑1 (MQ=38)
                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                           <  1:6807256/137‑1 (MQ=38)
                                    cGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgg                                                                                          <  1:14286930/137‑1 (MQ=38)
                                    cGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgg                                                                                          <  1:10242723/137‑1 (MQ=38)
                                     gggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTggg                                                                                         <  2:15133697/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTtggg                                                                                        <  2:7912130/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                        <  2:14266961/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                        <  1:11222225/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                        <  2:8493595/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                        <  1:2960239/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                        <  1:13136081/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                        <  1:4618868/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                        <  2:14891476/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                        <  1:15323056/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                        <  1:8170131/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                        <  1:5042615/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                        <  1:8869143/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                        <  1:925344/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                        <  2:12050163/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                        <  1:16579188/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                        <  2:10665600/137‑1 (MQ=38)
                                      ggCGGACATGAGCAACTCCCTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                        <  2:1636941/137‑1 (MQ=38)
                                           aCATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGgctgct                                                                                   <  1:5963711/137‑1 (MQ=38)
                                           aCATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGgctgct                                                                                   <  1:13092374/137‑1 (MQ=38)
                                               gAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGc                                                                               <  2:5403227/137‑1 (MQ=38)
                                               gAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGc                                                                               <  1:9744382/137‑1 (MQ=38)
                                                      tCCTTTTGTCCTGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCTTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGGGTATGCAACAAACCTGGGGCTGCTATGCAGCCCtt                                                                        <  2:10480872/137‑1 (MQ=37)
                                                             gTCATGTGTCGCCGACAGCATAACAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTAGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGac                                                                 <  2:13745681/137‑1 (MQ=38)
                                                                 tttttCGCCGACAGCATAGCAACCCCGCTCGGTTTTGGGAGCGGCCTTGTGTCGCGAAGGGGCCGCA‑‑‑TAGTGG‑‑‑CC‑‑‑CC‑GGCAGTTTCCGCGTGTATGCAACAAACCTGGGGTTGCTGGGCAGCCCTTTCGCGACACaa                                                             <  2:7744202/133‑1 (MQ=25)
                                                                 tttgtCGCCGACAGCAGAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGGGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATAAAAAAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACaa                                                             <  1:18133880/135‑1 (MQ=37)
                                                                 tgtgTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCTTTGTGTCGCGAAAGGGCCGCA‑‑‑TG‑‑‑G‑‑‑CGTCCCC‑GGCAATTTCTGCGTGGATGAAAAAAACTTGGGGTTGCGATGCAGCCCTTTCGCGACACaa                                                             <  1:4356144/137‑1 (MQ=25)
                                                                 tgtgTCGCCGAAAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCACA‑‑‑TA‑‑‑G‑‑‑CGACCCC‑GACAATTCCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACaa                                                             <  2:7525298/137‑1 (MQ=38)
                                                                 tgtgTCGCCAACAGCATAGCAACACCTCCCAGTTGTGGGAGCGGCCTTGTGTCGCGAAGGGGCCGCG‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAACCCTGGGTCTGCTGTGCAGCCCTTTCGCGACACaa                                                             <  2:15358820/137‑1 (MQ=25)
                                                                 tgtgTCGCCAACAGCATAGCAACACCCCTCAGGTGTGGGAGCGGCCTTGTTTCACGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCACTTCCCGCGGGTATACAACAAACCTGGGGCTACTATGCAGCCCTTTCGAGACACaa                                                             <  2:1781412/137‑1 (MQ=255)
                                                                 tgtgCCGCCGACAACACAGCAACACCGCTCAGGTGGGGGAGCGGCTTTGTTTCGAGAAGGGGCCGCATGGTG‑‑‑G‑‑‑CC‑‑‑CC‑GGCAATTTCCGCGTGCAGGAAACAAACTTGGGGTTACTATGCAGCCTTTTCGCGAAACaa                                                             <  2:444414/137‑1 (MQ=255)
                                                                     tCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAACTTCCGCATGTACGCAAAAAACCCGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGcc                                                         >  1:1761064/1‑137 (MQ=25)
                                                                      cGCCGAAAGCAGAAAAACACTGCTGGGGGGGGGGGGTGGCTTTGTGTAGAGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGAAGTTTCTGTGGGGGGGAAACAAACCGGGGGTTGCGATGCATCCCTTTCGCGACACAAGGccc                                                        <  2:15883582/137‑1 (MQ=255)
                                                                           aCAGCATAGCAACACCTCGAAGGTTTGGGAGCGGCCTTGTTTCGCGAAAGGGCCACA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATGTCCGCGTGTATGCAACAAACCTGTGGCTGCTATGCAGCCCTTTCGAGACACAAGGCCCACaaa                                                   <  2:17682082/137‑1 (MQ=255)
                                                                                acaGCAACACCGCTCAGGTGTGGGATCGGCCTTGTGTCGCAAAAGGGCCGAA‑‑‑GA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGGATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGAcc                                              <  1:15650711/135‑1 (MQ=25)
                                                                                  aacAACACCGCTAAGGGGGGGGAGCGGCCGTGGGTCGAGAAGGGGCCGAA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGAAATTCCTGCGTGTATAAAAAAACCCTGGGCCTGCTATGCACCCCTTCCGAAACACAAGGCCCACAAAGGACcgc                                            <  1:12599707/135‑1 (MQ=255)
                                                                                  aGCAACACCGCTCAGGTGTGGGAGAGGCCTTGTGTCGCAAAAGGGACACA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCACACGTATGCAACAAAACTGGGGCTGCTATGCAACACTTTCGCGACACAAGGCACACAAAGGACcga                                            >  2:18378216/1‑136 (MQ=255)
                                                                                       ccccGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAGTTTCCGCGTGTGTGCAACAACCTTGGGGTTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAgg                                       <  2:3240678/135‑1 (MQ=25)
                                                                                        acCGCTCAGGTGGGGGAGTGGCCTTGTTTCGAAAAGGGGCCGCA‑‑‑GA‑‑‑G‑‑‑GGGC‑CCAGGCAGTTTCCGTGTGTAGGAAACAAACCTGGGGCTGCTATGAAGCCTTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGc                                      <  2:1052195/137‑1 (MQ=255)
                                                                                             tCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAAAAAACCTGGGGCGGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAg                                 <  2:15495377/137‑1 (MQ=37)
                                                                                             tCAGGTGTGGGAGCGACCTTGTGTCACAAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGAAATTTCCGCGTGTATGCAACAACCCTGGGGCTGCGAGGCATCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAg                                 <  2:14682564/137‑1 (MQ=255)
                                                                                               aGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAAAAAACCTGGGGCTGCTATGCAGCCCCTTCGAGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGc                               >  1:17917086/1‑137 (MQ=37)
                                                                                                  tgtgGGGGCGGCTTTGTGTCGCGAAAGGGCCGCA‑‑‑GA‑‑‑G‑‑‑CGGC‑CCTGGCAATTTCCGCGTGTATGCAACAAACCTGGGACTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                            <  2:12696649/137‑1 (MQ=37)
                                                                                                  tgtgGGAGTGGCCTTGTGTCGCAAAAGGGCCACA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTCCCGCGTGTATGAAACAAACCGGGGGCTGCTAGGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                            <  2:10075374/137‑1 (MQ=37)
                                                                                                   ttgGGAGCGGCCTTGTTTCGCGAAAGGGCCGCA‑‑‑GA‑‑‑G‑‑‑CGGC‑CCTGGCAATTTCCGCGGGGATGCAACAAACCTGGGGCTGCTAGACAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCgg                           <  2:1964872/136‑1 (MQ=25)
                                                                                                   gtgGGAGCGGCCTTGTTTCGCGAAAGGTCCGCG‑‑‑TG‑‑‑GTGTCC‑‑‑CC‑GGCACTTTCCGTGTGGCTGCAACAACCCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCgg                           <  1:17696115/137‑1 (MQ=255)
                                                                                                         gCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                             <  2:28261/65‑1 (MQ=17)
                                                                                                          cGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTTTATGCAACAACCCTGGTGCTGCTATGCATCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGc                    <  2:6725725/137‑1 (MQ=37)
                                                                                                                  tttCGAGAAGGGGCCGCA‑‑‑GA‑‑‑G‑‑‑CGCCCCC‑GGCACTTTCCGTGTGTGTAAAACAACCCTGGGTTTGCTGTACATCCTTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTg            <  1:1090449/135‑1 (MQ=255)
                                                                                                                  tgtCGAGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAGTTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTg            <  2:8561502/137‑1 (MQ=255)
                                                                                                                   ttCGCGAAAGGGCCGCA‑‑‑TG‑‑‑G‑‑‑CGGC‑CCTGGCAGTTTCCGTGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCTTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTgg           <  2:805771/136‑1 (MQ=255)
                                                                                                                         aaaGGTCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCACTTTCCGTGTGTATGCAACAAACTTGGGGCTGTTATGCATCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCa     <  1:10569640/137‑1 (MQ=255)
                                                                                                                          aaGGGCCGCA‑‑‑TC‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGATTCAGGCCCGGAACTGGCCCAGGGTGGCGTGCAg    >  2:18094729/1‑137 (MQ=255)
                                                                                                                          aaGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGTCCC‑GGCAATTTCCTCATGTATGCAACAAACCTGGGGCTGCTCTGCAGCCCTTTCGCGACACAAGGCCCACAAAGGAACGCACAGGCTTCCGGACCGGAACTGGCCCAGGCTGGCGCGCAg    >  2:4642990/1‑137 (MQ=255)
                                                                                                                          aaGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGATATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAAATGGCCCAGGCTGGAGCGCAg    >  2:3112967/1‑137 (MQ=255)
                                                                                                                          aaGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAATTTCCGCGTGTATGCAAAAAACCTGGGGCTGCTCTGCAGCCCTTTCGCGACACAAGGCCCAAAAAGGACCGCGCAGGATTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAg    >  2:8789610/1‑137 (MQ=255)
                                                                                                                          aaGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑GGCAAGTTCCGCATGTATGCAACAAAACTGGGGCTGCTATGCAGCCCTTTTGCGACACAAGGCCCACAAAAGACAGCGTAGGCTTCAGGCCCAGAAATGGACCAAGCTGGCGCGTAg    >  2:3679956/1‑137 (MQ=255)
                                                                                                                          aaGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CGGCCCC‑CGCAATTTCCGCGTGTATGCAACAAAACTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCCGGCTTCCGGCCCGGAACTGGCCCCGGCTTGCGCGCAg    >  2:12957103/1‑137 (MQ=37)
                                                                                                                                                  |||                                                                                                                          
TGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCA‑‑‑TA‑‑‑G‑‑‑CG‑‑‑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 < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑