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→GGG) rhsDRhs 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

CTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGAC‑TTC‑GGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCT  >  NC_000913/523186‑523469
                                                                                                                                             |                                                                                                                                                 
cTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGGTTGGCAgggg                                                                                                                                         >  1:982904/1‑149 (MQ=255)
  ggTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCCGGGGTCGGCAgggg                                                                                                                                         >  1:582134/1‑147 (MQ=255)
   gTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑CTTGTCCAGGGGTCgggcgggg                                                                                                                                         >  1:799452/1‑142 (MQ=255)
         tgtgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGGTTGGGCGGGGTAAGaaat                                                                                                                                 >  1:871153/1‑148 (MQ=255)
          gtgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGTTCGGCAGGTGTAAGAATTg                                                                                                                                <  1:373051/150‑1 (MQ=255)
          gtgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGTTCGGCAGGTGTAAGAATTg                                                                                                                                <  1:300363/150‑1 (MQ=255)
            gtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGTTCGGGAGGTGTAAGaaattg                                                                                                                               >  1:372336/1‑145 (MQ=255)
              cATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTTTCCAGGGGTCGGCAGGGGTAAGAATTgg                                                                                                                               >  1:1163338/1‑147 (MQ=255)
              cATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGTTCGGCAGGTGTAAGAATTGgcgcg                                                                                                                           <  1:98534/151‑1 (MQ=255)
                 tttattttTTTTTGTTGCTAACTCAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGGCCC‑AATGTGCAGGGGTTGGCAGGGGGTAGAATTTGGGCGgcc                                                                                                                        >  1:1133731/1‑151 (MQ=255)
                        ttttttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGAGCGTCAGGGAGATATGACTTAGTATGGCGGTCCC‑ATTTTCCAAGGTTGGGGAGGGGTTAAAATTTGGGCGaccacc                                                                                                                     >  1:193385/1‑147 (MQ=255)
                          ttttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCggcgt                                                                                                                <  1:459755/150‑1 (MQ=255)
                            ttgttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCAATTCCCCAGGGTGC‑GCAGGGGGAAGACATCGCGCGCCCAc                                                                                                                      >  2:1201541/1‑142 (MQ=255)
                               ttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑CCTTTCCAGGGGGTGGCCGGAGTTAGAAATGGGGCCGCCCCCcgcgcggcgt                                                                                                           >  1:227348/1‑150 (MQ=255)
                               ttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGATCCC‑ATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGgcgtggcgtg                                                                                                          >  2:29341/1‑151 (MQ=255)
                                tgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTTCC‑ACTGTCCAGGGGTTGGGAGGGGTAAGAAATGGTGCGGCCCCCggccgggcc                                                                                                            >  1:858529/1‑142 (MQ=255)
                                tgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑CTTGTCCAGGGGTCGGCAGGTGTAAGAATTGGCGCGCCCACCGgcggggcgtg                                                                                                          >  2:769375/1‑150 (MQ=255)
                                tgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGGTCGGCAGGTGTAAGAATTGGCGCGCCCACCcgcgggggggg                                                                                                          >  1:328921/1‑146 (MQ=255)
                                tgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCggcgtggcg                                                                                                            >  2:273898/1‑148 (MQ=255)
                                 gCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGc                                                                                                         <  1:1189314/150‑1 (MQ=255)
                                   tAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑CATGGTCAGGGGTCGGCAGGTGGAAGAAATGGGGCGCCCACCAgcgttgggtggt                                                                                                        >  1:569919/1‑147 (MQ=255)
                                    aaCACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTc                                                                                                       <  1:145326/149‑1 (MQ=255)
                                    aaCACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCg                                                                                                      <  1:840461/150‑1 (MQ=255)
                                    aaCACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGGGGGCTCgg                                                                                                     >  2:1050659/1‑151 (MQ=255)
                                     acacAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑CTTGTCCCGGGGTCGGCAGGTGTAAGAAATGGCGCGCCCACCCGCGGGGGGGGCTCgg                                                                                                     >  1:470839/1‑150 (MQ=255)
                                      cacaGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGgtg                                                                                                   <  1:119819/151‑1 (MQ=255)
                                      cacaGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGGTCGGCAGGGGGAAGAAATGGCGCGCCCCCCcgcgtggggggccccg                                                                                                     >  1:1080036/1‑140 (MQ=255)
                                       acaGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCCCCGGCGTGGCGTGCTCgggg                                                                                                   >  2:751786/1‑148 (MQ=255)
                                        caGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCCCGgggt                                                                                                  >  2:558336/1‑150 (MQ=255)
                                         cgggATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑CTTGTCCAGGGTTCGGCAGGGGTTAGAATTGGCGCGCCCACCGGCGGGGCGGGGTCGgggtg                                                                                                 >  1:369006/2‑150 (MQ=255)
                                         aGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGAGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑CTTGTCCAGGGGGCTGCAGGGGGAAAAATTTCGGCGCCCCCCCGCGTGGCGGGCGCGgtg                                                                                                   >  1:866019/1‑148 (MQ=255)
                                             atatGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCcgg                                                                                            <  2:142547/151‑1 (MQ=255)
                                                 gAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCcgggcgg                                                                                        <  2:1058577/151‑1 (MQ=255)
                                                      aTAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGGTCGGCAGGGGTAAGAATTTGCGCGCCCACCGGCGTGGCGGGCTCGGGGGGTCCGGGGGGGGTg                                                                                    >  2:610591/1‑150 (MQ=255)
                                                         aCTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGGTCGGCAGGGGTAAGAAATGGCGCGCCCCCCGGCGGGGGGGGCCCGGGGGGTCccggggggg                                                                                      >  1:141736/1‑144 (MQ=255)
                                                          cTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGAGCGTCAGGGAGATATGACTCAGTATGGCGGTACC‑ATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGGGCGCCCACCCGCGTGGCGGGCGCGGGGGGTTcggggggg                                                                                       >  1:122694/1‑143 (MQ=255)
                                                            aaaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑CTTGTCCAGGGGTCGGCAGGGGGAAGAATTTGGGCGCCCACCGGCGTGGCGTGCTCGGGGTGTCCGGGGGGGGAGAC‑TTC‑gg                                                                           >  2:961160/1‑151 (MQ=255)
                                                               aGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGTTCGGCCGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTCCGGGCGGGGTGAC‑TTC‑ggg                                                                          >  2:269339/1‑149 (MQ=255)
                                                               aGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGAC‑TTC‑GGGCa                                                                        <  2:334735/151‑1 (MQ=255)
                                                                gggCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGGTCGGGAGGGGTAAGAATTGGTGCGCCCACCGGCGTGGCGGGCTTGGTGGGTCCGGGGGGGAGGAC‑TTC‑gggga                                                                        >  2:538727/1‑148 (MQ=255)
                                                                 ggCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGAC‑TTC‑GGGCAAc                                                                      <  2:433694/151‑1 (MQ=255)
                                                                      tttATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGTTCGGCCGGGGTAAAAATTGGCGCGCCCCCCGGCCGTGCGTGCGCGGGGGGTTCGGGCGGGATGAC‑TTC‑CGGCAAccccg                                                                  >  2:968216/1‑148 (MQ=255)
                                                                         atatGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGGTCGGCAGGGGTAAGAAATGGCGCGCCCCCCGGCGGGGCGTGCTCGGGGGGTCCGGGGGGGGTGAC‑TTCGGGGCACCCCggaa                                                                >  1:1150950/1‑148 (MQ=255)
                                                                              aGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGAC‑TTC‑GGGCAACCCGGTAAATCCgc                                                         <  2:1077243/151‑1 (MQ=255)
                                                                              aGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCTCGGGGTGTCCCGGCGGGGGGAC‑TTC‑GGGGAACCCGGGAAAACCgc                                                         >  2:862061/1‑151 (MQ=255)
                                                                                    aaaCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑CTTGTCCAGGGGTCGGCAGGGGTAAGAAATGGCGCGCCCACCCGCGGGGGGGGGTCGGGGGGGCCGGGGGGGGGGAA‑CTC‑GGGGAACCCCGGaaaaaccc                                                         >  1:917911/1‑140 (MQ=255)
                                                                                    aaaCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGAC‑TTC‑GGGCAACCCGGTAAATCCGCTGCTgg                                                   <  1:287747/151‑1 (MQ=255)
                                                                                        cAGAGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGGTCGGCAGGGGGAAGAATTTGCGCGCCCACCGGCGTGGCGGTGTCGGGGTTGCCGGGGGGGGTGAC‑CTC‑GGGGAACCCCGGAAAACCCCTGCTgggggg                                               >  1:265480/1‑150 (MQ=255)
                                                                                                 gTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCCGCGGGGCGGGCTCGGTGTGTCCGGGGGGGGGGAC‑CTC‑GGGCAACCCCGGAAAACCCCCGCTGGGGGGGGAgggg                                        >  1:749229/1‑147 (MQ=255)
                                                                                                   cAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGAC‑TTC‑GGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGtgctg                                     <  2:294599/150‑1 (MQ=255)
                                                                                                   cAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGGTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTTCGGGGGGGGTGGCTTTC‑GGGGACCCCGG‑AAAACCCCCGCTGGGGGGGGAgggg                                        >  1:877422/1‑145 (MQ=255)
                                                                                                       gagaTATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGAC‑TTC‑GGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCgg                                <  1:996074/151‑1 (MQ=255)
                                                                                                              gACTCAGTATGGCGGTCCC‑ATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGAC‑TTC‑GGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACg                         <  1:462068/151‑1 (MQ=255)
                                                                                                                  cAGTATGGCGGTCCC‑ATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGAC‑TTC‑GGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGAcc                     <  2:247319/151‑1 (MQ=255)
                                                                                                                   aGTATGGCGGTCCC‑ATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGAC‑TTC‑GGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCt                    <  2:632221/151‑1 (MQ=255)
                                                                                                                    gTATGGCGGTCCC‑ATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGAC‑TTC‑GGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCtt                   <  1:105020/151‑1 (MQ=255)
                                                                                                                    gTATGGCGGTCCC‑ATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGAC‑TTC‑GGGCAACCCGGGAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCtt                   >  2:229458/1‑151 (MQ=255)
                                                                                                                         gCGGTCCC‑ATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGAC‑TTC‑GGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCt              <  2:255245/151‑1 (MQ=255)
                                                                                                                          cGGTCCC‑ATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTCCGGGCGGGATGAC‑TTC‑GGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCt              >  1:303265/1‑150 (MQ=255)
                                                                                                                            gTCCC‑ATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGAC‑TTC‑GGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcc           >  2:911075/1‑151 (MQ=255)
                                                                                                                             tCCC‑ATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCCCGGCGTGCCCGGGCGGGACGAC‑TTC‑CGGCAACCCGGCAAATCCGCCGCTGGGGGCGCAGGCGCTGCCCGGCGAGACGGACCTTGCGCTgccc          >  1:15133/1‑151 (MQ=255)
                                                                                                                                       ccAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGAC‑TTC‑GGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCt  <  2:694215/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
CTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCC‑ATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGAC‑TTC‑GGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCT  >  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: ‑