Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913523,3260TG45.2% 86.7 / 12.8 62G22G (GGT→GGGrhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base T (12/22);  new base G (28/0);  total (40/22)
Fisher's exact test for biased strand distribution p-value = 1.72e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCT  >  NC_000913/523186‑523469
                                                                                                                                             |                                                                                                                                                 
cTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTTGGCAgggg                                                                                                                                         >  1:982904/1‑149 (MQ=255)
  ggTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGGTCGGCAgggg                                                                                                                                         >  1:582134/1‑147 (MQ=255)
   gTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCgggcgggg                                                                                                                                         >  1:799452/1‑142 (MQ=255)
         tgtgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTTGGGCGGGGTAAGaaat                                                                                                                                 >  1:871153/1‑148 (MQ=255)
          gtgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTg                                                                                                                                <  1:373051/150‑1 (MQ=255)
          gtgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTg                                                                                                                                <  1:300363/150‑1 (MQ=255)
            gtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGGAGGTGTAAGaaattg                                                                                                                               >  1:372336/1‑145 (MQ=255)
              cATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTTTCCAGGGGTCGGCAGGGGTAAGAATTgg                                                                                                                               >  1:1163338/1‑147 (MQ=255)
              cATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGgcgcg                                                                                                                           <  1:98534/151‑1 (MQ=255)
                 tttattttTTTTTGTTGCTAACTCAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGGCCCAATGTGCAGGGGTTGGCAGGGGGTAGAATTTGGGCGgcc                                                                                                                        >  1:1133731/1‑151 (MQ=255)
                        ttttttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGAGCGTCAGGGAGATATGACTTAGTATGGCGGTCCCATTTTCCAAGGTTGGGGAGGGGTTAAAATTTGGGCGaccacc                                                                                                                     >  1:193385/1‑147 (MQ=255)
                          ttttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCggcgt                                                                                                                <  1:459755/150‑1 (MQ=255)
                            ttgttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCAATTCCCCAGGGTGCGCAGGGGGAAGACATCGCGCGCCCAc                                                                                                                      >  2:1201541/1‑142 (MQ=255)
                               ttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCCTTTCCAGGGGGTGGCCGGAGTTAGAAATGGGGCCGCCCCCcgcgcggcgt                                                                                                           >  1:227348/1‑150 (MQ=255)
                               ttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGATCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGgcgtggcgtg                                                                                                          >  2:29341/1‑151 (MQ=255)
                                tgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTTCCACTGTCCAGGGGTTGGGAGGGGTAAGAAATGGTGCGGCCCCCggccgggcc                                                                                                            >  1:858529/1‑142 (MQ=255)
                                tgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCAGGTGTAAGAATTGGCGCGCCCACCGgcggggcgtg                                                                                                          >  2:769375/1‑150 (MQ=255)
                                tgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGTGTAAGAATTGGCGCGCCCACCcgcgggggggg                                                                                                          >  1:328921/1‑146 (MQ=255)
                                tgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCggcgtggcg                                                                                                            >  2:273898/1‑148 (MQ=255)
                                 gCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGc                                                                                                         <  1:1189314/150‑1 (MQ=255)
                                   tAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCATGGTCAGGGGTCGGCAGGTGGAAGAAATGGGGCGCCCACCAgcgttgggtggt                                                                                                        >  1:569919/1‑147 (MQ=255)
                                    aaCACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTc                                                                                                       <  1:145326/149‑1 (MQ=255)
                                    aaCACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCg                                                                                                      <  1:840461/150‑1 (MQ=255)
                                    aaCACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGGGGGCTCgg                                                                                                     >  2:1050659/1‑151 (MQ=255)
                                     acacAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCCGGGGTCGGCAGGTGTAAGAAATGGCGCGCCCACCCGCGGGGGGGGCTCgg                                                                                                     >  1:470839/1‑150 (MQ=255)
                                      cacaGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGgtg                                                                                                   <  1:119819/151‑1 (MQ=255)
                                      cacaGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGGAAGAAATGGCGCGCCCCCCcgcgtggggggccccg                                                                                                     >  1:1080036/1‑140 (MQ=255)
                                       acaGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCCCCGGCGTGGCGTGCTCgggg                                                                                                   >  2:751786/1‑148 (MQ=255)
                                        caGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCCCGgggt                                                                                                  >  2:558336/1‑150 (MQ=255)
                                         cgggATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGTTCGGCAGGGGTTAGAATTGGCGCGCCCACCGGCGGGGCGGGGTCGgggtg                                                                                                 >  1:369006/2‑150 (MQ=255)
                                         aGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGAGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGGCTGCAGGGGGAAAAATTTCGGCGCCCCCCCGCGTGGCGGGCGCGgtg                                                                                                   >  1:866019/1‑148 (MQ=255)
                                             atatGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCcgg                                                                                            <  2:142547/151‑1 (MQ=255)
                                                 gAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCcgggcgg                                                                                        <  2:1058577/151‑1 (MQ=255)
                                                      aTAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTTGCGCGCCCACCGGCGTGGCGGGCTCGGGGGGTCCGGGGGGGGTg                                                                                    >  2:610591/1‑150 (MQ=255)
                                                         aCTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAAATGGCGCGCCCCCCGGCGGGGGGGGCCCGGGGGGTCccggggggg                                                                                      >  1:141736/1‑144 (MQ=255)
                                                          cTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGAGCGTCAGGGAGATATGACTCAGTATGGCGGTACCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGGGCGCCCACCCGCGTGGCGGGCGCGGGGGGTTcggggggg                                                                                       >  1:122694/1‑143 (MQ=255)
                                                            aaaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCAGGGGGAAGAATTTGGGCGCCCACCGGCGTGGCGTGCTCGGGGTGTCCGGGGGGGGAGACTTCgg                                                                           >  2:961160/1‑151 (MQ=255)
                                                               aGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCCGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTCCGGGCGGGGTGACTTCggg                                                                          >  2:269339/1‑149 (MQ=255)
                                                               aGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCa                                                                        <  2:334735/151‑1 (MQ=255)
                                                                gggCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGGAGGGGTAAGAATTGGTGCGCCCACCGGCGTGGCGGGCTTGGTGGGTCCGGGGGGGAGGACTTCgggga                                                                        >  2:538727/1‑148 (MQ=255)
                                                                 ggCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAAc                                                                      <  2:433694/151‑1 (MQ=255)
                                                                      tttATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCCGGGGTAAAAATTGGCGCGCCCCCCGGCCGTGCGTGCGCGGGGGGTTCGGGCGGGATGACTTCCGGCAAccccg                                                                  >  2:968216/1‑148 (MQ=255)
                                                                         atatGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAAATGGCGCGCCCCCCGGCGGGGCGTGCTCGGGGGGTCCGGGGGGGGTGACTTCGGGGCACCCCggaa                                                                >  1:1150950/1‑148 (MQ=255)
                                                                              aGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCgc                                                         <  2:1077243/151‑1 (MQ=255)
                                                                              aGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCTCGGGGTGTCCCGGCGGGGGGACTTCGGGGAACCCGGGAAAACCgc                                                         >  2:862061/1‑151 (MQ=255)
                                                                                    aaaCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCAGGGGTAAGAAATGGCGCGCCCACCCGCGGGGGGGGGTCGGGGGGGCCGGGGGGGGGGAACTCGGGGAACCCCGGaaaaaccc                                                         >  1:917911/1‑140 (MQ=255)
                                                                                    aaaCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTgg                                                   <  1:287747/151‑1 (MQ=255)
                                                                                        cAGAGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGGAAGAATTTGCGCGCCCACCGGCGTGGCGGTGTCGGGGTTGCCGGGGGGGGTGACCTCGGGGAACCCCGGAAAACCCCTGCTgggggg                                               >  1:265480/1‑150 (MQ=255)
                                                                                                 gTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCCGCGGGGCGGGCTCGGTGTGTCCGGGGGGGGGGACCTCGGGCAACCCCGGAAAACCCCCGCTGGGGGGGGAgggg                                        >  1:749229/1‑147 (MQ=255)
                                                                                                   cAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGtgctg                                     <  2:294599/150‑1 (MQ=255)
                                                                                                   cAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTTCGGGGGGGGTGGCTTTCGGGGACCCCGGAAAACCCCCGCTGGGGGGGGAgggg                                        >  1:877422/1‑145 (MQ=255)
                                                                                                       gagaTATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCgg                                <  1:996074/151‑1 (MQ=255)
                                                                                                              gACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACg                         <  1:462068/151‑1 (MQ=255)
                                                                                                                  cAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGAcc                     <  2:247319/151‑1 (MQ=255)
                                                                                                                   aGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCt                    <  2:632221/151‑1 (MQ=255)
                                                                                                                    gTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCtt                   <  1:105020/151‑1 (MQ=255)
                                                                                                                    gTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGGAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCtt                   >  2:229458/1‑151 (MQ=255)
                                                                                                                         gCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCt              <  2:255245/151‑1 (MQ=255)
                                                                                                                          cGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCt              >  1:303265/1‑150 (MQ=255)
                                                                                                                            gTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcc           >  2:911075/1‑151 (MQ=255)
                                                                                                                             tCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCCCGGCGTGCCCGGGCGGGACGACTTCCGGCAACCCGGCAAATCCGCCGCTGGGGGCGCAGGCGCTGCCCGGCGAGACGGACCTTGCGCTgccc          >  1:15133/1‑151 (MQ=255)
                                                                                                                                       ccAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCt  <  2:694215/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
CTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCT  >  NC_000913/523186‑523469

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