Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,053,516:1 (G)6→7 100% coding (2663/2745 nt) torS ← hybrid sensory histidine kinase in two‑component regulatory system with TorR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,053,5101.G91.3% 253.8 / 18.3 80coding (2669/2745 nt)torShybrid sensory histidine kinase in two‑component regulatory system with TorR
Reads supporting (aligned to +/- strand):  ref base . (6/1);  new base G (36/37);  total (42/38)
Fisher's exact test for biased strand distribution p-value = 1.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

GTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGG  >  NC_000913/1053411‑1053606
                                                                                                     |                                                                                                 
gttCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGaa                                                                                                    >  2:840040/1‑100 (MQ=255)
  tCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑gg                                                                                                >  1:865859/1‑100 (MQ=255)
  tCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑gg                                                                                                >  1:1403405/1‑100 (MQ=255)
    ttAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑gggg                                                                                              <  2:134660/100‑1 (MQ=255)
      aaCTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTTTTATGCAACCATGCTTCCAGATCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑gggggg                                                                                            <  2:1513268/100‑1 (MQ=255)
        cTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑gggggg                                                                                            <  1:2345289/98‑1 (MQ=255)
        cTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTGATGCAACCATGCTTCCAGTGCGGCAACACTGCGTGTAATTTCTTCGTGGGGAA‑‑ggggg                                                                                             <  1:2317452/97‑1 (MQ=255)
          ggCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑gggggggcag                                                                                        >  2:1028453/1‑96 (MQ=255)
          ggCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑gggggggcag                                                                                        <  2:433012/100‑5 (MQ=255)
          ggCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑gggggggcag                                                                                        >  2:2477416/1‑96 (MQ=255)
           gCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTCTGCAACCATGCTTCCAGAGCGGCAACCCTGCGTGTAATTTCTTCGTGTGGAA‑‑gggggg                                                                                            >  1:464583/1‑95 (MQ=255)
           gCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGAAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑ggggggggca                                                                                        <  1:140079/99‑5 (MQ=255)
            cTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGAt                                                                                       <  1:593498/100‑1 (MQ=255)
            cTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGAt                                                                                       >  2:2012036/1‑100 (MQ=255)
                tatTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATg                                                                                   <  1:1450215/100‑1 (MQ=255)
                tatTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATg                                                                                   <  1:1699167/100‑1 (MQ=255)
                tatTGATTCAAAGCGCGTTCAGGTCTTTCTGATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAAtn                                                                                   <  1:1411636/100‑2 (MQ=255)
                 atTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATgg                                                                                  <  1:2508158/100‑1 (MQ=255)
                 atTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATgg                                                                                  <  2:384966/100‑1 (MQ=255)
                  ttGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGgc                                                                                 <  1:1206591/100‑1 (MQ=255)
                   tgcttCAAAGCGCGTTTAGGTCTTTCGGATGGAACCATGCTTCCAGAGCGGCAGCACTGCGTGTAAGTTCTTCGTGTGGGAGGGGGGGCAGATAATGgct                                                                                <  1:1176646/97‑1 (MQ=255)
                   tGATTCAAATCGCGTTCAGGTCTTTTTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGgcn                                                                                <  1:1366600/100‑2 (MQ=255)
                   tGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGgct                                                                                >  1:1352030/1‑100 (MQ=255)
                   tGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGgct                                                                                <  1:2316757/100‑1 (MQ=255)
                    gATTCAAGTGGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGGTCGGCAACACTGCGGGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGgttg                                                                               >  1:361523/1‑100 (MQ=255)
                    gATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGgctg                                                                               >  1:1348865/1‑100 (MQ=255)
                     aTTCGAATCGCGGTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGgctgc                                                                              <  1:2289425/100‑1 (MQ=255)
                        cAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGgctgctgc                                                                           <  1:78067/100‑1 (MQ=255)
                         aaaTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGgctgctgct                                                                          >  2:2189283/1‑100 (MQ=255)
                         aaaTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGgctgctgct                                                                          >  1:2096989/1‑100 (MQ=255)
                          aaTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTc                                                                         >  2:2236840/1‑100 (MQ=255)
                           aTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTcc                                                                        <  1:554087/100‑1 (MQ=255)
                            tCGCGTTCAGGTCTTTGTTATGCAACCATGCTGCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCa                                                                       <  1:1516216/100‑1 (MQ=255)
                            tCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCa                                                                       <  2:2260431/100‑1 (MQ=255)
                                gTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTtg                                                                   <  2:797256/100‑1 (MQ=255)
                                gTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTtg                                                                   <  2:796713/100‑1 (MQ=255)
                                 ttCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTtgt                                                                  >  2:1997069/1‑100 (MQ=255)
                                  tCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTtgtg                                                                 <  1:909387/100‑1 (MQ=255)
                                    aGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTgcg                                                               >  2:2128226/1‑100 (MQ=255)
                                     ggTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTgcgc                                                              >  2:2571907/1‑100 (MQ=255)
                                     ggTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTgc                                                                >  1:1559207/1‑98 (MQ=255)
                                      gTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCa                                                             >  2:1335183/1‑100 (MQ=255)
                                       tctttcttATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCAt                                                            <  2:9139/100‑1 (MQ=255)
                                         tttcttATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgc                                                         <  1:1625774/100‑1 (MQ=255)
                                          ttcttATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGCGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAg                                                          >  1:1811721/1‑99 (MQ=255)
                                              tATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTCCGTGTGGCAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggc                                                     >  1:1186315/1‑100 (MQ=255)
                                               aTGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggct                                                    >  1:1316441/1‑100 (MQ=255)
                                                 gCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGctggctgg                                                  >  2:1887376/1‑100 (MQ=255)
                                                 gCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGctggctgg                                                  >  2:1402134/1‑100 (MQ=255)
                                                  cAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGctggctgg                                                  >  1:1846634/1‑99 (MQ=255)
                                                  cAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCtggctggc                                                 >  1:1979522/1‑100 (MQ=255)
                                                  cAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCtggctggc                                                 >  1:1978507/1‑100 (MQ=255)
                                                  cAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCtggctggc                                                 <  2:1176318/100‑1 (MQ=255)
                                                  cAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCtggctggc                                                 <  2:1098476/100‑1 (MQ=255)
                                                    aCCATGCTTCCAGAGCGGCCACACTGCGTGTAATTTCTCCGTGTGGAAGGGGGGGCAGATAATGGCTGCTACTCCAGTTGTGCGCATAGCtagctggc                                                 >  1:544750/1‑98 (MQ=255)
                                                    aCCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGa                                               >  2:289854/1‑100 (MQ=255)
                                                    aCCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGa                                               >  1:2453804/1‑100 (MQ=255)
                                                    aCCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGat                                              >  2:113773/1‑100 (MQ=255)
                                                     ccATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGat                                              <  2:233019/100‑1 (MQ=255)
                                                       aTGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGatat                                            <  2:1842875/100‑1 (MQ=255)
                                                         gCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATg                                           <  1:268271/100‑1 (MQ=255)
                                                          cTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCAt                                        >  1:1725364/1‑100 (MQ=255)
                                                            tCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTc                                      >  2:871461/1‑100 (MQ=255)
                                                             ccAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTc                                      <  2:1885322/100‑1 (MQ=255)
                                                              cAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTcc                                     >  2:1381359/1‑100 (MQ=255)
                                                                gagCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCa                                   <  2:1715226/100‑1 (MQ=255)
                                                                   cGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGAc                                <  2:1316448/100‑1 (MQ=255)
                                                                    ggCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACt                               >  1:119984/1‑100 (MQ=255)
                                                                    ggCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACt                               >  1:1824715/1‑100 (MQ=255)
                                                                     gCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTc                              >  1:2104178/1‑100 (MQ=255)
                                                                     gCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGACTGGCGATATGCATTCCCAGACTc                              >  1:254785/1‑100 (MQ=255)
                                                                       aaCACTGCGTGTAATTTCCTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGa                            <  2:1261430/100‑1 (MQ=255)
                                                                          acTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCa                         >  2:366084/1‑100 (MQ=255)
                                                                          acTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCa                         <  2:1358905/100‑1 (MQ=255)
                                                                             gCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACt                      <  2:202094/100‑1 (MQ=255)
                                                                               gtgtAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGc                    >  2:645492/1‑100 (MQ=255)
                                                                                  tAATTTATTCGTGTGGAAGGGGGGGCAGACAATGGCCGCTGCTCAAGTTGTCCGCATGGCTGTCTGGCGATATGCATTCCCAGACTCGAGCAACTGc                    <  1:747877/97‑1 (MQ=255)
                                                                                   aaTTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCtttt                <  1:81976/100‑1 (MQ=255)
                                                                                    aTTTCTTCGTGTGGAAGGGGGGGCAGATGATGGCTGCTGCTCCAGTTGTGCGCATGGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCttttt               <  1:2100110/100‑1 (MQ=255)
                                                                                    aTTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCttttt               >  2:2360054/1‑100 (MQ=255)
                                                                                    aTTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCttttt               <  2:2158262/100‑1 (MQ=255)
                                                                                    aTTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCttttt               >  2:674487/1‑100 (MQ=255)
                                                                                        cttcGTGTGGAAA‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCt           >  2:1698206/1‑100 (MQ=255)
                                                                                        cttcGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCt           <  2:984474/100‑1 (MQ=255)
                                                                                         ttcGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTg          >  1:134813/1‑100 (MQ=255)
                                                                                          tcGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTACTTTTTAGCTGa         <  1:2110070/100‑1 (MQ=255)
                                                                                            gtgtGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATg       <  1:2045738/100‑1 (MQ=255)
                                                                                            gtgtGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATg       >  2:2091767/1‑100 (MQ=255)
                                                                                              gtgGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATgcg     >  2:446022/1‑100 (MQ=255)
                                                                                              gtgGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATgcg     >  1:212650/1‑100 (MQ=255)
                                                                                               tgGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATgcgc    >  2:700406/1‑100 (MQ=255)
                                                                                                ggAA‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCgg  >  2:1294919/1‑100 (MQ=255)
                                                                                                     |                                                                                                 
GTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGG  >  NC_000913/1053411‑1053606

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