Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,056,531 A→G 7.1% intergenic (‑192/+94) 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,5310AG 7.1% 188.0 / 3.2 66intergenic (‑192/+94)glnL/glnAsensory histidine kinase NtrB/glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base A (33/28);  new base G (2/3);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.41e-01

AAAAAGCGCGAAGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCAT  >  NC_000913/4056396‑4056669
                                                                                                                                       |                                                                                                                                          
aaaaaGCGCGAAGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGGAGGCCGGATAAGACGCa                                                                                                                                     >  1:408519/1‑143 (MQ=255)
      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                                   <  1:115055/143‑1 (MQ=255)
                                                                                                  gCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTgta                                   <  2:227974/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)
                                                                                                                                       |                                                                                                                                          
AAAAAGCGCGAAGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCAT  >  NC_000913/4056396‑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: