Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,224,4830CT52.6% 53.8 / inf 613S219F (TCC→TTC) metHhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base C (165/125);  new base T (169/153);  total (335/278)
Fisher's exact test for biased strand distribution p-value = 2.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 620 total aligned reads displayed.

GCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGT  >  NC_000913/4224334‑4224633
                                                                                                                                                     |                                                                                                                                                      
  tggtggAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGCCGCCTTCgg                                                                                                                                                     >  1:274501/1‑151 (MQ=255)
    gtggAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGgc                                                                                                                                                   <  6:203059/151‑1 (MQ=255)
         aGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACg                                                                                                                                              <  6:242638/151‑1 (MQ=255)
           gcggCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGct                                                                                                                                            <  7:1901138/149‑1 (MQ=255)
               cgcgGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTcc                                                                                                                                        >  8:735762/1‑151 (MQ=255)
                gcgGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCg                                                                                                                                       <  3:204608/151‑1 (MQ=255)
                   gATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGc                                                                                                                                    >  7:870997/1‑151 (MQ=255)
                   gATCTGATCATGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGAGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGAGTTGAGCTGCAGATTATTATCTCCGGCACCATCACCGACGCCTCCGCGAGCACGATCTCCGGGc                                                                                                                                    >  5:1278148/1‑151 (MQ=255)
                    aTCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGc                                                                                                                                    >  5:309235/1‑150 (MQ=255)
                       tGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGac                                                                                                                                >  8:1670327/1‑151 (MQ=255)
                        gATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAg                                                                                                                                  <  1:1306547/148‑1 (MQ=255)
                        gATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAg                                                                                                                                  <  1:1022015/148‑1 (MQ=255)
                        gATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCa                                                                                                                                   <  6:471920/147‑1 (MQ=255)
                        gATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCa                                                                                                                                   <  7:1289058/147‑1 (MQ=255)
                         aTCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGacca                                                                                                                              <  2:1459578/151‑1 (MQ=255)
                          tCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGaccac                                                                                                                             >  7:1057667/1‑151 (MQ=255)
                           ccTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGaccacc                                                                                                                            >  4:1812373/1‑151 (MQ=255)
                             tGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCg                                                                                                                           >  6:1203521/1‑150 (MQ=255)
                             tGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCg                                                                                                                           >  8:1830092/1‑150 (MQ=255)
                              gATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGa                                                                                                                          <  6:798088/150‑1 (MQ=255)
                              gATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGaa                                                                                                                         <  5:561434/151‑1 (MQ=255)
                                 tGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGc                                                                                                                       >  7:739472/1‑150 (MQ=255)
                                 tGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGc                                                                                                                       >  7:1541118/1‑150 (MQ=255)
                                 tGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGc                                                                                                                       >  2:794830/1‑150 (MQ=255)
                                  gAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCAt                                                                                                                     >  1:1419330/1‑151 (MQ=255)
                                     aCCGTTTTCGACACCCTTAACGCGAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCAtttt                                                                                                                  <  2:1829664/151‑1 (MQ=255)
                                     aCCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCAtttt                                                                                                                  <  6:866929/151‑1 (MQ=255)
                                        gTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTAc                                                                                                                >  2:1709286/1‑150 (MQ=255)
                                          tttCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAAc                                                                                                             >  4:1931972/1‑151 (MQ=255)
                                          tttCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAAc                                                                                                             >  3:2160979/1‑151 (MQ=255)
                                             cGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCa                                                                                                          <  6:532530/151‑1 (MQ=255)
                                               acacCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCAtt                                                                                                        >  3:106050/1‑151 (MQ=255)
                                                cacCCTTCACGCCAAAGCGGCGGTATTTGCGCTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTg                                                                                                       <  5:2171845/151‑1 (MQ=255)
                                                cacCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTg                                                                                                       <  4:991820/151‑1 (MQ=255)
                                                 acCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTgc                                                                                                      >  7:362371/1‑151 (MQ=255)
                                                  cccTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTgcg                                                                                                     <  1:869801/151‑1 (MQ=255)
                                                  cccTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTgcg                                                                                                     >  1:1927098/1‑151 (MQ=255)
                                                  cccTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTgcg                                                                                                     >  1:1751096/1‑151 (MQ=255)
                                                   ccTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTgcgc                                                                                                    >  1:485351/1‑151 (MQ=255)
                                                     ttAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGcc                                                                                                   >  2:553579/1‑150 (MQ=255)
                                                      tAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGCCGCCTCCGGGCGCCCGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCAc                                                                                                 >  3:843807/1‑151 (MQ=255)
                                                       aaCGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACg                                                                                                >  8:475577/1‑151 (MQ=255)
                                                        aCGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGc                                                                                               >  2:2104698/1‑151 (MQ=255)
                                                          gCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCg                                                                                             >  8:602658/1‑151 (MQ=255)
                                                           ccAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGTCAGACCACCGAAGCATTTTACAACTCATTGCGACacgacg                                                                                             >  7:1722201/1‑150 (MQ=255)
                                                              aaGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCACCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCAATTTACAACTCATTGCGCCACGCCGAAGc                                                                                         >  1:760548/1‑151 (MQ=255)
                                                               aGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGct                                                                                        <  8:919287/151‑1 (MQ=255)
                                                                 cggcggTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGctct                                                                                      <  4:674289/151‑1 (MQ=255)
                                                                 cggcggTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGGTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGctct                                                                                      <  3:1160185/151‑1 (MQ=255)
                                                                  ggcggTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTg                                                                                     <  5:903548/151‑1 (MQ=255)
                                                                  ggcggTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTg                                                                                     >  3:501478/1‑151 (MQ=255)
                                                                    cggTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGAc                                                                                   <  3:1623333/151‑1 (MQ=255)
                                                                      gTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCt                                                                                 <  7:979338/151‑1 (MQ=255)
                                                                      gTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCt                                                                                 >  5:2052502/1‑151 (MQ=255)
                                                                      gTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCt                                                                                 <  4:2138679/151‑1 (MQ=255)
                                                                       tATTTGCGGTGAAAACGGAGTTTGCAGCGCTGGGCGTTGAGTTGCCGAGTATGCTCTCCGGCCCCATCACCGCCGCCTCCGGGCGCACGCTCTCCGGGCAGCCGACCGCCGGATTTTACAACTCATTGCGCCACGCCGAAGCTCTGcccgt                                                                                >  7:1399430/1‑149 (MQ=255)
                                                                       tATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCCCCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCtt                                                                                >  5:1938583/1‑151 (MQ=255)
                                                                           tGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGc                                                                            >  5:1144552/1‑151 (MQ=255)
                                                                              ggTGAAAACGGAGTTTGAAGGGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTg                                                                         >  4:1105599/1‑151 (MQ=255)
                                                                                  aagaCGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACt                                                                     <  2:463972/148‑1 (MQ=255)
                                                                                     aCGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACtgtg                                                                  <  6:1428713/151‑1 (MQ=255)
                                                                                     aCGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACtgtg                                                                  >  4:1004673/1‑151 (MQ=255)
                                                                                       ggAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTgcg                                                                <  4:1705179/151‑1 (MQ=255)
                                                                                        gAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTgcgc                                                               >  4:1569891/1‑151 (MQ=255)
                                                                                         aGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCt                                                              >  8:1132905/1‑151 (MQ=255)
                                                                                               aaGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGcc                                                        <  8:1748703/151‑1 (MQ=255)
                                                                                                 gcgcTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCg                                                      <  4:1465648/151‑1 (MQ=255)
                                                                                                   gcTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGAt                                                    <  5:634852/151‑1 (MQ=255)
                                                                                                     ggggCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGa                                                  <  4:2042587/150‑1 (MQ=255)
                                                                                                      gggCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGaa                                                 <  3:1473004/151‑1 (MQ=255)
                                                                                                        gCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACt                                               >  2:89127/1‑151 (MQ=255)
                                                                                                         cGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTg                                              >  5:1537355/1‑151 (MQ=255)
                                                                                                         cGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTg                                              >  7:2131506/1‑151 (MQ=255)
                                                                                                          gTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTgc                                             <  3:336259/151‑1 (MQ=255)
                                                                                                           ttGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTgcg                                            >  5:1779382/1‑151 (MQ=255)
                                                                                                               gCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGcc                                          >  2:124458/1‑149 (MQ=255)
                                                                                                                 tGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGt                                       >  4:1329836/1‑150 (MQ=255)
                                                                                                                 tGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGt                                       >  7:505737/1‑150 (MQ=255)
                                                                                                                  gCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGt                                       <  5:818335/149‑1 (MQ=255)
                                                                                                                   ccGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACg                                    <  8:144776/151‑1 (MQ=255)
                                                                                                                    cGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGt                                   >  6:669755/1‑151 (MQ=255)
                                                                                                                        tATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCa                                <  8:2154619/150‑1 (MQ=255)
                                                                                                                         aTGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAg                               <  3:1337095/150‑1 (MQ=255)
                                                                                                                           gATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGa                             <  4:409131/150‑1 (MQ=255)
                                                                                                                            aTCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGc                           >  2:279993/1‑151 (MQ=255)
                                                                                                                            aTCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGc                           >  8:328914/1‑151 (MQ=255)
                                                                                                                               tcCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGt                        >  4:831073/1‑151 (MQ=255)
                                                                                                                                ccGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTc                       <  2:671503/151‑1 (MQ=255)
                                                                                                                                ccGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGc                           <  1:625956/147‑1 (MQ=255)
                                                                                                                                ccGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTc                       <  4:11447/151‑1 (MQ=255)
                                                                                                                                     aCCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACgg                   >  2:192281/1‑150 (MQ=255)
                                                                                                                                     aCCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGa                  <  6:1852620/151‑1 (MQ=255)
                                                                                                                                           aCCGCCGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCCCGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCgg            >  1:2112910/1‑151 (MQ=255)
                                                                                                                                            ccGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCgg            >  6:1593284/1‑150 (MQ=255)
                                                                                                                                            ccGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCgg            >  6:1024787/1‑150 (MQ=255)
                                                                                                                                             cGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGaa          <  1:1505634/151‑1 (MQ=255)
                                                                                                                                             cGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTCTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGa           >  7:620672/1‑150 (MQ=255)
                                                                                                                                                cGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACTCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGc       >  3:1964556/1‑151 (MQ=255)
                                                                                                                                                 gCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGACTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCt      >  2:1061756/1‑151 (MQ=255)
                                                                                                                                                     ccGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGt  >  4:299414/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGT  >  NC_000913/4224334‑4224633

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