Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,332,7910AC40.0% 103.9 / 15.9 60V161G (GTG→GGG) basSsensory histidine kinase in two‑component regulatory system with BasR
Reads supporting (aligned to +/- strand):  ref base A (36/0);  new base C (0/24);  total (36/24)
Fisher's exact test for biased strand distribution p-value = 2.77e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGC  >  NC_000913/4332654‑4332940
                                                                                                                                         |                                                                                                                                                     
tACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTgg                                                                                                                                          >  1:787390/1‑151 (MQ=255)
  ccGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTGTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCg                                                                                                                                        >  1:774134/1‑151 (MQ=255)
   cGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCg                                                                                                                                        >  2:593506/1‑150 (MQ=255)
   cGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGt                                                                                                                                       >  1:43414/1‑151 (MQ=255)
     gaagaaAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTc                                                                                                                                     >  2:877963/1‑151 (MQ=255)
        gaaAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGc                                                                                                                                   >  1:851665/1‑150 (MQ=255)
           aaTGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCCTCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGtt                                                                                                                               >  2:111760/1‑151 (MQ=255)
           aaTGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGTCGTGCCACTAACGGTGCTACATCAATGTGATGCGTATTCGCCAGCAGTTGCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAg                                                                                                                                 >  2:1042737/1‑149 (MQ=255)
             tgccTGCCCGGCACGCGCCGGTGGCAGCGGTTGGGAGCCCCTCCCCATCATCTGAACAGGCCGTCCCACAAGCGTTGCTACATAAAGGTATTGCGTTTCCCCCAGCTGTCCCGGATGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTc                                                                                                                              <  2:1168936/147‑1 (MQ=255)
             taactacccgccccgcgcgATTTACAGGAGGTGGGAGAGCCTCTCCACCTTCTAACCAACCCGTGCCCCTATGGGGGTTACAACAAGGTGGGTCTTTTCCCCCAGCTTTTCCAGATCCAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTc                                                                                                                              <  1:218025/137‑1 (MQ=255)
                  gTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTgcgc                                                                                                                        >  2:306676/1‑151 (MQ=255)
                   tCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGcgcg                                                                                                                       >  2:201149/1‑151 (MQ=255)
                    ccgccccgcgCCGTTGGCAGCAGCGGGGAGCCCCTCCCCTTCATTTGATCAAGCCGGGCCACTACGGGTGCTACATAAATGTGATGCGTTTTCGCCAGCAGTCCCAGATGCAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTg                                                                                                                           <  1:357289/140‑1 (MQ=255)
                         ccgcgCCGGTTGCGGCGGCGGGGAGCCCCCCTCCACCTCCGAACCAAGCCGCCCCACTACCGGTGCTACATCAAGGTGGGGCTTTTCCCCCAGCAGTCCCAGATGAAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTc                                                                                                                  <  2:760396/149‑1 (MQ=255)
                           gcgcCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGc                                                                                                               >  1:1007092/1‑151 (MQ=255)
                            cgcCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAg                                                                                                                >  2:183389/1‑149 (MQ=255)
                              ccAGTGGCGGCGGCGGGGAGCCCCTCCCCATCTTCTAATCAGCCCGTGCCACAACGGGTGCTACATCAATGTGATGCGTTTTCGCAAGCAGTCCCAGATGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGt                                                                                                            <  2:1064134/151‑1 (MQ=255)
                                aGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTa                                                                                                           >  2:277010/1‑150 (MQ=255)
                                aGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTa                                                                                                           >  2:92195/1‑150 (MQ=255)
                                  ttGCACCAGCGGGGAACCCCCCCCCTTCTTCTGATCAGGCGGCGCCACTAACGGTCCAACATCAATGTGAGGCGTTTTCACCAGCAGTTCCAGAGGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTaaa                                                                                                         <  1:770206/150‑1 (MQ=255)
                                   tGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAaac                                                                                                        >  1:303295/1‑150 (MQ=255)
                                   tGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAaac                                                                                                        >  1:897006/1‑150 (MQ=255)
                                            tGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGt                                                                                              >  1:578687/1‑151 (MQ=255)
                                             gggAGCCCCTCCCCATCATCTGATCAAGCCTTCCCACTACCGTTGCTACATCAATGTGATGCTTTTCCGCAAGCATTTCCAGATGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGt                                                                                              <  1:637976/150‑1 (MQ=255)
                                             gggAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGtt                                                                                             >  1:378586/1‑151 (MQ=255)
                                              ggaccccctccCCTTCATCTACACAACCCGTCCCCCAACCGGTGCTACATCAATGTGTTCCTTTTTCGCCAGCATTCCCGGATGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTa                                                                                            <  2:461865/144‑1 (MQ=255)
                                                   ccccccccaccaTCGGAACAAGCGGGCCCACAAAGGGGGCTAAAAAAATGTGATGTTTTTCCGCCGGCTTTCCCAGAGGCAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATcc                                                                                         <  2:476925/143‑1 (MQ=255)
                                                              atcGGACCACGCCGTCCCACTAACGGTGCAACATCAATGTGATGCGTTTTCGCCGGCTGTTCCAGAAGAAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTc                                                                             <  2:1056321/150‑1 (MQ=255)
                                                               tcTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTc                                                                             >  2:1155712/1‑149 (MQ=255)
                                                                  gATCAAGCCGTCCCACTAACGGTGCAACATAAATGTGATGCGTTTTCGCCAGCATTCCAAGATGCAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGt                                                                                    <  1:1122477/139‑1 (MQ=255)
                                                                   aTCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCg                                                                        >  1:1005697/1‑150 (MQ=255)
                                                                   aTCAAGCCGTCCCACAAGGGGTGCAAAATAAATGTGATGCTTTTTCGCCAGAGGTTCCAGATAAAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAg                                                                           <  1:664791/147‑1 (MQ=255)
                                                                      agcccgccccACTACCGGTGCTACATAATTGTGTTGGTTTTCCGCCGGCAGTCCCAGATAAAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGc                                                                         <  2:796029/138‑1 (MQ=255)
                                                                      aaGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTg                                                                    >  2:100886/1‑151 (MQ=255)
                                                                      aaGCCGTGCCACTAACGGGGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGCCAGGCGAc                                                                      >  2:337383/1‑149 (MQ=255)
                                                                      aaGCCGTCCCACAAACGGTGCTACATCAATGTGATGCTTTTCCGCCAGCATTCCCAGATGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGc                                                                         <  1:61515/146‑1 (MQ=255)
                                                                         ccgccccACTAACGGTCATCAAACAGTGTGTTGTTTTTTCCCCAGAATCCCCAGATGCAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACt                                                                     <  2:742146/142‑1 (MQ=255)
                                                                              caccaaaCGTTGCTACAAAATTGGGATGGTTTTCCGCAAGCAGTTCCAGATGAAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGa                                                              <  2:29212/144‑1 (MQ=255)
                                                                                   aaggttGCTACATAATTGTGATGCGTTTGCCCCGGCATTCCCAGAGGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGttt                                                        <  2:87226/145‑1 (MQ=255)
                                                                                     cGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAAc                                                     >  2:685937/1‑151 (MQ=255)
                                                                                          cTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCg                                                 >  2:1006629/1‑150 (MQ=255)
                                                                                             cATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAac                                             >  2:762992/1‑151 (MQ=255)
                                                                                                cAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTAGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAaccac                                          >  2:96550/1‑151 (MQ=255)
                                                                                                 aaTGTGATGCTTTTCCGCCAGCAGTTCAGGATGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAacc                                            <  1:384387/148‑1 (MQ=255)
                                                                                                  aTGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCg                                        >  1:121279/1‑151 (MQ=255)
                                                                                                     tgATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCtt                                     >  2:773226/1‑151 (MQ=255)
                                                                                                      gATGCGTTTCCGCCAGCAGTTCCAGAGGCAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGc                                       <  1:871503/148‑1 (MQ=255)
                                                                                                        tgttttttCCCCCGCAGCCCCAAAGACAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCtt                                     <  2:1123128/144‑1 (MQ=255)
                                                                                                                cGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAgg                          >  2:1122376/1‑151 (MQ=255)
                                                                                                                cGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAgg                          >  2:344605/1‑151 (MQ=255)
                                                                                                                cGCCAGCAGTCCGAGAACAAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAgg                          <  2:138310/151‑1 (MQ=255)
                                                                                                                 gCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAggg                         >  1:310278/1‑151 (MQ=255)
                                                                                                                       agTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGgcgc                   >  1:11722/1‑151 (MQ=255)
                                                                                                                       agTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGgcgc                   >  2:902590/1‑151 (MQ=255)
                                                                                                                         tccgaGATCCAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGgc                     <  1:404035/143‑1 (MQ=255)
                                                                                                                           ccAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTCAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtgtg               >  2:227752/1‑151 (MQ=255)
                                                                                                                           ccAGATGAAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtg                 <  1:490770/149‑1 (MQ=255)
                                                                                                                           caaataccAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTg                       <  2:832350/136‑1 (MQ=255)
                                                                                                                               aTGCAAACGAACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGaa             >  2:3874/1‑149 (MQ=255)
                                                                                                                                        cACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGc  >  2:297242/1‑151 (MQ=255)
                                                                                                                                         |                                                                                                                                                     
TACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGC  >  NC_000913/4332654‑4332940

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: