Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,056,539 T→A 13.2% intergenic (‑200/+86) glnL ← / ← glnA sensory histidine kinase in two‑component regulatory system with GlnG/glutamine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,056,5390TA13.2% 199.2 / 21.8 68intergenic (‑200/+86)glnL/glnAsensory histidine kinase in two‑component regulatory system with GlnG/glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base T (31/28);  new base A (3/6);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.98e-01

GGATAAAAAGCGCGAAGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCATACGCACGCGGTCATCTT  >  NC_000913/4056392‑4056686
                                                                                                                                                   |                                                                                                                                                   
ggATAAAAAGCGCGAAGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCAtt                                                                                                                                                    >  2:61845/1‑149 (MQ=255)
       aaGCGCGAAGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCg                                                                                                                                             >  3:201750/1‑149 (MQ=255)
           gcgAAGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATGCGACGCAAATGCGTcttatc                                                                                                                                         <  4:92666/149‑6 (MQ=255)
            cgAAGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATcc                                                                                                                                        >  2:219443/1‑149 (MQ=255)
             gAAGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATGCGACGCAAATGCGTcttatccg                                                                                                                                       <  2:25563/149‑8 (MQ=255)
              aaGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATGCGACGCAAATGCGTCTTATCCgg                                                                                                                                      <  3:374164/149‑1 (MQ=255)
                   tCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATGCGACGCAAATGCGTcttatccg                                                                                                                                       >  2:115742/1‑136 (MQ=255)
                   tCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATGCGACGCAAATGCGTCTTATCCg                                                                                                                                       <  1:115742/143‑1 (MQ=255)
                   tCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGa                                                                                                                                 >  3:17480/1‑149 (MQ=255)
                     agagAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGAtg                                                                                                                               >  1:260444/1‑149 (MQ=255)
                      gagaATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGAtgt                                                                                                                              >  1:327941/1‑149 (MQ=255)
                        gaATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGc                                                                                                                            <  2:380133/149‑1 (MQ=255)
                         aaTTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCt                                                                                                                           <  1:133797/149‑1 (MQ=255)
                              aCGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCa                                                                                                                      >  3:288461/1‑149 (MQ=255)
                                ggAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTc                                                                                                                                              <  1:292323/123‑1 (MQ=255)
                                ggAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTc                                                                                                                                              >  2:292323/1‑123 (MQ=255)
                                   gaAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGAt                                                                                                                 >  3:417391/1‑149 (MQ=255)
                                   gaAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGAt                                                                                                                 <  2:331217/149‑1 (MQ=255)
                                      aaaaGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATccc                                                                                                              >  3:222758/1‑149 (MQ=255)
                                        aaGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAg                                                                                                            <  4:288461/149‑1 (MQ=255)
                                         aGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACACCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAgg                                                                                                           <  1:75669/149‑1 (MQ=255)
                                           cccATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTTTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATCTGCtt                                                                                                                          >  3:210393/1‑132 (MQ=255)
                                           cccATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCtt                                                                                                                          <  4:210393/132‑1 (MQ=255)
                                                      aTGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACa                                                                                              >  1:469159/1‑149 (MQ=255)
                                                         ggCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGc                                                                                           <  1:61845/149‑1 (MQ=255)
                                                         ggCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGc                                                                                           >  3:451375/1‑149 (MQ=255)
                                                           cTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGa                                                                                         <  1:325556/149‑1 (MQ=255)
                                                            tACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTACGCCGGATGCGACGCAAATGCGTcttat                                                                                                                                          <  1:196667/99‑5 (MQ=37)
                                                            tACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTACGCCGGATGCGACGCAAATGCGTcttat                                                                                                                                          >  2:196667/1‑95 (MQ=37)
                                                               aGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAg                                                                                     >  3:362859/1‑149 (MQ=255)
                                                                gATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGt                                                                                    <  4:17480/149‑1 (MQ=255)
                                                                  tAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGttt                                                                                  <  3:319680/149‑1 (MQ=255)
                                                                   aGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTc                                                                                 <  2:260444/149‑1 (MQ=255)
                                                                     cTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCa                                                                               <  1:427461/149‑1 (MQ=255)
                                                                      tGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCAc                                                                              >  3:238467/1‑149 (MQ=255)
                                                                         cAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGc                                                                           <  4:451375/149‑1 (MQ=255)
                                                                              ttCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAgg                                                                                                           >  2:378881/1‑112 (MQ=255)
                                                                              ttCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAgg                                                                                                           <  1:378881/112‑1 (MQ=255)
                                                                                 aCGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTaaaa                                                                   <  1:243984/149‑1 (MQ=255)
                                                                                      ggAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACtt                                                              <  4:238467/149‑1 (MQ=255)
                                                                                         gagCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGa                                                           >  2:48003/1‑149 (MQ=255)
                                                                                             gACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCt                                                       >  4:210351/1‑149 (MQ=255)
                                                                                             gACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCt                                                       <  4:417391/149‑1 (MQ=255)
                                                                                              aCTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTg                                                      <  2:327941/149‑1 (MQ=255)
                                                                                                  aCCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTgtagt                                                  >  2:411156/1‑149 (MQ=255)
                                                                                                   ccTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTgtagta                                                 >  3:367736/1‑149 (MQ=255)
                                                                                                    cTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAg                                                                                                            >  3:167264/1‑89 (MQ=255)
                                                                                                    cTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAg                                                                                                            <  4:167264/89‑1 (MQ=255)
                                                                                                          cGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCa                                          <  2:469159/149‑1 (MQ=255)
                                                                                                             ccTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAAc                                       <  1:411156/149‑1 (MQ=255)
                                                                                                             ccTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAAc                                       >  4:329289/1‑149 (MQ=255)
                                                                                                             ccTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAAc                                       <  4:367736/149‑1 (MQ=255)
                                                                                                              cTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAAct                                      >  1:251672/1‑149 (MQ=255)
                                                                                                              cTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAAct                                      >  4:462987/1‑149 (MQ=255)
                                                                                                                aCCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGc                                             <  1:181092/140‑1 (MQ=255)
                                                                                                                aCCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGc                                             >  2:181092/1‑140 (MQ=255)
                                                                                                                 ccACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTa                                   >  3:326654/1‑149 (MQ=255)
                                                                                                                       catcaCCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAAc                                       >  1:418944/1‑139 (MQ=255)
                                                                                                                       catcaCCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAAc                                       <  2:418944/139‑1 (MQ=255)
                                                                                                                          caCCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCg                          >  4:155677/1‑149 (MQ=255)
                                                                                                                            ccGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGa                        <  1:219443/149‑1 (MQ=255)
                                                                                                                                  aGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCATa                  <  4:222758/149‑1 (MQ=255)
                                                                                                                                   ggCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCAtgc                 <  4:326654/149‑3 (MQ=255)
                                                                                                                                    gCCGGATGCGACGCAAATGCGTCTTATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTaaa                                                                    >  2:140573/1‑97 (MQ=25)
                                                                                                                                    gCCGGATGCGACGCAAATGCGTCTTATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTaaa                                                                    <  1:140573/97‑1 (MQ=37)
                                                                                                                                             gACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCATACGCACGCGGTc       >  4:116977/1‑149 (MQ=255)
                                                                                                                                             gACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCATACGCACGCGGTc       >  4:90139/1‑149 (MQ=255)
                                                                                                                                             gACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCCACCGGATGCGGAGTCATACGCACGCGGTc       >  3:122974/1‑149 (MQ=255)
                                                                                                                                                  aTTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCATACGCACGCGGTCAtctt  >  2:318137/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGATAAAAAGCGCGAAGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCATACGCACGCGGTCATCTT  >  NC_000913/4056392‑4056686

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