Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,056,539 T→A 6.4% intergenic (‑200/+86) glnL ← / ← glnA sensory histidine kinase NtrB/glutamine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,056,5390TA 6.4% 184.7 / 4.2 65intergenic (‑200/+86)glnL/glnAsensory histidine kinase NtrB/glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base T (31/29);  new base A (2/2);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.01e-01

CGCGAAGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCAT  >  NC_000913/4056402‑4056669
                                                                                                                                         |                                                                                                                                  
cgcgAAGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCg                                                                                                                               <  2:441486/143‑1 (MQ=255)
cgcgAAGCATAAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGTCTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCg                                                                                                                               <  1:483490/143‑1 (MQ=255)
 gcgAAGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGt                                                                                                                              <  1:585703/143‑1 (MQ=255)
 gcgAAGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGt                                                                                                                              <  1:507850/143‑1 (MQ=255)
   gAAGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCg                                                                                                                            >  1:153086/1‑143 (MQ=255)
    aaGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGc                                                                                                                           >  1:477437/1‑143 (MQ=255)
     aGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCa                                                                                                                          >  1:121981/1‑143 (MQ=255)
        aTCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGCCGCACTTGCGTCGCATcc                                                                                                                       >  1:632872/1‑143 (MQ=255)
         tCAGAGAATGGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATGCGACGCAAATGCGTcttatccg                                                                                                                      <  1:348591/143‑8 (MQ=255)
           agagAATTGACGGAGAATAAAGCCCATGCAGGGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATGCGACGCAAATGCGTCTTATCCGGc                                                                                                                    <  2:684266/143‑1 (MQ=255)
             agaAGTGACGGAGAAACAAGCCCAGGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGCCTCACCTGCTCCGGCCTACCACATCAGCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGt                                                                                                                  <  2:581350/143‑1 (MQ=255)
                   ggcGGAGAAAAAAGCCCATGCAGAGAGCAGCAACCCATTAAAGAAATACCGCTCGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGAtgtg                                                                                                            <  1:180882/141‑1 (MQ=255)
                    aCGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGc                                                                                                           <  2:668822/143‑1 (MQ=255)
                      ggAGAAAAAAGCCCATGCAGAGATGGGCTACAGACAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCtt                                                                                                         <  1:624681/143‑1 (MQ=255)
                         gaAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGcc                                                                                                      <  1:301266/143‑1 (MQ=255)
                            aaaaGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCAGGCATACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGGGATGTGATTGcacac                                                                                                   >  1:199795/1‑139 (MQ=255)
                                gCCCATGCAGCGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGGCGGATAAGACGCATTTGCATCGCATCCGGCGTGATGTGCTTGCCACCGATc                                                                                               <  2:420419/143‑1 (MQ=255)
                                gCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATc                                                                                               <  2:413992/143‑1 (MQ=255)
                                  ccATGCAGAGATGGGCTACAGATAGCTGACAAACATCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCCTGTAGGCCGGCTAAGACGCATTTGCGTCACATACGGCGTGATGTGACTGCCACCggcacc                                                                                             >  2:578452/1‑138 (MQ=255)
                                    aTGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTCCCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGGGCTTGGCACCGaacgacg                                                                                           >  1:773133/1‑137 (MQ=255)
                                     tGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAgg                                                                                          <  1:707987/143‑1 (MQ=255)
                                      gCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGc                                                                                         >  1:792894/1‑143 (MQ=255)
                                      gCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGGGATGTGCTTGCCACCGATCCCAGGc                                                                                         >  1:122577/1‑143 (MQ=255)
                                       cAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGcc                                                                                        >  1:532156/1‑143 (MQ=255)
                                          agaCGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCCTCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGc                                                                                     <  1:614144/143‑1 (MQ=255)
                                                ccTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAgaga                                                                               <  1:441487/142‑1 (MQ=255)
                                                  tACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCCGCGATCCCAGGCCTGCCAGAGAcc                                                                             >  2:802946/1‑142 (MQ=255)
                                                   aCAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGCCGTGATGTGCTTGCCACCGATCCCAAGGCTGCCAGCGACAg                                                                            >  1:251169/1‑143 (MQ=255)
                                                    cAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAgg                                                                           <  1:393312/143‑1 (MQ=255)
                                                     aGATAGCTGACAAACTTCACGTTGGAGCGCGACTCACCTGCGCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGGGGTTGCCAGCGAGACGAGGCCTGCCAGAGACAggg                                                                          >  1:44707/1‑142 (MQ=255)
                                                          ggtGACAAACTTCGCGTTGGCGAGCGACTCGCCTGCTCCGGCCTACCACGTCCTCCCCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGaaaa                                                                     <  2:90315/141‑1 (MQ=255)
                                                          gCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGaaaa                                                                     <  1:753094/143‑1 (MQ=255)
                                                            tGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGGTGGGCTTGCCACCGGTCCCAGGCCTGGCAGAGGCAGGCGGACAGt                                                                   >  1:592542/1‑143 (MQ=255)
                                                            tGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGt                                                                   >  1:784903/1‑143 (MQ=255)
                                                             gACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGtt                                                                  >  1:310559/1‑143 (MQ=255)
                                                              acaacgttCACGTTGGAGAGCGCCTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGttt                                                                 <  2:684955/137‑1 (MQ=255)
                                                              aCAAACTTCCCGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGttt                                                                 <  1:647551/143‑1 (MQ=255)
                                                                     tCATGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGc                                                          >  1:386206/1‑143 (MQ=255)
                                                                     tCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGc                                                          >  2:150515/1‑143 (MQ=255)
                                                                         gTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGGTCCCAGGCCTGCCGTAGACAGGCGAACAGTTTCCACGGCAACt                                                      >  1:759945/1‑143 (MQ=255)
                                                                                  cGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCCCGGCACCCACAacaccc                                             >  2:148170/1‑141 (MQ=255)
                                                                                    aCTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCGACTAAAACGCTTAg                                           >  1:606080/1‑143 (MQ=255)
                                                                                        aCCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGCTGTCCTTGCCACCGATCCCAGGCCTCCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGc                                       >  1:412124/1‑143 (MQ=255)
                                                                                        aCCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGc                                       >  1:435303/1‑143 (MQ=255)
                                                                                         ccTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCt                                      >  1:765316/1‑143 (MQ=255)
                                                                                            gCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTgta                                   <  2:227974/143‑1 (MQ=255)
                                                                                            gCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTgta                                   <  1:115055/143‑1 (MQ=255)
                                                                                              tCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGCGACAGGCGGAAAGTTTCCCCGGCAACCAAAACACTTAGACGCTgtggg                                 >  1:450469/1‑140 (MQ=255)
                                                                                                cGGCCTAGCAGATCATCGCCGTGTAGGCTGGATAAGACGCATTTGCGTCGCAGCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCGGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTAc                               <  1:348483/143‑1 (MQ=255)
                                                                                                  gCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAg                             <  2:342809/143‑1 (MQ=255)
                                                                                                   cttACCCCATCCCCGCCGTGTAGGCCGGATATGACGCATGTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGc                            <  2:73209/141‑1 (MQ=255)
                                                                                                   ccTACCACATCATCGCCGGGTAGGCCGTATAAGCCGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGc                            <  2:73260/143‑1 (MQ=255)
                                                                                                   ccTACCACATCATCACCGTGTCGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGc                            <  2:480456/143‑1 (MQ=255)
                                                                                                       ccACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACGGCTCaa                        >  1:279266/1‑143 (MQ=255)
                                                                                                                  ccGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGa             >  1:724232/1‑143 (MQ=255)
                                                                                                                        aGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCCGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGACGCggg       >  2:798778/1‑142 (MQ=255)
                                                                                                                        aGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGa       <  2:386206/143‑1 (MQ=255)
                                                                                                                         ggCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTGGACGGTGTAGTACAGCTCAAACTCTACCGGATGCGGAg      <  1:553811/143‑1 (MQ=255)
                                                                                                                         cgCCGGATGCGACGCAAATGCGTCTTATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGAACCTCAGTCGGTGGAACACGGTAGTACGGCACAAACTCTACCGGATGCGGAg      >  1:583670/2‑143 (MQ=255)
                                                                                                                          gCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGCCAAGTTTCCACGGCAACTAACACACTTCGCCGCTGTAGTCCCGCCCACGCTCTACCGGACGCggcgc     >  2:236826/1‑140 (MQ=255)
                                                                                                                          gCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCCACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGt     >  1:177151/1‑143 (MQ=255)
                                                                                                                           ccGGATGCGACGCAAATGCGTCTTATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACGCTACTGGATGCGGAgag    >  2:118605/1‑141 (MQ=255)
                                                                                                                             ggATGAGAAGCATCTGCGTCGCATCCGGCGTGATGTGCTGGCCACCGAGCCTAGGCCTGGCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCGCAAACTCTACCGGATGCGGAGTCAt  <  2:483135/143‑1 (MQ=255)
                                                                                                                             ggATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCCACTCCGACACTTAGACCGTGTACAACATCCCAAACGCCAACTGATGCGGAGTAAt  >  1:298905/1‑143 (MQ=255)
                                                                                                                             ggATAAGACGCATTTGCGTCGAATCCGGCGTGATGTGCTTGCCACCGATCCCAGCCCAGTCAGAGACAGGCGAAAAGATACCACAGCAACTAAAACACCTCGACCCCGTAGCAACGCTCAAACTCCACCGGATGCgcacgcaa  >  2:497667/1‑136 (MQ=255)
                                                                                                                                         |                                                                                                                                  
CGCGAAGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCAT  >  NC_000913/4056402‑4056669

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: