Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,510,3950CA40.0% 89.8 / 32.8 56S131R (AGC→AGA) ydcRputative DNA‑binding transcriptional regulator and putative aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (28/5);  new base A (0/22);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 6.21e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAGGCATTA  >  NC_000913/1510249‑1510542
                                                                                                                                                     |                                                                                                                                                    
cGAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCaaaa                                                                                                                                                  <  1:99567/151‑1 (MQ=255)
 gagaTGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑c                                                                                                                                                >  2:469308/1‑151 (MQ=255)
 gagaTGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑c                                                                                                                                                <  2:953835/151‑1 (MQ=255)
  agaTGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑cc                                                                                                                                               >  1:707828/1‑151 (MQ=255)
    aTGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGc                                                                                                                                             <  1:305398/151‑1 (MQ=255)
       aaGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCcac                                                                                                                                          >  1:680238/1‑151 (MQ=255)
         gCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCcac                                                                                                                                          >  2:757277/1‑149 (MQ=255)
           aGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCg                                                                                                                                      >  1:38617/1‑151 (MQ=255)
                  atCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCg                                                                                                                               >  1:740570/1‑151 (MQ=255)
                    cAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTg                                                                                                                             >  2:855929/1‑151 (MQ=255)
                    cAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTC‑CCCGA‑CCCGCGACTTTCCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGAAAAG‑CCGCCACAGCGATGAGCg                                                                                                                               <  1:137881/149‑1 (MQ=255)
                        acTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTg                                                                                                                         >  2:130434/1‑151 (MQ=255)
                          ttttttttttGATTTGTTGCGGCCCAGTCGCGACCCGTGGGTC‑TTCCCTTTTCCCCCGCCTTTCCCCGC‑CCCGCTATTTTCCCCCCCCC‑AAAAAAAAACCGGCTCGCGGGCGGGGGTAAAAAAAA‑CCGCCACAGCGATGAGCGt                                                                                                                              <  1:601257/139‑1 (MQ=255)
                          ttATTTTTTTGATAGGTGGCAGCCCGGCCGCGCCCCGTGTGTCG‑TCCGTTTTCCCCCGCCTTTCCCCGC‑CCCGCGATTTCCCCCCCCCC‑AACAACAAACCCGCTCGCTGGCGAAGGAAAGAAAAA‑CCGCCACAGCGATGAg                                                                                                                                 <  1:886062/141‑1 (MQ=255)
                          ttATATTTTTAATATGTTGCAGGCCAGTCGCGATCCGTCGGTCTTTCCTTTTCCCCCGGCCTTC‑CCCGC‑CCCGCGATTTTCCCCCCCCC‑AACAACAAACCCGCTCGCTGGCGCAGGTAAGAAAAA‑CCGCCACAGCGATGAGCGTGATTg                                                                                                                         <  2:247929/149‑1 (MQ=255)
                             tttttttGATTTGTTGCGGGCCGGTCGCGACCCGTCGGTCTTTCCTTTTGCCTCGGCTTTC‑CCCGC‑CCCGCGCTTTTCCCCCCCCC‑AACAACAAACCCGCTCGTTGGCGCAGGAAAGAAAAA‑CCGCCACAGCGATGAGCGTGATTg                                                                                                                         <  1:659192/144‑1 (MQ=255)
                               tttttGATTTGTGGCAGGCCAGGCGCGCCCCGCGGTTCGTTCGGTTTCCCGCGCCTTTC‑CCCGA‑CCCGCG‑CTTTCCCCCCCCCAAACAACAAACCCGCTCGCTGGCGCAGGAAACCAAAT‑CCGCCACAGCGATGAGCGTGATTGaaaa                                                                                                                     <  2:613704/148‑1 (MQ=255)
                               ttttgGATTTGTGGCGGGCCAGGCGCGTCCCGTCGTTCGTCCCGTTGCCCCCGCCCTTC‑CCCGC‑CCCGCGACTTTCCCCCCTCC‑AACAACAAACCCGCTCGCGGGCGCGGGAAAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGaaaa                                                                                                                     <  2:382237/148‑1 (MQ=255)
                                ttttGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACtt                                                                                                                  >  2:943703/1‑150 (MQ=255)
                                ttttGATATGTTGCAGCCCGGCCGCGATCCGGGGGTCTTTCTTTTTCCCCGGCCCTTC‑CCCGC‑CCCGCCTCTTTCCCCCCCCC‑AAAAAAAAACCCCCCCGCTGGCGCAGGTAAGAAAAC‑CCGCCACAGCGATGAGCGTGATTGAAAACtt                                                                                                                  <  1:316134/150‑1 (MQ=255)
                                       atGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGcc                                                                                                            >  2:782535/1‑149 (MQ=255)
                                        tGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAgg                                                                                                         >  2:822317/1‑151 (MQ=255)
                                        tGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGAGTGATTGAAAACTTACCGCCAgg                                                                                                         >  1:850575/1‑151 (MQ=255)
                                           tGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGaa                                                                                                       >  1:263547/1‑150 (MQ=255)
                                              aGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGc                                                                                                   >  2:678032/1‑151 (MQ=255)
                                               ggccggccGCGACCCGCGGGTC‑TTCCCTTTTCCCCCGCCTTTCCCCGC‑CCCGCGTTTTTCCCCCCCCA‑AAAAACAAACCCGCCGGCGGGCGAAGGAAAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACg                                                                                                    <  2:21365/149‑1 (MQ=255)
                                                 ccacccgcgccccGTCGTTCGTCCCTTTGCCCTCGGCTTTC‑CCCGC‑CCCGCGTTTTTCCCCCCCCC‑ACCAACAAACCCGCTCGCTGGCGGAGGAAAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGc                                                                                                   <  2:241028/137‑1 (MQ=255)
                                                   aGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAAc                                                                                              >  1:680196/1‑151 (MQ=255)
                                                       gcgATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCg                                                                                          >  1:715351/1‑151 (MQ=255)
                                                       gccacccGTCGGTCTTTCGTTTGGCCTCGGCTTTC‑CCCCC‑CCCCCGCTTTTCCCCCCCCC‑AACAACAAACCCCCTCGCGGGCGCGGGAAAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGc                                                                                           <  1:830892/145‑1 (MQ=255)
                                                        gcgccccttggttcttCCCTTTGCCCCCGCCTTTCCCCGC‑CCCGCGTTTTTCCCCCCCCC‑AACAACAAACCCGCTCGCGGGCGAAGGAAAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACt                                                                                             <  2:409808/134‑1 (MQ=255)
                                                         gATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTc                                                                                        >  1:219915/1‑151 (MQ=255)
                                                          gtcccttcgttcttCCCTTTGCCCCCGCCTTT‑CCCCACCCCGCGCTTTTCCCCCCCCC‑AACAACAAACCCCCCCGCGGGCGCAGGAAAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACt                                                                                             <  1:761777/134‑1 (MQ=255)
                                                            ccGTGGGCCTTTCCGTTT‑CCCCCGCCTTT‑CCCCG‑CCCCCCGCTTTTCCCCCCCCAAACAACAAACCCCCCCGGGGGCGCAGGAAAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTg                                                                                            <  1:993116/144‑1 (MQ=255)
                                                              gtcggtcgTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCt                                                                                   >  1:743935/1‑151 (MQ=255)
                                                                 ggtcgTTCCGTTTGCCTCGGCCTTT‑CCTGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTAtt                                                                                >  1:761692/1‑151 (MQ=255)
                                                                         cctttcgccgcgCCTTTCCCCGC‑CCCGCGCCTTTCCCCCCCCA‑AAAAACAAACCCGCTCGCGGGCGCAGGAAAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTa                                                                                  <  2:671055/131‑1 (MQ=255)
                                                                                   gcccTTC‑CCCGC‑CCCGCG‑CTTTTCCCCCCCCAAAAAACAAACCCGCCCGGGGGCGAGGAAAAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCt                                                                <  2:97755/147‑1 (MQ=255)
                                                                                      ctttccccgc‑cccgcGTTTTTCCCCCCCCC‑AAAAACAAACCCGCTCGCGGGCGGAGGAAAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGcc                                                                 <  1:831445/141‑1 (MQ=255)
                                                                                      ctttccccgc‑cccgcGCTTTTCCCCCCCCC‑AACAACAAACCCCCTCGCGGGCGCAGAAAAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCtt                                                               <  1:757156/143‑1 (MQ=255)
                                                                                        ttccccgc‑cccgcGTTTTTCCCCCCCCA‑AAAAACAAACCCGCTCGCGGGCGGGGGAAAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACa                                                            <  2:606687/146‑1 (MQ=255)
                                                                                                 cccGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCAcc                                                  >  2:368598/1‑151 (MQ=255)
                                                                                                     cgACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCAttt                                              >  1:245437/1‑151 (MQ=255)
                                                                                                        cTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTtctc                                           >  1:326265/1‑151 (MQ=255)
                                                                                                         tttcCCCCCCCA‑AACAACAAACCCCCCCGCGGGGGCGGGTAAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCAt                                                <  2:126953/145‑1 (MQ=255)
                                                                                                                ccccc‑AACAACAAACCCCCCCGCTGCCG‑AGGGAAACAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGaa                                   <  2:931363/148‑1 (MQ=255)
                                                                                                                ccTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCtgatga                                    >  1:175967/1‑150 (MQ=255)
                                                                                                                    c‑aacaaCTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTg                               >  1:439043/1‑151 (MQ=255)
                                                                                                                             aacccccccgctGGCGAAGAAAAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACt                       <  1:782414/143‑1 (MQ=255)
                                                                                                                              aaCCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTg                      >  1:973437/1‑151 (MQ=255)
                                                                                                                                 ccctcgggggCGAAGGAAAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCg                   <  2:570332/143‑1 (MQ=255)
                                                                                                                                     ccctGCCGCAGGTCAATAAAT‑CCGCCACAGCGATGAGCGTGATTGAAATCTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGc               <  1:382226/149‑1 (MQ=255)
                                                                                                                                         ggCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGtt            >  1:303650/1‑150 (MQ=255)
                                                                                                                                         ggCGAAGGTAAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGtt            <  2:1069621/150‑1 (MQ=255)
                                                                                                                                            gcAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTag          >  2:144598/1‑149 (MQ=255)
                                                                                                                                                 aaaacAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAGGCAtt   <  2:884168/147‑1 (MQ=255)
                                                                                                                                                   aGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAGGCATTa  >  2:377316/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAGGCATTA  >  NC_000913/1510249‑1510542

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