Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,671,4150CT41.4% 81.7 / 29.6 58K160K (AAG→AAA) yhjBputative DNA‑binding transcriptional response regulator
Reads supporting (aligned to +/- strand):  ref base C (30/4);  new base T (0/24);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 1.60e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCTGCTGGCTGTAAAATTATCCTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATT‑CCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTC‑AGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATCCCCATCGCCGCACTATTTACCGCC  >  NC_000913/3671266‑3671559
                                                                                                                                                        |                                                                                                                                                
tCCTGCTGGCTGTAAAATTATCCTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATT‑CCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTC‑AGCGCTCTGCCGATCTCCt                                                                                                                                                 >  2:11029/1‑151 (MQ=255)
    gctggctgTAAAATTATCCTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATT‑CCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTC‑AGCGCTCTGCCGATCTCCTTAt                                                                                                                                              <  2:427808/150‑1 (MQ=255)
    gctggctgTAAAATTATCCTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATT‑CCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTC‑AGCGCTCTGCCGATCTCCTTAt                                                                                                                                              <  1:716989/150‑1 (MQ=255)
      tggctgTAAAATTATCCTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATT‑CCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTC‑AGCGCTCTGCCGATCTCCTTATTTg                                                                                                                                           >  2:1086532/1‑151 (MQ=255)
       ggctgTAAAATTATCCTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATT‑CCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTC‑AGCGCTCTGCCGATCTCCTTAttt                                                                                                                                            <  1:388969/149‑1 (MQ=255)
       ggctgTAAAATTATCCTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATT‑CCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTC‑AGCGCTCTGCCGATCTCCTTATTTg                                                                                                                                           >  1:1170268/1‑150 (MQ=255)
           gTAAAATTATCCTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATT‑CCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTC‑AGCGCTCTGCCGATCTCCTTATTTGATTcg                                                                                                                                      >  2:598947/1‑151 (MQ=255)
             aaaaTTATCCTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATT‑CCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTC‑AGCGCTCTGCCGATCTCCTTATTTGATTcgcc                                                                                                                                    >  1:754962/1‑151 (MQ=255)
             aaaaTTATCCTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATT‑CCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTC‑AGCGCTCTGCCGATCTCCTTATTTGATTcgcc                                                                                                                                    >  2:811624/1‑151 (MQ=255)
                  tATCCTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATT‑CCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTC‑AGCGCTCTGCCGATCTCCTTATTTGATTCGCccgccg                                                                                                                               >  2:982293/1‑151 (MQ=255)
                    tCCTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATT‑CCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTC‑AGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCt                                                                                                                             >  2:269701/1‑151 (MQ=255)
                     ccTCTTCAGGAGGGGATGTTTAACATCCTTGCGGCCTGGGAGCGATTTTTCACTTCCAGACGGCGATACAGGGATT‑CCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTC‑GGCGCTCTGCCGCGCTCCTTACTTGATTcgccc                                                                                                                                   >  2:700224/1‑144 (MQ=255)
                      ctctTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATT‑CCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTC‑AGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTa                                                                                                                            >  2:138359/1‑150 (MQ=255)
                        ctTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATT‑CCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTC‑AGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAAc                                                                                                                          >  1:334757/1‑150 (MQ=255)
                              gaagatatatTTAACATCTTTGCCGCTTGGTGGTTTTTTTTCACTCCCAGACGGCGAAACAGGGCTT‑CAAGATGCTTTTTTCCTGTTCCGG‑TACTGATATTC‑AGCGCTCTGCCGATCTCTTTTTTTGATTCGCCCGCCGCTAACATGGtt                                                                                                                    <  2:34218/144‑1 (MQ=255)
                                     tatTTACCATCTTTGCCGCCTGGGTGCGTTTTTTCCCTCCCAGACGGCGATACGGGGATT‑CCAGATGCGTTTTTCCCTTCCGGG‑TACTGATATTC‑GTCGCTCTGCCGATCTCTTTATTTTTTTCGCCCGCCGCTAACATGGTTaaaa                                                                                                                <  1:1071701/147‑1 (MQ=255)
                                            catcatTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATT‑CCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTC‑AGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTgg                                                                                                     >  2:502074/1‑151 (MQ=255)
                                                 ttGCCCCCTGGGTGCTTTTTTCCACTCCCAGAGGGCGATACGGGGATC‑CCAGATGCGTTTTTCCCTTCCCGG‑TACTGATATAC‑AGCGCTCTCCCGATCTCTTTTTTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTggcgg                                                                                                  <  1:141149/149‑1 (MQ=255)
                                                   gccgctTGGGGGCGTTTTTCCCCTTCCAGCCGGCGATAGAGGGATC‑CCAGTTGCTTTTTTCCCTTTCCGT‑TATTAATATTC‑AGCGTTCTCCCAATCTCTTTTTTTGTTCCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTggcgg                                                                                                  <  2:870990/147‑1 (MQ=255)
                                                      gccTGGGTGCTTTTTTTCCCTCCCGGAGGGCGATACGGGAATT‑CCAGATGCGTTTTTCCCGTTCCGT‑TACTATTATCC‑AGCGCTCTCCCATTCTCTTTTTTTGTTTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTa                                                                                           <  1:261095/151‑1 (MQ=255)
                                                        cTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATT‑CCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTC‑AGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAAc                                                                                         >  2:915146/1‑151 (MQ=255)
                                                          gggTGCTTTTTTTCCCTTCCAGACGGCGATACAGGGATT‑CCAGTTGCTTTTTTCCCTTTCCGT‑TACTAATATCC‑AGCGTTCTGCCGATCTCTTTTTTTGATCCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGa                                                                                       <  2:203089/151‑1 (MQ=255)
                                                          gggTGCGTTTTTCCCCTCCCAGAGGGCGATACAGGGATC‑CCAGATGCGTTTTTACCTTTCCGT‑TACTGATTTCC‑AGCGCTCTGCAGTTCTCTTTTTTTGATCCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTggcgggcg                                                                                               <  1:65223/143‑1 (MQ=255)
                                                           ggtgttttttttCCCTCCCAGGGGGGGATAGGGGGATT‑CAAGTTGC‑TTTTTTCCCTTCCGGTTATTGATATTC‑GGCGCTCTGCCGATCTCTTTTTTTGTTTCCCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGAt                                                                                      <  2:51311/144‑1 (MQ=255)
                                                           ggTGCGTTTTTTCACTCCCAGACGGCAATACAGGGATT‑CCAGATGCGTTTTTCCCGTCCCGG‑TACTGATATTC‑AGCGCTCTGCCGTTCTCTTTATTTTATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGAt                                                                                      <  1:910369/151‑1 (MQ=255)
                                                                  tttttCACTTCCAGACGGCGATACAGGGATT‑CCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTC‑AGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTg                                                                                   >  2:990653/1‑147 (MQ=255)
                                                                  ttttccccccccAGCCGGAGATACAGGGCTT‑CTAGATCCGCTTTTCCCTTACTGT‑TTCTGATTGTC‑AGCGCCCTCCCGTTCTCTTTATTGGTTTCTCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTgaga                                                                                <  1:461682/140‑1 (MQ=255)
                                                                       ccctcccAGACGGCAATACAGGTTTC‑CCAGTTGC‑TTTTTTCCCCTCCGG‑TTATTATTTTCTGGCGCTCTCCCTTCCTCTTTTTTTGTTTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCttt                                                                           <  2:239704/145‑1 (MQ=255)
                                                                        aCTTCCAGACGGCGATACAGGGATT‑CCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTC‑AGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCttt                                                                           >  2:173016/1‑149 (MQ=255)
                                                                        aCTTCCAGACGGAGAAAGGGGGATT‑CAAGATGCCTTTTTCCCGTTCGGG‑TACTGATATAC‑AGCGTTTTGCGGATCTTTTTATTTGATCCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCttt                                                                           <  1:168154/149‑1 (MQ=255)
                                                                           cccAGACGGCGATACAGGGATT‑CCAGATGCGTTTTTCCCGTCCCGG‑TATTGATATTC‑GGCCCTCTCCCGATCTCTTTTTTTTATCCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAAt                                                                        <  1:834327/148‑1 (MQ=255)
                                                                             cAGACGGCGATACAGGGATT‑CCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTC‑AGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTcc                                                                    >  2:587730/1‑151 (MQ=255)
                                                                              aGACGGCGATACAGGGATT‑CCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTC‑AGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCt                                                                   >  2:95991/1‑151 (MQ=255)
                                                                                   gcaaTAAAGGTATT‑CCAGATGTCTTTTTCCCTTCCCGT‑TCCTGTTATCC‑AGCGCTCTCCCGATCTCTTTATTTGCTTCCCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCtt                                                                  <  2:223464/144‑1 (MQ=255)
                                                                                    caaTACGGGGATC‑CCAGATCCG‑TTTTTCCCTTCCGGTTACTGATATCC‑AGCGCTCTGCCGTTCTCTTTTTTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCg                                                             <  1:126520/149‑1 (MQ=255)
                                                                                         cAGGGATT‑CCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTC‑AGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGcgtc                                                        >  2:782847/1‑151 (MQ=255)
                                                                                         cAGGGATT‑CCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTC‑AGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGcgtc                                                        >  2:299279/1‑151 (MQ=255)
                                                                                         acgggtttcccAGATGCTTTTTTACCTTTCCGT‑TCCTGATATTC‑AGCGCTCTCCCGATCTTTTTATTTGTTCCCCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCg                                                             <  1:487835/138‑1 (MQ=255)
                                                                                           gggATT‑CCAGTTGCGTTTTTCCCGTTCCGG‑TACTGATTTTC‑AGCGCTCTGCCGTTCTCTTTTTTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGcgt                                                         <  1:131247/148‑1 (MQ=255)
                                                                                           gggATT‑CCAGATGCTTTTTTACCGTTCCGT‑TACTGATATTC‑AGCGCTCTGCCGATCTCTTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGcgt                                                         <  1:133871/148‑1 (MQ=255)
                                                                                           gggATC‑CAAGATGCTTTTTTCCCGTCCCGG‑TACTAATATTC‑AGCCCTCTGCCGATCTCTTTTTTTGATCCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGcgtcgt                                                      <  1:36997/151‑1 (MQ=255)
                                                                                            ggATTCCCAGTTGCTTTTTTCCCTTCCCGG‑TACTGATATTC‑ACCGCTCTGCTTTTCTCTTTTTTTGTTCCCCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCgg                                                            <  2:210058/145‑1 (MQ=255)
                                                                                                t‑CCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTC‑AGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGcc                                                 >  1:352582/1‑151 (MQ=255)
                                                                                                  caagttGCGTTTTTCCCTTTCCGG‑TACTAATTCCC‑AGCGCTCTGCCGATCTCTTTATTTGCTCCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGc                                                  <  2:147310/144‑1 (MQ=255)
                                                                                                   cAGATGCGCTTTTACCGTTCCGG‑TACTGATATTC‑AGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGcc                                                 >  2:899/1‑149 (MQ=255)
                                                                                                     gATGCTTTTTTCCCGTTCCGG‑TATTGATATTC‑GGCCCCCTGCCTTTCTCTTTATTTGTTTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAg                                               <  2:719865/149‑1 (MQ=255)
                                                                                                       tGCGCTTTTACCGTTCCGG‑TACTGATATTC‑AGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGtctc                                           >  1:275788/1‑151 (MQ=255)
                                                                                                            ttttCCCTTCCCGT‑TACTGATATCC‑ACCCCTCTCCCGATCTCTTTATTTGTTCCCCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCa                                         <  1:785273/148‑1 (MQ=255)
                                                                                                               tACCGTTCCGG‑TACTGATATTC‑AGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTCATGTCCTTTTCCGGCGCCGTCCGCCAGTCTCCAGGc                                      >  2:677340/1‑148 (MQ=255)
                                                                                                                 ccttccggg‑TATTAATATCC‑GGCGCTCTCCCTTTCTCTTTATTTGATCCCCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAg                                    <  1:581350/141‑1 (MQ=255)
                                                                                                                  cGTTCCGG‑TACTGATATTC‑AGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAAc                                >  2:250778/1‑151 (MQ=255)
                                                                                                                           tACTAATATCC‑AGCCTTGGGCCGTTCTCCTTTTTTGCTCCTCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAAca                               <  1:475270/144‑1 (MQ=255)
                                                                                                                                 tatTC‑AGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATCCCCATcgc                  >  1:671783/1‑151 (MQ=255)
                                                                                                                                          gcTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATCCCCATCGCCGCACTAt          >  1:198911/1‑151 (MQ=255)
                                                                                                                                              tGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATCCCCATCGCCGCACTATTTAc      >  2:1171021/1‑151 (MQ=255)
                                                                                                                                               gCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATCCCCATCGCCGCACTATTTAcc     >  1:15955/1‑151 (MQ=255)
                                                                                                                                                   aTCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATCCCCATCGCCGCACTATTTACCGcc  >  1:24326/1‑150 (MQ=255)
                                                                                                                                                   aTCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATCCCCATCGCCGCACTATTTACCGcc  >  1:727481/1‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                                
TCCTGCTGGCTGTAAAATTATCCTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATT‑CCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTC‑AGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATCCCCATCGCCGCACTATTTACCGCC  >  NC_000913/3671266‑3671559

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: ‑