Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913523,3350TG42.7% 123.7 / 19.0 82G25G (GGT→GGG) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base T (21/26);  new base G (35/0);  total (56/26)
Fisher's exact test for biased strand distribution p-value = 7.85e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGTGTGTCCGGGCGGGATGACTT‑C‑GGG‑CAA‑CCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCT  >  NC_000913/523191‑523483
                                                                                                                                                |                                                                                                                                                         
cTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGgtg                                                                                                                                                          >  1:779737/1‑146 (MQ=255)
     gtgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTg                                                                                                                                                 <  3:387477/150‑1 (MQ=255)
     gtgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTg                                                                                                                                                 <  1:28211/150‑1 (MQ=255)
      tgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTTTCCAGGGGTCGGCAGGGGgaagaaat                                                                                                                                                  >  1:1195040/1‑146 (MQ=255)
      tgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTgg                                                                                                                                                >  4:243807/1‑150 (MQ=255)
         cATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGgcgc                                                                                                                                             <  1:673977/150‑1 (MQ=255)
               attttTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCAcc                                                                                                                                      <  2:567079/151‑1 (MQ=255)
               attttTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCAcc                                                                                                                                      <  1:547881/151‑1 (MQ=255)
                tttttttttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTTCCCTTTTCCAGGGGTCGGCCGGGGGAAGAATTGGAGCGCCCAcc                                                                                                                                      >  1:651392/1‑150 (MQ=255)
                 ttttttttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCg                                                                                                                                     <  2:831995/150‑1 (MQ=255)
                  tttttttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCgg                                                                                                                                    <  1:327942/150‑1 (MQ=255)
                      tttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGTTCGGCAGGGGTAAGAATTTGCGCCCCCACCgg                                                                                                                                    >  1:654/1‑146 (MQ=255)
                        tgttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGGCGGGGTAAGAATTTGCGCGCCCACCggcgggg                                                                                                                               >  2:625619/1‑149 (MQ=255)
                             cTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAAATGGGGCGCCCACCggcgggggggg                                                                                                                           >  2:182548/1‑144 (MQ=255)
                                 cacaGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCT‑CGgtg                                                                                                                   >  1:625116/1‑151 (MQ=255)
                                 cacaGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGGAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑Cgg                                                                                                                     >  4:40979/1‑149 (MQ=255)
                                  acaGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGTTCGGCAGGGGTAAGAAATGGCGCGCCCACCCGCGTGGCGTGGT‑Cgg                                                                                                                     >  2:257398/1‑148 (MQ=255)
                                   caGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCCGGGGTCGGCAGGGGGAAGAATTTGCGCGCCCCCCCGCGGGGCGTGGT‑Cgggg                                                                                                                   >  2:838100/1‑147 (MQ=255)
                                   caGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGGCGGGGTAAGAATTTGCGCGCCCACCCGCGTGGGGTGCT‑CGgggt                                                                                                                  >  4:223719/1‑150 (MQ=255)
                                   caGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGgtgtg                                                                                                                 <  1:1043649/151‑1 (MQ=255)
                                   caGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCT‑Cgg                                                                                                                     >  4:56482/1‑147 (MQ=255)
                                    aGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGGAAGAATTGGCGTGCCCACCGGCGGGGCGGGCT‑CGgtg                                                                                                                   >  1:207190/1‑148 (MQ=255)
                                    aGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTTAGAATTGGCGCGCCCCCCGGGCGGGCGTGCT‑Cgg                                                                                                                     >  3:390452/1‑146 (MQ=255)
                                      ggATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGTGTGTca                                                                                                              <  2:331740/151‑2 (MQ=255)
                                         tatGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTTCCCTTTTCCAGGGGTCGGCCGGGGGAAGAATTTGCGGGGCCCCCGGCGTGGCGGGCT‑CCGGGTGTTcggg                                                                                                           >  2:638558/1‑151 (MQ=255)
                                         tatGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGTTCGGCAGGTGTAAGAAATGGCGCGCCCACCCGCGTGGGGTGGT‑CGGTGTGTCccgg                                                                                                           >  1:172773/1‑151 (MQ=255)
                                          atGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGGAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGGGTGTCcggg                                                                                                           >  4:80027/1‑150 (MQ=255)
                                           tGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGGGGCGTGCT‑CGGGGTGTCcgggcg                                                                                                         >  4:328857/1‑151 (MQ=255)
                                            gAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGGGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGTGTGTCcgggcg                                                                                                         <  1:181371/150‑1 (MQ=255)
                                            gAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTTTCCAGGGGTTTGCAGGGGTAAATAATTGTGCCCCCCCCGGCCGGGCGGGGT‑CCGTGTGTTccggcgg                                                                                                        >  2:145125/1‑151 (MQ=255)
                                               cAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGTGTGTCCGGGCGGGAt                                                                                                     <  2:311623/151‑1 (MQ=255)
                                                  tAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCAGGGGTAAGGATTTGGGCGCCCCCCCGCGTGGGGCGGT‑CGGGGTGGTcgggggggggg                                                                                                    >  2:688077/1‑146 (MQ=255)
                                                  tAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCAATGTTCAGGGGTCGGCAGGGGTAAAAATTGGCGCGCCCACCGGCGGGGCGGGCT‑CGGGGTGTGcgggggggggg                                                                                                    >  2:948886/1‑146 (MQ=255)
                                                   aaCTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGGGTGCT‑CGGGGGGTCCGGGGGGGGTGACt                                                                                                 >  1:858547/1‑151 (MQ=255)
                                                     cTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGGGTGTCCGGGCGGGGTGACtt                                                                                                >  4:149188/1‑150 (MQ=255)
                                                      tAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGGTCGGGAGGGGGAAGAATTTGCGCGGCCCCCCGCGGGGCGGGGT‑CGGGGGTGCCGGGGGGGGTGAGtt                                                                                                >  2:337042/1‑149 (MQ=255)
                                                       aaaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTTTTCCGGGGTCGGCAGGTGTAAGAATTGGCGCGCCCCCCCGCGGGGCGGGCT‑CGGGGTGTCCGGGCGGGGTGGCTT‑ctg                                                                                            >  4:336250/1‑149 (MQ=255)
                                                       aaaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGTGTGTCCGGGCGGGATGACTT‑c                                                                                              <  1:1187779/149‑1 (MQ=255)
                                                       aaaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGGGTGGT‑CGGGGTGTCCGGGGGGGAGGAATT‑C‑gg                                                                                           >  1:873068/1‑151 (MQ=255)
                                                         aaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCAGGGGTAAGAATTTGCGCGCCCCCCCGCGGGGCGGGGT‑CGGGGGGGCCGGGGGGGGGGGCTT‑C‑ggg                                                                                          >  2:467452/1‑150 (MQ=255)
                                                         aaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGGGGCGTGCT‑CGGGGTGTCCGGGCGGGGTGACCT‑C‑ggg                                                                                          >  1:456623/1‑150 (MQ=255)
                                                              cACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGTGTGTCCGGGCGGGATGACTT‑C‑GGG‑CAA‑ccc                                                                                  <  3:404872/151‑1 (MQ=255)
                                                               aCTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGGGTGTCCGGGCGGGGTGACTT‑C‑GGG‑CAA‑acc                                                                                  >  4:286874/1‑150 (MQ=255)
                                                                 tttATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTTCCATTGTCCAGGGGTCGGCAGGTGTGAGAATTGGGGCGCCCCCCGGCGCGGCGTGCTCCGGGGGTTCCGGGGGGGTGGATTC‑G‑GGG‑CAA‑cc                                                                                   >  2:595992/1‑148 (MQ=255)
                                                                   tataTGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGGGGGCC‑CGGTGGGTCCGGGCGGGGTGGCTT‑T‑CGG‑CAA‑acc                                                                                  >  1:497811/1‑146 (MQ=255)
                                                                   tataTGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCT‑CGGGGTGTCCGGGCGGGGTGACTT‑C‑GGG‑GAA‑ccccg                                                                                >  1:670549/1‑146 (MQ=255)
                                                                   tataTGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTATCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGGGTGTCCGGGGGGGGTGACTT‑T‑GGG‑GAA‑ccccgga                                                                              >  1:1210174/1‑146 (MQ=255)
                                                                      atGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGTGTGTCCGGGCGGGATGACTT‑C‑GGG‑CAA‑CCCGGTAAAt                                                                           <  2:540184/150‑1 (MQ=255)
                                                                        gAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGTGTGTCCGGGCGGGATGACTT‑C‑GGG‑CAA‑CCCGGTAAATCCg                                                                        <  2:343006/151‑1 (MQ=255)
                                                                            ggAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCCTGTCCCGGGGTCGGGAGGTGTAAAAAATTGGGCCCCCCCCCGCGTGGGGCGCT‑CGGGGTGTCCGCGCGGGGGTGCTT‑CTGGG‑CAA‑CCCGGCaaactcc                                                                        >  2:868577/1‑144 (MQ=255)
                                                                            ggAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGTGTGTCCGGGCGGGATGACTT‑C‑GGG‑CAA‑CCCGGTAAATCCgctgc                                                                    <  1:12690/151‑1 (MQ=255)
                                                                              aaaaCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCAGGGGGAAGAATTGGCGCGCCCACCGGCGTGGGGTGCT‑CGGGGTGTTCGGGGGGGGTGAATTCC‑GGG‑CAA‑CCCGGAAAAcccc                                                                        >  1:704658/1‑145 (MQ=255)
                                                                                 aCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCT‑CCGTGTGGCCGGGGGGGGTGACTT‑T‑GGG‑GAA‑CCCGGGAAAACCGCCGCTgggg                                                               >  1:467722/1‑151 (MQ=255)
                                                                                   cAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCT‑CGGGGGGTCCCGGGGGGGTGACTT‑C‑GGG‑CAA‑CCCGGGAAATCCCCCGCTggggg                                                              >  1:313957/1‑150 (MQ=255)
                                                                                    aGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGTGTGTCCGGGCGGGATGACTT‑C‑GGG‑CAA‑CCCGGTAAATCCGCTGCTGGGGGCg                                                            <  2:1136268/151‑1 (MQ=255)
                                                                                      cggcgCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCT‑CGGTGTGTTCGGGCGGGGTGACTT‑T‑CGG‑CAA‑CCCGGTAAAACCCCTGCTggggg                                                              >  2:160000/1‑147 (MQ=255)
                                                                                      cggcgCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGTGTGTCCGGGCGGGATGACTT‑C‑GGG‑CAA‑CCCGGTAAATCCGCTGCTGGGGGCGaa                                                          <  1:477205/151‑1 (MQ=255)
                                                                                           cgTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGTGTAAGAATTGGCGCGCCCCCCGGCGTGGCGTGCT‑CGGGGTGTCCGGGCGGGGTGACTT‑C‑GGG‑CAA‑CCCGGGAAATCCCCTGCTGGGGGCGAAgggg                                                      >  4:110336/1‑148 (MQ=255)
                                                                                            gTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGTGTGTCCGGGCGGGATGACTT‑C‑GGG‑CAA‑CCCGGTAAATCCGCTGCTGGGGGCGAAGGtgct                                                    <  2:1194825/151‑1 (MQ=255)
                                                                                             tCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGGAAGAATTGGCGCGCCCACCGGCGTGGCGGGCT‑CGGGGTGTCCGGGAGGGATGGCTT‑C‑GGG‑CAA‑CCCGGTAAATCCCCTGCTCGGGGCGAAGGtgct                                                    >  1:751449/1‑150 (MQ=255)
                                                                                               aGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGGGTGTCCGGGCGGGGTGACTT‑C‑GGG‑CAA‑CCCGGGAAATCCGCTGCTGGGGGCGAAGGGGCTGcc                                                 >  4:90278/1‑151 (MQ=255)
                                                                                                gggAGATATGACTCAGTATGGCGGTCCCATTGTCCAAGGTTCGGGGGGGGTAAGAAATGGCGCGCCCACCGGCGCGGGGTGCT‑CGGGGCGTCCGGGGGGGGGGACCT‑C‑GGGCCAA‑CCGGGTAAA‑CCCCTGCTGGGGGCGGAGGGGCTGgcc                                                >  3:135707/1‑151 (MQ=255)
                                                                                                  gagaTATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGGGTGTTCGGGCGGGATGACTT‑C‑GGG‑CAA‑CCCGGTAAATCCGCTGCTGGGGGCGAAGGGGCTGccc                                                >  4:142770/1‑149 (MQ=255)
                                                                                                      tatGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCCGGGGTAAGAATTTGCGCGCCCACCGGCGTGGCGTGCT‑CGGGGTGTCCGGGGGGGGTGACTT‑C‑CGG‑CAA‑CCCGGGAAAACCGCTGCTGGGGGGGGAGGGGCTGCCCcgcg                                            >  2:1199569/1‑149 (MQ=255)
                                                                                                       atGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCCGGGGTAAGAATTGGCGCGCCCACCGGCGTGGGGGGCT‑CGGGGTGTCCCGGCGGGGTGACCT‑C‑GGG‑GAACCCCGGAAAATCC‑CTGCTGGGGGGGGAGGGGCTGccccggggg                                          >  1:703752/1‑144 (MQ=255)
                                                                                                           cTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCT‑CGGGGTGTTCGGGCGGGGTGACTT‑C‑GGG‑CAA‑CCCCGGAAATCCCCTGCTGGGGGCGGAGGTGGTGGCCCGCGAGACgg                                      >  2:536656/1‑150 (MQ=255)
                                                                                                             cAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGTGTGTCCGGGCGGGATGACTT‑C‑GGG‑CAA‑CCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGAcc                                   <  2:835643/151‑1 (MQ=255)
                                                                                                              aGTATGGCGGTCCCCTTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCC‑CGGGGGGGCCGGGGGGGGGGACTT‑T‑GGG‑GAA‑CCCCGGAAAAACCCCGCTGGGGGGGGAGGGGCTGCCCGGCGGGACGGAccc                                  >  2:1005901/1‑150 (MQ=255)
                                                                                                                 aTGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGTGTGTCCGGGGGGGGTGGCTT‑C‑GGG‑CAA‑CCCGGGAAATCCGCTGCTGGGGGGGGAGGGGGTGCCCCGCGAGACGGGccct                                 >  2:563829/1‑147 (MQ=255)
                                                                                                                   ggCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGTGTGTCCGGGCGGGATGACTT‑C‑GGG‑CAA‑CCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTgcg                              <  4:437818/150‑1 (MQ=255)
                                                                                                                            aTTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGTGTGTCCGGGCGGGATGACTT‑C‑GGG‑CAA‑CCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggcc                    >  1:847785/1‑151 (MQ=255)
                                                                                                                            aTTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGTGTGTCCGGGCGGGATGACTT‑C‑GGG‑CAA‑CCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggcc                    >  4:255270/1‑151 (MQ=255)
                                                                                                                              tGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGTGTGTCCGGGCGGGATGACTT‑C‑GGG‑CAA‑CCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggccc                   >  1:1117582/1‑150 (MQ=255)
                                                                                                                               gTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGTGTGTCCGGGCGGGATGACTT‑C‑GGG‑CAA‑CCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGc                 >  2:1160170/1‑151 (MQ=255)
                                                                                                                                 ccAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGTGTGTCCGGGCGGGATGACTT‑C‑GGG‑CAA‑CCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCt                <  1:1007972/150‑1 (MQ=255)
                                                                                                                                    gggTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGTGTGTCCGGGCGGGATGACTT‑C‑GGG‑CAA‑CCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGcc             >  2:613456/1‑150 (MQ=255)
                                                                                                                                      gTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGTGTGTCCGGGCGGGATGACTT‑C‑GGG‑CAA‑CCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGtt          <  1:504746/151‑1 (MQ=255)
                                                                                                                                        tCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGTGTGTCCGGGCGGGATGACTT‑C‑GGG‑CAA‑CCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTc         <  2:408514/150‑1 (MQ=255)
                                                                                                                                        tCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGTGTGTCCGGGCGGGATGACTT‑C‑GGG‑CAA‑CCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTc         >  1:670155/1‑150 (MQ=255)
                                                                                                                                         cGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGTGTGTCCGGGCGGGATGACTT‑C‑GGG‑CAA‑CCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCAt       >  1:778871/1‑151 (MQ=255)
                                                                                                                                         cGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGTGTGTCCGGGCGGGATGACTT‑C‑GGG‑CAA‑CCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCAt       >  3:220128/1‑151 (MQ=255)
                                                                                                                                           gCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGTGTGTCCGGGCGGGATGACTT‑C‑GGG‑CAA‑CCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATtc     <  1:864604/151‑1 (MQ=255)
                                                                                                                                              ggTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGTGTGTCCGGGCGGGATGACTT‑C‑GGG‑CAA‑CCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATtctct  <  2:1006573/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                         
CTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCT‑CGGTGTGTCCGGGCGGGATGACTT‑C‑GGG‑CAA‑CCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCT  >  NC_000913/523191‑523483

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