Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,053,516 (G)6→7 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.G98.9% 331.1 / ‑2.6 90coding (2669/2745 nt)torShybrid sensory histidine kinase in two‑component regulatory system with TorR
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (47/42);  total (48/42)
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.30e-01

Only 100 of 108 total aligned reads displayed.

GTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGC  >  NC_000913/1053411‑1053604
                                                                                                    |                                                                                               
gttCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGaa                                                                                                  >  1:2559458/1‑100 (MQ=255)
gttCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGaa                                                                                                  >  1:2384019/1‑100 (MQ=255)
 ttCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑g                                                                                               <  2:1946389/100‑1 (MQ=255)
  tCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑gg                                                                                              <  2:2070758/100‑1 (MQ=255)
   cTTAACTGGCTGATATTGACTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑ggg                                                                                             >  1:2146468/1‑100 (MQ=255)
     tAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑ggggg                                                                                           <  2:2595199/100‑1 (MQ=255)
     tAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑ggggg                                                                                           >  2:755488/1‑100 (MQ=255)
     tAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑ggggg                                                                                           >  1:1539667/1‑100 (MQ=255)
     tAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑ggggg                                                                                           >  2:687487/1‑100 (MQ=255)
     tAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑ggggg                                                                                           <  1:2062603/100‑1 (MQ=255)
       aCTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑ggggggg                                                                                         >  1:2185354/1‑99 (MQ=255)
       aCTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑ggggggg                                                                                         <  2:1922226/100‑2 (MQ=255)
        cTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑gggggggc                                                                                        <  2:1540386/100‑3 (MQ=255)
        cTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑gggggggc                                                                                        <  1:202987/100‑3 (MQ=255)
        cTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGCGGAA‑‑gggggggc                                                                                        <  2:1855706/100‑3 (MQ=255)
         tGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑gggggggca                                                                                       <  1:1635940/100‑4 (MQ=255)
            cTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCCACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGAt                                                                                     >  1:647385/1‑100 (MQ=255)
             tGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATa                                                                                    <  2:2571798/100‑1 (MQ=255)
             tGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATa                                                                                    >  2:1127186/1‑100 (MQ=255)
             tGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATa                                                                                    >  2:708374/1‑100 (MQ=255)
              gATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATaa                                                                                   >  2:2545018/1‑100 (MQ=255)
              gATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATaa                                                                                   >  1:1967219/1‑100 (MQ=255)
              gATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATaa                                                                                   <  1:1353744/100‑1 (MQ=255)
                  ttGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGgc                                                                               >  1:800158/1‑100 (MQ=255)
                  ttGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGgc                                                                               >  2:2117060/1‑100 (MQ=255)
                   tGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGgct                                                                              <  2:339362/100‑1 (MQ=255)
                   tGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGgct                                                                              <  2:2292186/100‑1 (MQ=255)
                   tGATTCAAATCGCGTTCAGGTCTTTCGTATGCAACCATGCTTCCAGGGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGgct                                                                              <  1:2583581/100‑1 (MQ=255)
                    gATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGgctg                                                                             >  1:2404842/1‑100 (MQ=255)
                    gATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑GGGGGGCAGATAATGgctgc                                                                            >  2:2091006/1‑100 (MQ=255)
                     aTTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGgctgc                                                                            >  2:568968/1‑100 (MQ=255)
                       tCAAATCTCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGgctgctg                                                                          >  2:1490262/1‑100 (MQ=255)
                         aaaTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGgctgctgct                                                                        <  1:2445009/100‑1 (MQ=255)
                           aTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTcc                                                                      >  2:2697181/1‑100 (MQ=255)
                           aTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTcc                                                                      <  2:1889005/100‑1 (MQ=255)
                            tCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCa                                                                     <  2:2265928/100‑1 (MQ=255)
                              gcgTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGt                                                                   <  1:520616/100‑1 (MQ=255)
                              gcgTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGt                                                                   >  2:267803/1‑100 (MQ=255)
                               cgTTGAGGTCTTTCTTATGCAACCATGCATCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAAATAATGGCTGCTGCTCCAGtt                                                                  <  1:195205/100‑1 (MQ=255)
                                gTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTtg                                                                 >  2:1891932/1‑100 (MQ=255)
                                gTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTtg                                                                 <  2:71942/100‑1 (MQ=255)
                                gTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTtg                                                                 >  2:898534/1‑100 (MQ=255)
                                gTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTtg                                                                 <  1:1106112/100‑1 (MQ=255)
                                  tCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTtgtg                                                               <  2:271427/100‑1 (MQ=255)
                                  tCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTtgtg                                                               >  1:513978/1‑100 (MQ=255)
                                   cAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTgc                                                              <  1:2330402/100‑1 (MQ=255)
                                   cAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTgc                                                              >  1:2346153/1‑100 (MQ=255)
                                     ggTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTgcgc                                                            <  2:2071366/100‑1 (MQ=255)
                                     ggTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTgcgc                                                            >  2:520369/1‑100 (MQ=255)
                                     ggTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTgcgc                                                            >  2:632816/1‑100 (MQ=255)
                                      gTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCa                                                           >  2:1411280/1‑100 (MQ=255)
                                        ctttcttATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATa                                                         >  2:2324730/1‑100 (MQ=255)
                                         tttcttATGCAACCATGCTTCCAGAGCGGGAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATa                                                         <  1:526029/99‑1 (MQ=255)
                                           tcttATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgct                                                      >  2:1407363/1‑100 (MQ=255)
                                            cttATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctg                                                     >  2:927146/1‑100 (MQ=255)
                                            cttATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctg                                                     >  1:1791863/1‑100 (MQ=255)
                                             ttATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctgg                                                    <  1:755482/100‑1 (MQ=255)
                                              tATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCCGTTGTGCGCATAgctggc                                                   >  2:1863945/1‑100 (MQ=255)
                                               aTGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggct                                                  >  1:799383/1‑100 (MQ=255)
                                                tGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggctg                                                 >  2:188076/1‑100 (MQ=255)
                                                tGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggctg                                                 >  2:1593631/1‑100 (MQ=255)
                                                  cAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCtggctggc                                               >  2:1234551/1‑100 (MQ=255)
                                                  cAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCtggctggc                                               >  2:2280747/1‑100 (MQ=255)
                                                    aCCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGa                                             >  2:1177537/1‑100 (MQ=255)
                                                     ccATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGat                                            >  2:1007435/1‑100 (MQ=255)
                                                       aTGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGatat                                          <  1:2171823/100‑1 (MQ=255)
                                                       aTGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGatat                                          <  2:613794/100‑1 (MQ=255)
                                                        tGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGatat                                          <  1:71903/99‑1 (MQ=255)
                                                        tGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATg                                         <  2:2445870/100‑1 (MQ=255)
                                                        tGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATg                                         <  2:1788186/100‑1 (MQ=255)
                                                           ttCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCAt                                      <  2:2346172/100‑1 (MQ=255)
                                                            nccAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCAtt                                     >  1:2019739/2‑100 (MQ=255)
                                                              cAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTcc                                   >  2:2240037/1‑100 (MQ=255)
                                                                gagCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCa                                 >  1:1061289/1‑100 (MQ=255)
                                                                  gCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGa                               >  1:2422657/1‑100 (MQ=255)
                                                                   cGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGAc                              <  2:1729423/100‑1 (MQ=255)
                                                                     gCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGGCTc                            >  1:2596304/1‑100 (MQ=255)
                                                                       aaCACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGa                          <  1:2345707/100‑1 (MQ=255)
                                                                       aaCACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGa                          <  1:445048/100‑1 (MQ=255)
                                                                        acacTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAg                         <  1:27620/100‑1 (MQ=255)
                                                                        acacTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAg                         <  2:1967234/100‑1 (MQ=255)
                                                                         cacTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAg                         >  1:968679/1‑100 (MQ=255)
                                                                         cacTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGc                        <  1:1277211/100‑1 (MQ=255)
                                                                          ccgatctgtAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCa                       <  2:901106/94‑1 (MQ=255)
                                                                          acTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCa                       >  2:2291498/1‑100 (MQ=255)
                                                                          acTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCa                       <  2:1010434/100‑1 (MQ=255)
                                                                              cGTGTAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTg                   <  2:211033/100‑1 (MQ=255)
                                                                               gtgtAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGc                  >  2:2778890/1‑100 (MQ=255)
                                                                               gtgtAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGc                  <  1:907683/100‑1 (MQ=255)
                                                                                 gtAATTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCaagatcgg                >  1:901099/1‑94 (MQ=255)
                                                                                   aaTTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCttt               >  2:2326261/1‑99 (MQ=255)
                                                                                   aaTTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCtt                <  1:927139/98‑1 (MQ=255)
                                                                                    aTTTCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCttttt             <  2:2456046/100‑1 (MQ=255)
                                                                                     tttCTTCGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTa            <  2:865713/100‑1 (MQ=255)
                                                                                         ttcGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTg        >  2:2524195/1‑100 (MQ=255)
                                                                                         ttcGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTg        <  1:1234543/100‑1 (MQ=255)
                                                                                          tcGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGa       >  1:1860393/1‑100 (MQ=255)
                                                                                           cGTGTGGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGAt      <  2:913505/100‑1 (MQ=255)
                                                                                             tgtgGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATgc    >  1:1150856/1‑100 (MQ=255)
                                                                                              gtgGAAGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATgcg   >  2:1017334/1‑100 (MQ=255)
                                                                                                    |                                                                                               
GTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGC  >  NC_000913/1053411‑1053604

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