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→GGGrhsDRhs 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

CTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCT  >  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)
                                 cacaGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCTCGgtg                                                                                                                   >  1:625116/1‑151 (MQ=255)
                                 cacaGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGGAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCgg                                                                                                                     >  4:40979/1‑149 (MQ=255)
                                  acaGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGTTCGGCAGGGGTAAGAAATGGCGCGCCCACCCGCGTGGCGTGGTCgg                                                                                                                     >  2:257398/1‑148 (MQ=255)
                                   caGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCCGGGGTCGGCAGGGGGAAGAATTTGCGCGCCCCCCCGCGGGGCGTGGTCgggg                                                                                                                   >  2:838100/1‑147 (MQ=255)
                                   caGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGGCGGGGTAAGAATTTGCGCGCCCACCCGCGTGGGGTGCTCGgggt                                                                                                                  >  4:223719/1‑150 (MQ=255)
                                   caGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGgtgtg                                                                                                                 <  1:1043649/151‑1 (MQ=255)
                                   caGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCTCgg                                                                                                                     >  4:56482/1‑147 (MQ=255)
                                    aGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGGAAGAATTGGCGTGCCCACCGGCGGGGCGGGCTCGgtg                                                                                                                   >  1:207190/1‑148 (MQ=255)
                                    aGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTTAGAATTGGCGCGCCCCCCGGGCGGGCGTGCTCgg                                                                                                                     >  3:390452/1‑146 (MQ=255)
                                      ggATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTca                                                                                                              <  2:331740/151‑2 (MQ=255)
                                         tatGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTTCCCTTTTCCAGGGGTCGGCCGGGGGAAGAATTTGCGGGGCCCCCGGCGTGGCGGGCTCCGGGTGTTcggg                                                                                                           >  2:638558/1‑151 (MQ=255)
                                         tatGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGTTCGGCAGGTGTAAGAAATGGCGCGCCCACCCGCGTGGGGTGGTCGGTGTGTCccgg                                                                                                           >  1:172773/1‑151 (MQ=255)
                                          atGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGGAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTCcggg                                                                                                           >  4:80027/1‑150 (MQ=255)
                                           tGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGGGGCGTGCTCGGGGTGTCcgggcg                                                                                                         >  4:328857/1‑151 (MQ=255)
                                            gAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGGGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCcgggcg                                                                                                         <  1:181371/150‑1 (MQ=255)
                                            gAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTTTCCAGGGGTTTGCAGGGGTAAATAATTGTGCCCCCCCCGGCCGGGCGGGGTCCGTGTGTTccggcgg                                                                                                        >  2:145125/1‑151 (MQ=255)
                                               cAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAt                                                                                                     <  2:311623/151‑1 (MQ=255)
                                                  tAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCAGGGGTAAGGATTTGGGCGCCCCCCCGCGTGGGGCGGTCGGGGTGGTcgggggggggg                                                                                                    >  2:688077/1‑146 (MQ=255)
                                                  tAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCAATGTTCAGGGGTCGGCAGGGGTAAAAATTGGCGCGCCCACCGGCGGGGCGGGCTCGGGGTGTGcgggggggggg                                                                                                    >  2:948886/1‑146 (MQ=255)
                                                   aaCTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGGGTGCTCGGGGGGTCCGGGGGGGGTGACt                                                                                                 >  1:858547/1‑151 (MQ=255)
                                                     cTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTCCGGGCGGGGTGACtt                                                                                                >  4:149188/1‑150 (MQ=255)
                                                      tAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGGTCGGGAGGGGGAAGAATTTGCGCGGCCCCCCGCGGGGCGGGGTCGGGGGTGCCGGGGGGGGTGAGtt                                                                                                >  2:337042/1‑149 (MQ=255)
                                                       aaaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTTTTCCGGGGTCGGCAGGTGTAAGAATTGGCGCGCCCCCCCGCGGGGCGGGCTCGGGGTGTCCGGGCGGGGTGGCTTctg                                                                                            >  4:336250/1‑149 (MQ=255)
                                                       aaaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTc                                                                                              <  1:1187779/149‑1 (MQ=255)
                                                       aaaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGGGTGGTCGGGGTGTCCGGGGGGGAGGAATTCgg                                                                                           >  1:873068/1‑151 (MQ=255)
                                                         aaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCAGGGGTAAGAATTTGCGCGCCCCCCCGCGGGGCGGGGTCGGGGGGGCCGGGGGGGGGGGCTTCggg                                                                                          >  2:467452/1‑150 (MQ=255)
                                                         aaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGGGGCGTGCTCGGGGTGTCCGGGCGGGGTGACCTCggg                                                                                          >  1:456623/1‑150 (MQ=255)
                                                              cACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAAccc                                                                                  <  3:404872/151‑1 (MQ=255)
                                                               aCTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTCCGGGCGGGGTGACTTCGGGCAAacc                                                                                  >  4:286874/1‑150 (MQ=255)
                                                                 tttATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTTCCATTGTCCAGGGGTCGGCAGGTGTGAGAATTGGGGCGCCCCCCGGCGCGGCGTGCTCCGGGGGTTCCGGGGGGGTGGATTCGGGGCAAcc                                                                                   >  2:595992/1‑148 (MQ=255)
                                                                   tataTGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGGGGGCCCGGTGGGTCCGGGCGGGGTGGCTTTCGGCAAacc                                                                                  >  1:497811/1‑146 (MQ=255)
                                                                   tataTGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCTCGGGGTGTCCGGGCGGGGTGACTTCGGGGAAccccg                                                                                >  1:670549/1‑146 (MQ=255)
                                                                   tataTGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTATCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTCCGGGGGGGGTGACTTTGGGGAAccccgga                                                                              >  1:1210174/1‑146 (MQ=255)
                                                                      atGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAAt                                                                           <  2:540184/150‑1 (MQ=255)
                                                                        gAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCg                                                                        <  2:343006/151‑1 (MQ=255)
                                                                            ggAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCCTGTCCCGGGGTCGGGAGGTGTAAAAAATTGGGCCCCCCCCCGCGTGGGGCGCTCGGGGTGTCCGCGCGGGGGTGCTTCTGGGCAACCCGGCaaactcc                                                                        >  2:868577/1‑144 (MQ=255)
                                                                            ggAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCgctgc                                                                    <  1:12690/151‑1 (MQ=255)
                                                                              aaaaCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCAGGGGGAAGAATTGGCGCGCCCACCGGCGTGGGGTGCTCGGGGTGTTCGGGGGGGGTGAATTCCGGGCAACCCGGAAAAcccc                                                                        >  1:704658/1‑145 (MQ=255)
                                                                                 aCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCTCCGTGTGGCCGGGGGGGGTGACTTTGGGGAACCCGGGAAAACCGCCGCTgggg                                                               >  1:467722/1‑151 (MQ=255)
                                                                                   cAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCTCGGGGGGTCCCGGGGGGGTGACTTCGGGCAACCCGGGAAATCCCCCGCTggggg                                                              >  1:313957/1‑150 (MQ=255)
                                                                                    aGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCg                                                            <  2:1136268/151‑1 (MQ=255)
                                                                                      cggcgCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCTCGGTGTGTTCGGGCGGGGTGACTTTCGGCAACCCGGTAAAACCCCTGCTggggg                                                              >  2:160000/1‑147 (MQ=255)
                                                                                      cggcgCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGaa                                                          <  1:477205/151‑1 (MQ=255)
                                                                                           cgTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGTGTAAGAATTGGCGCGCCCCCCGGCGTGGCGTGCTCGGGGTGTCCGGGCGGGGTGACTTCGGGCAACCCGGGAAATCCCCTGCTGGGGGCGAAgggg                                                      >  4:110336/1‑148 (MQ=255)
                                                                                            gTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGtgct                                                    <  2:1194825/151‑1 (MQ=255)
                                                                                             tCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGGAAGAATTGGCGCGCCCACCGGCGTGGCGGGCTCGGGGTGTCCGGGAGGGATGGCTTCGGGCAACCCGGTAAATCCCCTGCTCGGGGCGAAGGtgct                                                    >  1:751449/1‑150 (MQ=255)
                                                                                               aGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTCCGGGCGGGGTGACTTCGGGCAACCCGGGAAATCCGCTGCTGGGGGCGAAGGGGCTGcc                                                 >  4:90278/1‑151 (MQ=255)
                                                                                                gggAGATATGACTCAGTATGGCGGTCCCATTGTCCAAGGTTCGGGGGGGGTAAGAAATGGCGCGCCCACCGGCGCGGGGTGCTCGGGGCGTCCGGGGGGGGGGACCTCGGGCCAACCGGGTAAACCCCTGCTGGGGGCGGAGGGGCTGgcc                                                >  3:135707/1‑151 (MQ=255)
                                                                                                  gagaTATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTTCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGGGCTGccc                                                >  4:142770/1‑149 (MQ=255)
                                                                                                      tatGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCCGGGGTAAGAATTTGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTCCGGGGGGGGTGACTTCCGGCAACCCGGGAAAACCGCTGCTGGGGGGGGAGGGGCTGCCCcgcg                                            >  2:1199569/1‑149 (MQ=255)
                                                                                                       atGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCCGGGGTAAGAATTGGCGCGCCCACCGGCGTGGGGGGCTCGGGGTGTCCCGGCGGGGTGACCTCGGGGAACCCCGGAAAATCCCTGCTGGGGGGGGAGGGGCTGccccggggg                                          >  1:703752/1‑144 (MQ=255)
                                                                                                           cTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCTCGGGGTGTTCGGGCGGGGTGACTTCGGGCAACCCCGGAAATCCCCTGCTGGGGGCGGAGGTGGTGGCCCGCGAGACgg                                      >  2:536656/1‑150 (MQ=255)
                                                                                                             cAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGAcc                                   <  2:835643/151‑1 (MQ=255)
                                                                                                              aGTATGGCGGTCCCCTTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCCCGGGGGGGCCGGGGGGGGGGACTTTGGGGAACCCCGGAAAAACCCCGCTGGGGGGGGAGGGGCTGCCCGGCGGGACGGAccc                                  >  2:1005901/1‑150 (MQ=255)
                                                                                                                 aTGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGGGGGGTGGCTTCGGGCAACCCGGGAAATCCGCTGCTGGGGGGGGAGGGGGTGCCCCGCGAGACGGGccct                                 >  2:563829/1‑147 (MQ=255)
                                                                                                                   ggCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTgcg                              <  4:437818/150‑1 (MQ=255)
                                                                                                                            aTTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggcc                    >  1:847785/1‑151 (MQ=255)
                                                                                                                            aTTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggcc                    >  4:255270/1‑151 (MQ=255)
                                                                                                                              tGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggccc                   >  1:1117582/1‑150 (MQ=255)
                                                                                                                               gTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGc                 >  2:1160170/1‑151 (MQ=255)
                                                                                                                                 ccAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCt                <  1:1007972/150‑1 (MQ=255)
                                                                                                                                    gggTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGcc             >  2:613456/1‑150 (MQ=255)
                                                                                                                                      gTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGtt          <  1:504746/151‑1 (MQ=255)
                                                                                                                                        tCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTc         <  2:408514/150‑1 (MQ=255)
                                                                                                                                        tCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTc         >  1:670155/1‑150 (MQ=255)
                                                                                                                                         cGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCAt       >  1:778871/1‑151 (MQ=255)
                                                                                                                                         cGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCAt       >  3:220128/1‑151 (MQ=255)
                                                                                                                                           gCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATtc     <  1:864604/151‑1 (MQ=255)
                                                                                                                                              ggTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATtctct  <  2:1006573/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                         
CTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCT  >  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: