Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 276,899 G→A G62E (GGG→GAG)  mmuM → CP4‑6 prophage; S‑methylmethionine:homocysteine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913276,8990GA96.4% 1366.3 / 24.4 412G62E (GGG→GAG) mmuMCP4‑6 prophage; S‑methylmethionine:homocysteine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (12/3);  new base A (201/196);  total (213/199)
Fisher's exact test for biased strand distribution p-value = 3.35e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 427 total aligned reads displayed.

AAACAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAACCC  >  NC_000913/276751‑277046
                                                                                                                                                    |                                                                                                                                                   
aaaCAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGGCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCggagg                                                                                                                                                   <  4:611954/151‑1 (MQ=255)
 aaCAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGgc                                                                                                                                                  <  5:492155/151‑1 (MQ=255)
  aCAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAATGCATCTTGATTACTACCGGGCGGAGgcg                                                                                                                                                 <  1:1007846/151‑1 (MQ=255)
        tatCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATgc                                                                                                                                           >  7:677449/1‑151 (MQ=255)
         atCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATgcg                                                                                                                                          >  3:1378573/1‑151 (MQ=255)
         atCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATgc                                                                                                                                           <  5:529386/150‑1 (MQ=255)
         atCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATgc                                                                                                                                           <  4:1124072/150‑1 (MQ=255)
            ctgctgctgGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGcgcg                                                                                                                                        >  7:546893/1‑150 (MQ=255)
            ctgctgctgGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGcgcg                                                                                                                                        >  3:649796/1‑150 (MQ=255)
             tgctgctgGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGAt                                                                                                                                      >  6:680327/1‑151 (MQ=255)
              gctgctgGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATc                                                                                                                                     <  6:419614/151‑1 (MQ=255)
                 gctgGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCAc                                                                                                                                   <  2:1009856/150‑1 (MQ=255)
                  ctgGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTg                                                                                                                                 >  2:1114/1‑151 (MQ=255)
                   tgGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGc                                                                                                                                >  3:1282695/1‑151 (MQ=255)
                        ggCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCt                                                                                                                           <  1:846737/151‑1 (MQ=255)
                        ggCGCGATGGCGAAGGAGCTGGAAGCGCGAGGGTTTAACTTAGCCGACAGCCTGTGGTCAGCCAAACTGATGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCt                                                                                                                           <  4:262811/151‑1 (MQ=255)
                         gcgcGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTa                                                                                                                          <  8:1555142/151‑1 (MQ=255)
                            cgATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATc                                                                                                                        >  3:397339/1‑150 (MQ=255)
                             gATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATc                                                                                                                        >  7:1519834/1‑149 (MQ=255)
                             gATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATc                                                                                                                        >  6:1431593/1‑149 (MQ=255)
                             gATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAg                                                                                                                      <  7:74256/151‑1 (MQ=255)
                             gATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAg                                                                                                                      <  4:1495141/151‑1 (MQ=255)
                             gATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCAGGGGCGCAATGCGCGATCACTGCCAGCTATc                                                                                                                        >  5:709083/1‑149 (MQ=255)
                              aTGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAgg                                                                                                                     <  4:767906/151‑1 (MQ=255)
                              aTGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAgg                                                                                                                     <  4:904767/151‑1 (MQ=255)
                                ggCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGc                                                                                                                    <  6:833832/150‑1 (MQ=255)
                                  cGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGAc                                                                                                                 >  6:1593977/1‑151 (MQ=255)
                                       gAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCGAGCTATCAGGCGACGCcgg                                                                                                            >  6:780155/1‑151 (MQ=255)
                                         gCTGGAAGCGCGAGTGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCAGCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCcggcg                                                                                                          <  1:463482/151‑1 (MQ=255)
                                            ggAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCAGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCTGGCGGAGGCGCAATGCCCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGc                                                                                                       >  5:186373/1‑151 (MQ=255)
                                              aaGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCtt                                                                                                     >  4:102140/1‑151 (MQ=255)
                                               aGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTc                                                                                                    <  4:938951/151‑1 (MQ=255)
                                                gcgcGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTcg                                                                                                   <  4:1313879/151‑1 (MQ=255)
                                                gagcgAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTcg                                                                                                   <  1:1137803/149‑1 (MQ=255)
                                                   cgAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTcgccg                                                                                                >  3:148945/1‑151 (MQ=255)
                                                   cgAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTcgccg                                                                                                <  8:1308980/151‑1 (MQ=255)
                                                   cgAGGGTGTAACTTAGCCGACAGCCTGTGGGCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTcgccg                                                                                                <  2:401346/151‑1 (MQ=255)
                                                    gAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTcgccgc                                                                                               >  5:1067913/1‑151 (MQ=255)
                                                    gAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTcgccgc                                                                                               >  7:245840/1‑151 (MQ=255)
                                                    gAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTcgccgc                                                                                               >  6:349131/1‑151 (MQ=255)
                                                     aGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCcgcg                                                                                              >  8:1191494/1‑151 (MQ=255)
                                                       ggTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCcgcgcg                                                                                            <  2:217812/151‑1 (MQ=255)
                                                          gtAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAgg                                                                                         <  2:735291/151‑1 (MQ=255)
                                                            aaCTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGt                                                                                        <  6:1071799/150‑1 (MQ=255)
                                                            aaCTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGt                                                                                        <  5:1392262/150‑1 (MQ=255)
                                                               ttAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTgg                                                                                    >  3:967221/1‑151 (MQ=255)
                                                               ttAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTgg                                                                                    >  3:772106/1‑151 (MQ=255)
                                                                 aGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGAt                                                                                  >  7:1410848/1‑151 (MQ=255)
                                                                  gCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATg                                                                                 <  7:1607320/151‑1 (MQ=255)
                                                                  gCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATg                                                                                 <  1:1281948/151‑1 (MQ=255)
                                                                   ccGATAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATg                                                                                 >  5:853371/1‑150 (MQ=255)
                                                                   ccGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATg                                                                                 >  4:515785/1‑150 (MQ=255)
                                                                   ccGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATg                                                                                 >  8:447665/1‑150 (MQ=255)
                                                                    cGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGaa                                                                               <  6:153780/151‑1 (MQ=255)
                                                                     gACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGa                                                                                >  5:1210351/1‑149 (MQ=255)
                                                                      aCAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAgc                                                                             <  8:84360/151‑1 (MQ=255)
                                                                       cAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAgcg                                                                            >  7:1072255/1‑151 (MQ=255)
                                                                        aGCCTGTGGTCAGCAAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAACTGCATCTTGATTAATACCGGGCTGAGGCGCAATGCGCGATCACTGCAAGCTATCAGGCGACGACGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAgcgc                                                                           >  1:197554/1‑151 (MQ=255)
                                                                               ggTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATGTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGa                                                                     >  2:1478148/1‑150 (MQ=255)
                                                                               ggTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCTATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGa                                                                     <  8:688785/150‑1 (MQ=255)
                                                                               ggTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGa                                                                     <  3:1135505/150‑1 (MQ=255)
                                                                                gTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGaa                                                                    >  5:1486719/1‑150 (MQ=255)
                                                                                 tCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCCCGAGGTCTGGATGAAGCGCAGGCGaa                                                                    >  6:952722/1‑149 (MQ=255)
                                                                                  cAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAgc                                                                 <  7:438909/151‑1 (MQ=255)
                                                                                  cAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAgc                                                                 >  6:1355496/1‑151 (MQ=255)
                                                                                   aGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAgcg                                                                <  1:1529816/151‑1 (MQ=255)
                                                                                     ccAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCt                                                              <  7:353946/151‑1 (MQ=255)
                                                                                     ccAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCt                                                              <  2:489209/151‑1 (MQ=255)
                                                                                      cAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTg                                                             <  8:194919/151‑1 (MQ=255)
                                                                                      cAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTg                                                             <  7:1318201/151‑1 (MQ=255)
                                                                                      cAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTg                                                             <  3:1613095/151‑1 (MQ=255)
                                                                                        aaGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGAt                                                           >  6:2743/1‑151 (MQ=255)
                                                                                                 tAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCaaa                                                    >  1:1437647/1‑149 (MQ=255)
                                                                                                   gAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCg                                                <  1:950226/151‑1 (MQ=255)
                                                                                                   gAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCg                                                <  2:1645894/151‑1 (MQ=255)
                                                                                                      aaCCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTgg                                             <  2:1181331/151‑1 (MQ=255)
                                                                                                      aaCCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTgg                                             >  2:1373714/1‑151 (MQ=255)
                                                                                                       aCCCGGAGCTTATCCGCGAGGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTggc                                            >  4:353366/1‑150 (MQ=255)
                                                                                                        cccGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTgg                                             >  8:1122307/1‑149 (MQ=255)
                                                                                                         ccGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGc                                          >  8:1602701/1‑151 (MQ=255)
                                                                                                          cGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCt                                         >  3:579982/1‑151 (MQ=255)
                                                                                                             aGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGgc                                      >  6:516666/1‑151 (MQ=255)
                                                                                                              gCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGgcg                                     >  4:1598673/1‑151 (MQ=255)
                                                                                                              gCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGgcg                                     <  7:965510/151‑1 (MQ=255)
                                                                                                                 tATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGgcgc                                    <  8:1043823/149‑1 (MQ=255)
                                                                                                                    ccGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTaaa                               <  5:343713/151‑1 (MQ=255)
                                                                                                                    ccGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTaaa                               <  5:64203/151‑1 (MQ=255)
                                                                                                                       cgAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGcc                            >  3:278883/1‑151 (MQ=255)
                                                                                                                        gAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGccc                           >  2:1304483/1‑151 (MQ=255)
                                                                                                                            tGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTg                        >  7:164915/1‑150 (MQ=255)
                                                                                                                                   tGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCg                   <  2:890887/148‑1 (MQ=255)
                                                                                                                                     aTTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCt              >  3:107940/1‑151 (MQ=255)
                                                                                                                                     aTTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCt              <  7:1212202/151‑1 (MQ=255)
                                                                                                                                     aTTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCt              >  3:928917/1‑151 (MQ=255)
                                                                                                                                     aTTACTACCGGGCGGAGGCGCAATGCGCGATCACCGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCt              <  4:1176345/151‑1 (MQ=255)
                                                                                                                                       tactacCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTgg            >  1:1652504/1‑151 (MQ=255)
                                                                                                                                       tactacCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTgg            >  5:1601729/1‑151 (MQ=255)
                                                                                                                                              gggcggaGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAgaga      >  8:312752/1‑150 (MQ=255)
                                                                                                                                              gggcggaGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGaa     <  6:1195617/151‑1 (MQ=255)
                                                                                                                                                gcggaGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAAcc   <  1:298142/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AAACAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAACCC  >  NC_000913/276751‑277046

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑