Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,486,083 T→C 5.5% D570D (GAT→GAC evgS → hybrid sensory histidine kinase in two‑component regulatory system with EvgA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,486,0830TC 5.5% 502.7 / 12.2 148D570D (GAT→GACevgShybrid sensory histidine kinase in two‑component regulatory system with EvgA
Reads supporting (aligned to +/- strand):  ref base T (66/73);  new base C (4/4);  total (71/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 154 total aligned reads displayed.

TGCCCAATGTGACCATTGACACATGGGACCTATATAGCGAGCAATTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACATTATTTCACTGCGGATTACTACAAAAATGCAATGTTACC  >  NC_000913/2485935‑2486232
                                                                                                                                                    |                                                                                                                                                     
tGCCCAATGTGACCATTGACACATTGGACCTATATAGCGAGCAATTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGAtt                                                                                                                                                      <  2:1380867/150‑1 (MQ=255)
tGCCCAATGTGACCATTGACACATGGGACCTATATAGCGAGCAATTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGAtt                                                                                                                                                      <  2:638746/150‑1 (MQ=255)
tGCCCAATGTGACCATTGACACATGGGACCTATATAGCGAGCAATTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGAtt                                                                                                                                                      <  2:1221944/150‑1 (MQ=255)
   ccAATGTGACCATTGACACATGGGACCTATATAGCGAGCAATTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAg                                                                                                                                                   <  2:2874961/150‑1 (MQ=255)
   ccAATGTGACCATTGACACATGGGACCTATATAGCGAGCAATTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAg                                                                                                                                                   <  2:2332339/150‑1 (MQ=255)
   ccAATGTGACCATTGACACATGGGACCTATATAGCGAGCAATTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAg                                                                                                                                                   <  2:1534387/150‑1 (MQ=255)
    cAATGTGACCATTGACACATGGGACCTATATAGCGAGCAATTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGa                                                                                                                                                  >  1:2088972/1‑150 (MQ=255)
    cAATGTGACCATTGACACATGGGACCTATATAGCGAGCAATTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGa                                                                                                                                                  <  1:845469/150‑1 (MQ=255)
       tgtgACCATTGACACATGGGACCTATATAGCGAGCAATTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGaaaa                                                                                                                                               <  2:2516673/150‑1 (MQ=255)
            ccATTGACACATGGGACCTATATAGCGAGCAATTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCaaa                                                                                                                                          <  2:914787/150‑1 (MQ=255)
            ccATTGACACATGGGACCTATATAGCGAGCAATTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCaaa                                                                                                                                          <  1:818496/150‑1 (MQ=255)
             cATTGACACATGGGACCTATATAGCGAGCAATTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAat                                                                                                                                         >  1:211662/1‑150 (MQ=255)
              aTTGACACATGGGACCTATATAGCGAGCAATTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAata                                                                                                                                        >  2:1561517/1‑150 (MQ=255)
              aTTGACACATGGGACCTATATAGCGAGCAATTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAata                                                                                                                                        >  2:360376/1‑150 (MQ=255)
                tGACACATGGGACCTATATAGCGAGCAATTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATc                                                                                                                                      >  1:77799/1‑150 (MQ=255)
                   cacaTGGGACCTATATAGCGAGCAATTTTTTATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAATACCAAATATCAtt                                                                                                                                   <  1:1439104/150‑1 (MQ=255)
                    acaTGGGACCTATATAGCGAGCAATTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTc                                                                                                                                  >  1:2137118/1‑150 (MQ=255)
                       tGGGACCTATATAGCGAGCAATTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGa                                                                                                                               <  1:785903/150‑1 (MQ=255)
                          gACCTATATAGCGAGCAATTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGaaaa                                                                                                                            <  2:269024/150‑1 (MQ=255)
                                 ataGCGAGCAATTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCActct                                                                                                                     <  2:477528/150‑1 (MQ=255)
                                 ataGCGAGCAATTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCActct                                                                                                                     <  2:477883/150‑1 (MQ=255)
                                    gCGAGCAATTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCgg                                                                                                                  >  1:1911971/1‑150 (MQ=255)
                                    gCGAGCAATTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCgg                                                                                                                  >  1:1551665/1‑150 (MQ=255)
                                      gAGCATTTGTATCTTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCCAAGATCATTCCGAAAAGCACTCTCGGAt                                                                                                                <  2:2270020/150‑1 (MQ=255)
                                      gAGCAATTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGAt                                                                                                                >  1:2774549/1‑150 (MQ=255)
                                        gCAATTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTc                                                                                                              <  1:92267/150‑1 (MQ=255)
                                        gCAATTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTc                                                                                                              <  1:2412386/150‑1 (MQ=255)
                                        gCAATTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTc                                                                                                              <  1:1083746/150‑1 (MQ=255)
                                         cAATTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTcc                                                                                                             >  2:2576594/1‑150 (MQ=255)
                                          aaTTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGACTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCt                                                                                                            >  2:2623055/1‑150 (MQ=255)
                                          aaTTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGACTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCt                                                                                                            >  2:1412947/1‑150 (MQ=255)
                                            ttttATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTa                                                                                                          >  2:684200/1‑150 (MQ=255)
                                              ttATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTATATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTAcc                                                                                                        >  2:1470967/1‑150 (MQ=255)
                                              ttATATTGTTACGACATTATCCGGTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGATCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTAcc                                                                                                        >  2:2459792/1‑150 (MQ=255)
                                               tataTTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCg                                                                                                       >  2:2302104/1‑150 (MQ=255)
                                                  atTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAAt                                                                                                    <  2:1334961/150‑1 (MQ=255)
                                                   ttGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATc                                                                                                   <  2:326750/150‑1 (MQ=255)
                                                    tGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATcc                                                                                                  >  1:73740/1‑150 (MQ=255)
                                                      ttACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCaa                                                                                                <  1:2791994/150‑1 (MQ=255)
                                                      ttACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCaa                                                                                                <  1:338973/150‑1 (MQ=255)
                                                      ttACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGCTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCaa                                                                                                <  1:964598/150‑1 (MQ=255)
                                                       tACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCGGTAAGGAATccgag                                                                                               <  2:2929633/150‑4 (MQ=255)
                                                          gACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTa                                                                                            >  1:1356555/1‑150 (MQ=255)
                                                          gACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTa                                                                                            >  1:2765818/1‑150 (MQ=255)
                                                                 tCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTa                                                                                     >  1:1005664/1‑150 (MQ=255)
                                                                  ccGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTaa                                                                                    <  2:1599624/150‑1 (MQ=255)
                                                                   cGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTaaa                                                                                   >  2:2839933/1‑150 (MQ=255)
                                                                    gTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAAc                                                                                  >  1:1556058/1‑150 (MQ=255)
                                                                     ttttATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACt                                                                                 >  2:1718988/1‑150 (MQ=255)
                                                                      tttATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTg                                                                                <  2:2088972/150‑1 (MQ=255)
                                                                      tttATTAGTTGGCAGTAGCCATTTATGGGGATTCTACCTCTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTg                                                                                <  1:2333113/150‑1 (MQ=255)
                                                                       ttattaGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTgg                                                                               <  1:1099312/150‑1 (MQ=255)
                                                                          ttaGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCaa                                                                            >  1:1713588/1‑150 (MQ=255)
                                                                          ttaGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCaa                                                                            <  1:1827921/150‑1 (MQ=255)
                                                                           taGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAg                                                                           <  1:977399/150‑1 (MQ=255)
                                                                           taGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAg                                                                           <  1:1932392/150‑1 (MQ=255)
                                                                              ttGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTa                                                                        >  1:806925/1‑150 (MQ=255)
                                                                              ttGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTa                                                                        >  1:2702979/1‑150 (MQ=255)
                                                                                ggCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAAt                                                                      <  1:540223/150‑1 (MQ=255)
                                                                                ggCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAAt                                                                      <  1:1512821/150‑1 (MQ=255)
                                                                                    gTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCa                                                                  <  2:2659001/150‑1 (MQ=255)
                                                                                        ccTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAg                                                              >  1:691271/1‑150 (MQ=255)
                                                                                           tttATGGGGATTCTACCTTTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAAATTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCa                                                           <  2:1438285/150‑1 (MQ=255)
                                                                                           tttATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCa                                                           <  2:1387962/150‑1 (MQ=255)
                                                                                               tGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCataat                                                       >  1:2692063/1‑150 (MQ=255)
                                                                                                 gggATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAg                                                     >  2:1240175/1‑150 (MQ=255)
                                                                                                 gggATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAg                                                     >  2:707311/1‑150 (MQ=255)
                                                                                                     ttCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCt                                                 <  1:1044741/150‑1 (MQ=255)
                                                                                                       cTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCttt                                               >  2:1991695/1‑150 (MQ=255)
                                                                                                        tACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCtttt                                              >  1:2306030/1‑150 (MQ=255)
                                                                                                             gTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGAATTATAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACATATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACa                                         >  1:2464040/1‑150 (MQ=255)
                                                                                                                aCGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACatta                                      <  1:1322901/150‑1 (MQ=255)
                                                                                                                aCGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACatta                                      <  1:138891/150‑1 (MQ=255)
                                                                                                                aCGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACatta                                      <  1:1998862/150‑1 (MQ=255)
                                                                                                                aCGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACatta                                      <  1:1998833/150‑1 (MQ=255)
                                                                                                                 cGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACattat                                     <  1:943243/150‑1 (MQ=255)
                                                                                                                  gCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACattatt                                    >  1:2464876/1‑150 (MQ=255)
                                                                                                                    tCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACATTATTTc                                  >  2:868640/1‑150 (MQ=255)
                                                                                                                    tCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACATTATTTc                                  >  1:2242441/1‑150 (MQ=255)
                                                                                                                    tCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACATTATTTc                                  >  1:1788149/1‑150 (MQ=255)
                                                                                                                     cAGTTCGTCGTCGTAAAGTCATTCAGGTTGATTTAGAAAACCAAATATCATTCCTAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACATTATTTCa                                 <  2:190944/150‑1 (MQ=255)
                                                                                                                     cAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACATTATTTCa                                 <  2:1682091/150‑1 (MQ=255)
                                                                                                                     cAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACATTATTTCa                                 <  2:825060/150‑1 (MQ=255)
                                                                                                                     cAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACATTATTTCa                                 <  2:1911971/150‑1 (MQ=255)
                                                                                                                       gTTCGTCGTCGTAAAGTCATTCAGGGTGACTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACATTATTTCACt                               <  1:1412947/150‑1 (MQ=255)
                                                                                                                       gTTCGTCGTCGTAAAGTCATTCAGGGTGACTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACATTATTTCACt                               <  1:2623055/150‑1 (MQ=255)
                                                                                                                       gTTCGTCGTCGTAAAGTCATTCAGGGTGACTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTGACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACATTATTTCACt                               <  1:2623044/150‑1 (MQ=255)
                                                                                                                           gtcgtcgtAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACATTATTTCACTGCgg                           <  2:2583934/150‑1 (MQ=255)
                                                                                                                           gtcgtcgtAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACATTATTTCACTGCgg                           >  2:629914/1‑150 (MQ=255)
                                                                                                                                 gtAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACATTATTTCACTGCGGATtact                     <  1:2399536/150‑1 (MQ=255)
                                                                                                                                  tAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCCTTTGAACATTATTTCACTGCGGATtacta                    >  2:2285795/1‑150 (MQ=255)
                                                                                                                                          ttCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACATTATTTCACTGCGGATTACTACAAAAATg            <  1:2302104/150‑1 (MQ=255)
                                                                                                                                          ttCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACATTATTTCACTGCGGATTACTACAAAAATg            <  2:1556058/150‑1 (MQ=255)
                                                                                                                                          ttCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACATTATTTCACTGCGGATTACTACAAAAATg            >  2:2194940/1‑150 (MQ=255)
                                                                                                                                           tCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACATTATTTCACTGCGGATTACTACAAAAATGc           >  1:2577490/1‑150 (MQ=255)
                                                                                                                                           tCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACATTATTTCACTGCGGATTACTACAAAAATGc           >  1:2577484/1‑150 (MQ=255)
                                                                                                                                            cAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACATTATTTCACTGCGGATTACTACAAAAATGCa          <  1:459504/150‑1 (MQ=255)
                                                                                                                                                   aTTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACATTATTTCACTGCGGATTACTACAAAAATGCAATGTTAc   <  2:1723082/150‑1 (MQ=255)
                                                                                                                                                    tttAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACATTATTTCACTGCGGATTACTACAAAAATGCAATGTTAcc  <  1:2739316/150‑1 (MQ=255)
                                                                                                                                                    tttAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACATTATTTCACTGCGGATTACTACAAAAATGCAATGTTAcc  >  2:2739316/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGCCCAATGTGACCATTGACACATGGGACCTATATAGCGAGCAATTTTATATTGTTACGACATTATCCGTTTTATTAGTTGGCAGTAGCCTTTTATGGGGATTCTACCTGTTACGCTCAGTTCGTCGTCGTAAAGTCATTCAGGGTGATTTAGAAAACCAAATATCATTCCGAAAAGCACTCTCGGATTCCTTACCGAATCCAACTTATGTTGTAAACTGGCAAGGTAATGTCATTAGTCATAATAGTGCTTTTGAACATTATTTCACTGCGGATTACTACAAAAATGCAATGTTACC  >  NC_000913/2485935‑2486232

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