Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,622,8160GA44.3% 32.3 / 61.6 62intergenic (+75/+144)yneM/dgcZinner membrane‑associated protein/diguanylate cyclase, zinc‑sensing
Reads supporting (aligned to +/- strand):  ref base G (34/0);  new base A (0/27);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 6.33e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAAC  >  NC_000913/1622668‑1622954
                                                                                                                                                       |                                                                                                                                          
ttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGa                                                                                                                                           >  1:1017752/1‑150 (MQ=255)
     gCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaaaaaa                                                                                                                                      >  1:1154726/1‑150 (MQ=255)
       cGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAAc                                                                                                                                   >  2:931167/1‑151 (MQ=255)
       cGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAAc                                                                                                                                   >  2:853574/1‑151 (MQ=255)
             gggAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGcc                                                                                                                               >  1:455940/1‑149 (MQ=255)
                      tttattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcggc                                                                                                                    >  1:989775/1‑151 (MQ=255)
                         attttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcggcgg                                                                                                                  >  1:1109413/1‑150 (MQ=255)
                          ttttttgtttttctgcCCGCTTATTCCCCCAACAAAGGGGATGACAAATAAACGGAGATACCCCCTCACCTACCCGCCCCCTTTTAACGGTGGTGAACAGGGGGCTTTTTTTTTATGACGACAAAAAAAAACCGCCAGTAAACcggcgg                                                                                                                  <  1:925437/149‑1 (MQ=255)
                               tggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGc                                                                                                           >  1:85667/1‑151 (MQ=255)
                               tgggtttctggCCGCGCATTTCGGCTACAACTGGGCTGACTAATGAACGGAGATAATCCCTCATATAACCGGCACCTTGTTACTGCTGTGTAGAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGACTg                                                                                                            >  1:974403/1‑150 (MQ=255)
                                ggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGCGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCCGCGgagactttt                                                                                                          >  1:615889/1‑146 (MQ=255)
                                   tttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc                                                                                                       >  1:302945/1‑151 (MQ=255)
                                   tttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc                                                                                                       >  2:610450/1‑151 (MQ=255)
                                    ttctggCCGCTTATTTCAGCAAAAAAGGGGATAATAAAGGAAGGGAAATAACCCCCCCCCTACCCGCCCCTTTTTTACATTTGTGTAAAGGGGGCTTATTTTTTAGGGGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc                                                                                                       <  2:351407/150‑1 (MQ=255)
                                        ggccgtttttttCACCCAAAAGGGGGATAACAAATGAGGGGAAAAAACCCCCCCCCAACCCGCCCCTTTGTTACAGTTGTGTAAAGGGGGCCTGTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCa                                                                                                      <  2:139201/138‑1 (MQ=255)
                                        gccggctttttcccgCCAAAAAGGGGATGATAAAGAAAGGGAAATAACCCCCCCCCAACCCGGCCCTTTTTTACGGTTGTGAAAAGGGGGCTTTTTTTTTAGGCGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCAt                                                                                                     <  1:154331/134‑1 (MQ=255)
                                          ccGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGa                                                                                                  >  2:450685/1‑149 (MQ=255)
                                           cgcgTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTCATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt                                                                                                 >  1:276204/1‑149 (MQ=255)
                                            gcgTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt                                                                                                 >  1:161732/1‑148 (MQ=255)
                                              tttttccacccAAAAATGGTATGATAAATAAAGGAAAATACCCCCCCCCCAACCCGCCCCTTTTTTACAGTTGTGAAAAGGGGGCTTGTTTTTTAGGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCa                                                                                                      <  1:176294/132‑1 (MQ=255)
                                               tttttCAGCCAAAAGGGGGATGACTAATGAGGGGAGATACCCCCCCCCCTACCCGCCCCTTTTTTACAGTTGTGTACAGGGGGCTTTTTTTTTAGGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAtt                                                                                            <  1:68712/148‑1 (MQ=255)
                                                 tttCACCCAAAAAGGGGATGACAAATGAAGGAAGATACCCCCCCCCCAACCCGCCCCTTTTTTACATTTGTGAAAAGGGGCCTTTTTTTTTAGGGCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAg                                                                                               <  2:631995/145‑1 (MQ=255)
                                                  ttCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTtgt                                                                                         >  2:569836/1‑150 (MQ=255)
                                                     aGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGttt                                                                                     >  1:879543/1‑151 (MQ=255)
                                                     aGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAATTGTGttt                                                                                     >  1:436597/1‑151 (MQ=255)
                                                      gCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGtttt                                                                                    >  1:912871/1‑151 (MQ=255)
                                                         aaaaaGGGGATGAAAAAGGAAGGGAAATAACCCCCCCCCTCCCCGCCCCTTTTTTACATTTGTGTAAAGGGGGCTTTTTTTTTATGACGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGtttt                                                                                    <  2:557724/146‑1 (MQ=255)
                                                           aaaTGGGATGACTAATGAACGGAGATAATACCTCACCTAACCGGCCCCTTGTTACAGTTGTGCACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGAGGAGAATGATTGCATGGATAGATTTGTGATTTGCtt                                                                                >  1:1083116/1‑150 (MQ=255)
                                                           aaaGGGGATGATAAATGAAGGGAAAAACCCCCTCCCCAACCCGCCCCTTTTTTACAGTGGTTAAAAGGGGCCTTTTTTTTTATGCGGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtt                                                                                <  1:3543/150‑1 (MQ=255)
                                                             aTGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGGTGATGCTTGCATGGGTAGATTTTTGTTTTTCtttt                                                                              >  1:780176/1‑150 (MQ=255)
                                                                         aTGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAGAGATTTGTGTTTTGCTTTTACGCCAACAGGc                                                                  >  1:530438/1‑150 (MQ=255)
                                                                           gAACGGAGATAACCCCCCCCAAACCCGCCCCTTTTTCAGGGTGTATCGCAGGGGTGCTGTTTTTAATGACGCAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAgg                                                                   <  2:234801/147‑1 (MQ=255)
                                                                            aggggAAATAACCCCTCACCAACCCGCCCCTTTTTTCGAGTTGTGTAAAGGGGGCTTGTTTTTTAGGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAt                                                                <  1:938565/146‑1 (MQ=255)
                                                                             aCGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAtttt                                                             >  2:1011332/1‑151 (MQ=255)
                                                                                 aaaaacccccTCCCCTACCCGCCCCTTTTTTACGGTTGTGTAAAGGGGGCTTTTTTTTTATGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTc                                                            <  1:704510/141‑1 (MQ=255)
                                                                                     acccccTCACAAACCCGCCCCTTTGTTACGGTTGTGAACAGGGGGCTTTTTTTTTAGGCGGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCAc                                                      <  2:638596/147‑1 (MQ=255)
                                                                                        cccccaccaccccGCCCCTTTGTAACTGTTGTATAAAGGGGTCTTTTTTTTTATGGGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTg                                                    <  2:1037986/138‑1 (MQ=255)
                                                                                          cTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATa                                                 >  2:483876/1‑150 (MQ=255)
                                                                                            cACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGTCGATAAAAAACCGCCAGTAAACCGGCGGCTAATGCCTGCATGGATAGATTTGTGTTTTGCTTTTATGCTAACAGGCATTTTCCTGCACTGATaa                                                >  1:924160/1‑149 (MQ=255)
                                                                                              ccTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGa                                             >  1:1010631/1‑150 (MQ=255)
                                                                                               cTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAAt                                           >  2:533450/1‑151 (MQ=255)
                                                                                                 acccGCCCCTTTTTTACATTTGTGTAAAGGGGGTTTTTTTTTTATGGCGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACg                                              <  1:857023/144‑1 (MQ=255)
                                                                                                 acccGCCCCTTTTTTACAGTTGTGTAAAGGGGCCTTGTTTTTTATGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGa                                             <  2:122527/145‑1 (MQ=255)
                                                                                                 aaCCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCg                                         >  2:340550/1‑151 (MQ=255)
                                                                                                    cgcccctTTTTTACAGTTGTAAAAAGGGGCCTTGTTTTTTATGACGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATaa                                                <  1:470072/138‑1 (MQ=255)
                                                                                                       ccctTTTTTACGGTTGTGAAAAGGGGGCCTTTTTTTTATGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCg                                         <  1:397734/145‑1 (MQ=255)
                                                                                                        cccTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGacac                                  >  2:203307/1‑151 (MQ=255)
                                                                                                         ccTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGacac                                  >  1:103667/1‑150 (MQ=255)
                                                                                                              ttACGGTTGTGTACAGGGGGCTTTTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAg                             <  2:989796/150‑1 (MQ=255)
                                                                                                              taaCATTTGTGAACAGGGGGGTTTTTTTTTTGGGCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTa                              <  1:551009/147‑1 (MQ=255)
                                                                                                                     ttGTGTACAGGGGCCCTATTTTTTAGGGGGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCa                        <  1:690565/149‑1 (MQ=255)
                                                                                                                     ttGTGTAAAGGGGGCTTGTTTTTTATGACGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAg                       <  2:359776/150‑1 (MQ=255)
                                                                                                                      tgtgTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGtt                     >  1:48535/1‑151 (MQ=255)
                                                                                                                      tggtaaaCAGGGGGCCTGTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCa                        <  1:170442/143‑1 (MQ=255)
                                                                                                                            aaagggggCTTGTTTTTTAGGCGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtct                 <  2:977240/146‑1 (MQ=255)
                                                                                                                             cAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtctc                >  2:181182/1‑150 (MQ=255)
                                                                                                                                tggaccTGTTTTTTATAGGGGCTAAAAAAATCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACATTAGCATCAGtttt                   <  1:796616/140‑1 (MQ=255)
                                                                                                                                 gggCTTGTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatga           <  2:1044244/151‑1 (MQ=255)
                                                                                                                                   gCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaat         >  2:765196/1‑151 (MQ=255)
                                                                                                                                       tGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGtt      >  2:1136171/1‑151 (MQ=255)
                                                                                                                                          tttttAGGACGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTaa    <  1:1163109/150‑1 (MQ=255)
                                                                                                                                            tttATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAAc  <  2:630437/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                          
TTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAAC  >  NC_000913/1622668‑1622954

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