Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,056,539 T→A 6.8% 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,5390TA 6.8% 762.1 / 34.3 220intergenic (‑200/+86)glnL/glnAsensory histidine kinase in two‑component regulatory system with GlnG/glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base T (106/99);  new base A (5/10);  total (111/109)
Fisher's exact test for biased strand distribution p-value = 1.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.02e-01

Only 100 of 228 total aligned reads displayed.

CGGATAAAAAGCGCGAAGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCATACGCACGCGGTCATCTT  >  NC_000913/4056391‑4056686
                                                                                                                                                    |                                                                                                                                                   
 ggATAAAAAGCGCGAAGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCAtt                                                                                                                                                    <  4:724278/149‑1 (MQ=255)
 ggATAAAAAGCGCGAAGCAGCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCAtt                                                                                                                                                    >  2:1090282/1‑149 (MQ=255)
    tAAAAAGCGCGAAGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGc                                                                                                                                                 >  2:1540071/1‑149 (MQ=255)
    tAAAAAGCGCGAAGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATAACCGTGTAGGCCGGATAAGACGCATTTGc                                                                                                                                                 >  4:1462288/1‑149 (MQ=255)
       aaaGCGCGAAGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTc                                                                                                                                              <  3:457715/149‑1 (MQ=255)
       aaaGCGCGAAGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTc                                                                                                                                              <  1:844050/149‑1 (MQ=255)
              gAAGCATCAGAGAATTGCCGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCGACCACATCATCACCGTGTAGGCCGGATAAGACGCAAATGCGTctta                                                                                                                                           <  3:462297/145‑4 (MQ=255)
              gAAGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCg                                                                                                                                       <  4:626130/149‑1 (MQ=255)
              gAAGCATCAGAGAATTGACGGAGAAAAAAGCCCATGAAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCg                                                                                                                                       >  3:1144403/1‑149 (MQ=255)
                    tCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGa                                                                                                                                 >  1:1250584/1‑149 (MQ=255)
                     cAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGAt                                                                                                                                <  2:227093/149‑1 (MQ=255)
                       gagaATTGACGGAGAAAAAAGCCTATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCAttt                                                                                                                                                   >  2:838967/1‑128 (MQ=255)
                          aaTTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCt                                                                                                                           <  4:743421/149‑1 (MQ=255)
                           aTTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCtt                                                                                                                          >  3:1240585/1‑149 (MQ=255)
                             tGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGc                                                                                                                        >  1:1028837/1‑149 (MQ=255)
                               aCGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCa                                                                                                                      <  1:473874/149‑1 (MQ=255)
                               aCGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCa                                                                                                                      <  2:183026/149‑1 (MQ=255)
                                 ggAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGAt                                                                                                                                >  4:448382/1‑137 (MQ=255)
                                 ggAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCAcc                                                                                                                    >  2:512859/1‑149 (MQ=255)
                                  gagaAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCg                                                                                                                   <  4:1074668/149‑1 (MQ=255)
                                    gaAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGAt                                                                                                                 >  4:1504211/1‑149 (MQ=255)
                                      aaaaaGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATcc                                                                                                               <  3:128951/149‑1 (MQ=255)
                                       aaaaGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGAt                                                                                                                                <  1:99769/131‑1 (MQ=255)
                                       aaaaGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATccc                                                                                                              >  3:1303306/1‑149 (MQ=255)
                                        aaaGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCa                                                                                                             <  4:1167465/149‑1 (MQ=255)
                                        aaaGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCa                                                                                                             <  3:888498/149‑1 (MQ=255)
                                         aaGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAg                                                                                                            <  4:720072/149‑1 (MQ=255)
                                         aaGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAg                                                                                                            <  3:1462288/149‑1 (MQ=255)
                                              cATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTg                                                                                                       >  3:57989/1‑149 (MQ=255)
                                               aTGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGc                                                                                                      <  1:649394/149‑1 (MQ=255)
                                               aTGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGc                                                                                                      >  4:511028/1‑149 (MQ=255)
                                               aTGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGc                                                                                                      >  3:421545/1‑149 (MQ=255)
                                               aTGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCAGCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGc                                                                                                      <  4:1595324/149‑1 (MQ=255)
                                                 gCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATGCGACGCAAATGCGTcttat                                                                                                                                          <  4:1294574/111‑5 (MQ=255)
                                                 gCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCAAATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCa                                                                                                    >  3:1173987/1‑149 (MQ=255)
                                                     agaTGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAgaga                                                                                                <  2:1266714/149‑1 (MQ=255)
                                                     agaTGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAgaga                                                                                                <  1:497292/149‑1 (MQ=255)
                                                      gaTGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGAc                                                                                               <  4:538833/149‑1 (MQ=255)
                                                         gggCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAgg                                                                                            >  3:573738/1‑149 (MQ=255)
                                                            cTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGa                                                                                         >  1:1081344/1‑149 (MQ=255)
                                                            cTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGa                                                                                         <  2:605111/149‑1 (MQ=255)
                                                              aCAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGaaa                                                                                       <  1:1218880/149‑1 (MQ=255)
                                                               cAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGaaaa                                                                                      <  4:681070/149‑1 (MQ=255)
                                                                    aGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTc                                                                                 <  4:1139098/149‑1 (MQ=255)
                                                                       tGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCAc                                                                              <  1:310010/149‑1 (MQ=255)
                                                                           aaaCTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCa                                                                          >  1:567800/1‑149 (MQ=255)
                                                                             aCTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAAc                                                                        <  4:859263/149‑1 (MQ=255)
                                                                               ttCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTa                                                                      >  2:719483/1‑149 (MQ=255)
                                                                               ttCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTa                                                                      >  3:287832/1‑149 (MQ=255)
                                                                                  aCGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTaaaa                                                                   <  3:621979/149‑1 (MQ=255)
                                                                                  aCGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTaaaa                                                                   <  3:1667874/149‑1 (MQ=255)
                                                                                          gagCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTc                                                                                 >  3:1537326/1‑127 (MQ=255)
                                                                                              gACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCt                                                       <  4:363177/149‑1 (MQ=255)
                                                                                               aCTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTg                                                      >  4:249205/1‑149 (MQ=255)
                                                                                                 tCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAacac                                                                >  1:740479/1‑137 (MQ=255)
                                                                                                 tCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTgta                                                    >  4:1541004/1‑149 (MQ=255)
                                                                                                 tCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTgta                                                    >  1:1004998/1‑149 (MQ=255)
                                                                                                   aCCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAg                                                                                             >  2:206789/1‑106 (MQ=255)
                                                                                                   aCCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTgtagt                                                  >  1:1402362/1‑149 (MQ=255)
                                                                                                    ccTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAg                                                                                                            <  4:140420/90‑1 (MQ=255)
                                                                                                    ccTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTgtagta                                                 >  4:786610/1‑149 (MQ=255)
                                                                                                     cTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTAc                                                >  3:1389754/1‑149 (MQ=255)
                                                                                                      tGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCa                                                                                                             >  4:1143215/1‑87 (MQ=255)
                                                                                                         tCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGGAGTACAGCt                                            >  3:364406/1‑149 (MQ=255)
                                                                                                                tACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAActc                                     >  2:1047169/1‑149 (MQ=255)
                                                                                                                tACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAActc                                     >  2:568077/1‑149 (MQ=255)
                                                                                                                 aCCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAActct                                    <  3:786610/149‑1 (MQ=255)
                                                                                                                  ccACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTa                                   <  3:11410/149‑1 (MQ=255)
                                                                                                                     catcatCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAg                                                                                             <  2:1693504/88‑1 (MQ=255)
                                                                                                                     catcatCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAg                                                                                             >  1:1693504/1‑88 (MQ=255)
                                                                                                                      atcatcACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCgg                               >  1:1614521/1‑149 (MQ=255)
                                                                                                                        catcaCCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGAt                             <  4:1240585/149‑1 (MQ=255)
                                                                                                                         atcaCCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATg                            >  2:1466022/1‑149 (MQ=255)
                                                                                                                          tcaCCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGc                           <  3:307191/149‑1 (MQ=255)
                                                                                                                          tcaCCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGc                           >  1:512140/1‑149 (MQ=255)
                                                                                                                          tcaCCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGc                           >  4:1385976/1‑149 (MQ=255)
                                                                                                                           caCCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCg                          <  1:767618/149‑1 (MQ=255)
                                                                                                                           caCCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCg                          >  2:1171507/1‑149 (MQ=255)
                                                                                                                           caCCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGGACAGCTCAAACTCTACCGGATGCg                          >  2:1119896/1‑149 (MQ=255)
                                                                                                                            aCCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCgg                         >  2:692381/1‑149 (MQ=255)
                                                                                                                             ccGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGa                        <  1:684892/149‑1 (MQ=255)
                                                                                                                                tgtAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTc                     <  4:1181876/149‑1 (MQ=255)
                                                                                                                                    ggCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCATac                 <  3:1279589/149‑1 (MQ=255)
                                                                                                                                    ggCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCATac                 <  1:1058184/149‑1 (MQ=255)
                                                                                                                                     gCCGGATGCGACGCAAATGCGTCTTATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCATacg                >  2:441419/1‑149 (MQ=255)
                                                                                                                                         gATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCATacgcacg            >  2:988394/1‑149 (MQ=255)
                                                                                                                                         gATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCATacgcacg            >  2:324420/1‑149 (MQ=255)
                                                                                                                                            gcgACGCAAATGCGTCTTATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGc                                                                                           <  1:800585/65‑1 (MQ=25)
                                                                                                                                            gcgACGCAAATGCGTCTTATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGc                                                                                           >  2:800585/3‑67 (MQ=25)
                                                                                                                                            aaGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCATACGCACGCgg         <  4:1014679/149‑1 (MQ=255)
                                                                                                                                              gACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGaaaa                                                                                      >  3:582039/1‑70 (MQ=255)
                                                                                                                                              gACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGaaaa                                                                                      <  4:582039/70‑1 (MQ=255)
                                                                                                                                              gACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCATACGCACGCGGTc       >  3:636579/1‑149 (MQ=255)
                                                                                                                                              gACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCATACGCACGCGGTc       >  4:542713/1‑149 (MQ=255)
                                                                                                                                              gACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCATACGCACGCGGTc       >  3:873369/1‑149 (MQ=255)
                                                                                                                                               aCGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCATACGCACGCGGTCa      >  4:373125/1‑149 (MQ=255)
                                                                                                                                               aCGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCATACGCACGCGGTCa      >  1:1346919/1‑149 (MQ=255)
                                                                                                                                               aCGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCATACGCAAGCGGTCa      <  3:697047/149‑1 (MQ=255)
                                                                                                                                                 gCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTgta                                                    >  3:1082441/1‑101 (MQ=255)
                                                                                                                                                  cATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCATACGCACGCGGTCAtct   <  3:720679/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGGATAAAAAGCGCGAAGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCATACGCACGCGGTCATCTT  >  NC_000913/4056391‑4056686

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