Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,145 A→T 15.6% intergenic (+9/‑66) yejM → / → proL essential inner membrane DUF3413 domain‑containing protein; lipid A production and membrane permeability factor/tRNA‑Pro

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,286,1450AT15.6% 312.2 / 41.8 125intergenic (+9/‑66)yejM/proLessential inner membrane DUF3413 domain‑containing protein; lipid A production and membrane permeability factor/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base A (51/54);  new base T (7/13);  total (58/67)
Fisher's exact test for biased strand distribution p-value = 3.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01

Only 100 of 138 total aligned reads displayed.

CTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCAT  >  NC_000913/2286046‑2286244
                                                                                                   |                                                                                                   
cTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATtata                                                                                                     >  2:70631/1‑100 (MQ=255)
cTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATtata                                                                                                     >  1:3326336/1‑100 (MQ=255)
  tACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATaaa                                                                                                   <  1:957767/100‑1 (MQ=255)
  tACGTGGTGAAAGAGTGAAACATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATaaa                                                                                                   >  2:2396197/1‑100 (MQ=255)
    cGTGTTGAAAGAGTGAAAGATGAAAATCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATc                                                                                                 <  1:3122881/100‑1 (MQ=255)
    cGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTACTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATc                                                                                                 <  1:575912/100‑1 (MQ=255)
      ttgtgaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAg                                                                                               <  1:422542/98‑1 (MQ=39)
      tggtgAAATAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAg                                                                                               <  2:1884931/100‑1 (MQ=255)
      tggtgAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAg                                                                                               <  2:976719/100‑1 (MQ=255)
      tggtgAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATAGCTAACTGATTAATTATAAATCAg                                                                                               >  1:1252222/1‑100 (MQ=255)
       ggtgAAAGAGTGGAGGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATc                                                                                                 <  2:1305012/97‑1 (MQ=255)
        gtgaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGtt                                                                                             <  2:655668/100‑1 (MQ=39)
        gtgaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAg                                                                                               <  2:2368015/98‑1 (MQ=255)
         tgaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGATGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTa                                                                                            <  1:777440/100‑1 (MQ=255)
           aaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGc                                                                                          <  1:2905873/100‑1 (MQ=38)
           aaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGc                                                                                          <  1:3502702/100‑1 (MQ=38)
           aaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGc                                                                                          <  2:3823069/100‑1 (MQ=38)
           aaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGc                                                                                          <  1:1216473/100‑1 (MQ=255)
           aaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCATTTTATCTCTAACTGATTAATTATAAATCAGTTAGc                                                                                          <  2:258818/100‑1 (MQ=255)
            aagagtgaaagaTGAAAAACAACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGc                                                                                          >  2:1722267/1‑99 (MQ=255)
             agattgaaagaTCAAACACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGCg                                                                                         <  2:3452109/99‑1 (MQ=37)
             agagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGCGa                                                                                        <  2:197427/100‑1 (MQ=38)
                 tgaaagaTGAAAATCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAata                                                                                    <  1:1483098/100‑1 (MQ=255)
                 tgaaagaTGAAAAACCACGGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAata                                                                                    <  2:1395784/100‑1 (MQ=255)
                  gaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAatat                                                                                   >  2:1451501/1‑100 (MQ=255)
                  gaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAatat                                                                                   >  2:1639205/1‑100 (MQ=255)
                        tGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACt                                                                             >  1:10366/1‑100 (MQ=255)
                        tGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAAAAGCGTTTTATCGCTAACTGAATAATTATAAATCAGTTAGCGAAATATCTTACt                                                                             >  1:3115998/1‑100 (MQ=255)
                          aaaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTg                                                                           >  2:1514616/1‑100 (MQ=255)
                           aaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACATACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGc                                                                          <  1:128475/100‑1 (MQ=255)
                            aaaCCACAGTTAAGTTTATTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCCAAATATCTTACTTGCa                                                                         >  2:2280012/1‑100 (MQ=255)
                              aCCACAGTTAAGTTTGTTATTGCAAGTGTTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAAt                                                                       <  1:1646921/100‑1 (MQ=255)
                               ccACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATc                                                                      >  2:117660/1‑100 (MQ=255)
                                   attgaaGTTGGTTATTACAAGTGCTGAGAGACGAGAAGCGTTTTATCGATAACTGATTAATTATAAATCAGTTATCGAAATATCTTACTTGCAATCGgtg                                                                  <  2:1850421/96‑1 (MQ=255)
                                    gTTAAGTTTGTTATTGCAAGTGCTGATAGACGAGTTGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCACTCGgtgt                                                                 <  1:2047545/100‑1 (MQ=255)
                                    gTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGgtgt                                                                 >  2:3577212/1‑100 (MQ=255)
                                    gTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGgtgt                                                                 <  2:2672815/100‑1 (MQ=255)
                                    gTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGgtgt                                                                 <  1:1870149/100‑1 (MQ=255)
                                     ttAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGtgtg                                                                <  1:3182007/100‑1 (MQ=255)
                                        aGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGaa                                                             >  2:3294294/1‑100 (MQ=255)
                                         gTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGaaa                                                            >  2:77512/1‑100 (MQ=255)
                                          tttGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGaaaa                                                           >  2:1339175/1‑100 (MQ=255)
                                             gTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACgg                                                        <  2:2443259/100‑1 (MQ=255)
                                               tATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgta                                                      >  2:1690827/1‑100 (MQ=255)
                                               tATTGCAAGTGCTGACAGACGAGAAGCCTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgt                                                       >  2:3900753/1‑99 (MQ=255)
                                                nttGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgtag                                                     >  1:1363989/2‑100 (MQ=255)
                                                aTTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgtag                                                     <  1:3535741/100‑1 (MQ=255)
                                                aTTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTGGCAATCGGTGTGGAAAACGgtag                                                     <  1:3375296/100‑1 (MQ=255)
                                                aTTGCAAGTGATGACAGACGAGAAGCGTTTTATCGCTAACTGATGAATTATAAAGCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACgg                                                        >  2:2775022/1‑97 (MQ=255)
                                                  tGCAAGTGATGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgtagta                                                   <  1:1916665/100‑1 (MQ=255)
                                                   gCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTAt                                                  >  1:2422275/1‑100 (MQ=255)
                                                   gCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTAt                                                  >  2:347409/1‑100 (MQ=255)
                                                   gCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTAt                                                  >  2:3797430/1‑100 (MQ=255)
                                                   gCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTAt                                                  >  1:3012828/1‑100 (MQ=255)
                                                    cAAGTCCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTAtt                                                 >  2:1768464/1‑100 (MQ=255)
                                                     aaGTGCTGACAGACGAAAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTATCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTa                                                <  2:1855105/100‑1 (MQ=255)
                                                       gTGCTGCCAGACAAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGTTGTGTAAAACGGTAGTATTagc                                              <  2:1093005/100‑1 (MQ=255)
                                                        tGCTGACAGACGATAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTagca                                             <  1:1806635/100‑1 (MQ=255)
                                                        tGCTGACAGACGAGAAGCGTTTTATCGCTAACTCATTAATTATAAATCAGTTAGCGAAATATCTTACTTACAATCGGAGTGGCAAATGGCAGTAt                                                  >  2:2575665/1‑95 (MQ=255)
                                                         gCTGACAGACGAGAAGCGTTTTATCGGTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTACAATCGGTGTGGAAAACGGTAGTATTagcag                                            >  1:2377685/1‑100 (MQ=255)
                                                         gCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTagca                                             >  1:2397770/1‑99 (MQ=255)
                                                           tGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGTAAAACGGTAGTATTAGCAGcc                                          <  2:474537/100‑1 (MQ=255)
                                                             aCAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCAc                                        <  1:2336134/100‑1 (MQ=255)
                                                              cAGACGAGAATCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCATTCAGTGTGGAAAACGGTAGTATTAGCAGCCACg                                       <  2:2569674/100‑1 (MQ=255)
                                                              cAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCAc                                        <  1:2695584/99‑1 (MQ=255)
                                                              cAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACg                                       >  1:2274412/1‑100 (MQ=255)
                                                              cAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACg                                       <  2:2805887/100‑1 (MQ=255)
                                                                gACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAg                                     >  1:1420887/1‑100 (MQ=255)
                                                                gACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAg                                     >  2:3932571/1‑100 (MQ=255)
                                                                gACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGATGTGGAAAACGGTAGTATTAGCAGCCACGAg                                     >  2:92218/1‑100 (MQ=255)
                                                                 aCGAGAAGCGTTTTATCGCTAACTGATTAATTATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAg                                     <  1:221445/99‑1 (MQ=255)
                                                                    agaAGCGTTTTATCGCTAACTTATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCg                                  <  1:1451516/99‑1 (MQ=255)
                                                                    agaAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCgg                                 >  2:497827/1‑100 (MQ=255)
                                                                     gaAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGc                                >  1:3554760/1‑100 (MQ=255)
                                                                      aaGCTTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCa                               <  2:3427379/100‑1 (MQ=255)
                                                                           ttttATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAg                          <  2:1758167/100‑1 (MQ=255)
                                                                             ttATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTAATTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAgcg                        <  1:569459/100‑1 (MQ=255)
                                                                              tATCGCTAACTGATTAATTATAAATCAGTTAGCGCAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAgcgc                       >  1:167004/1‑100 (MQ=255)
                                                                               ttCGCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCa                      >  1:2529357/2‑100 (MQ=38)
                                                                                  gCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGcc                   >  2:2610092/1‑100 (MQ=38)
                                                                                  gCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGAAGCCACGAGTCGGCACGTAGCGCAGcc                   >  2:3492221/1‑100 (MQ=38)
                                                                                  gCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGATGTGGAAAACGGTAATATTAGCAGCCACGAGTCGGCaa                              >  2:2848772/1‑88 (MQ=37)
                                                                                  gCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGcc                   >  1:1031938/1‑100 (MQ=255)
                                                                                   ccaaCTGATTTATAATTAATCAGTTAGCGAAGTATCTTACTTACAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGc                    >  1:2648216/3‑98 (MQ=37)
                                                                                   cTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAATAGGCACGTAGCGCAGcc                   >  1:1690520/1‑99 (MQ=37)
                                                                                   cTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCt                  <  2:2253458/100‑1 (MQ=255)
                                                                                      aCTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGt               <  2:1031927/100‑1 (MQ=255)
                                                                                        tGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAg             >  2:3276811/1‑100 (MQ=255)
                                                                                        tGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAg             >  2:2960635/1‑100 (MQ=255)
                                                                                        tGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAg             >  1:1071934/1‑100 (MQ=255)
                                                                                         gATTAATTATAAATCAGTTAGCGAATTATCTTACTTGCAACCGTTGTGGAAAACGGTATTATTAGCAGCTACGAGTCTTCACGTAGTGCAGCCTGGTAgc            <  2:2420254/100‑1 (MQ=255)
                                                                                             aattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCAc        <  2:2827946/100‑1 (MQ=255)
                                                                                             aattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCAc        <  1:2221749/100‑1 (MQ=255)
                                                                                             aattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCAc        >  2:2512166/1‑100 (MQ=255)
                                                                                              attaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAATCGGCACGTAGCGCAGCCTGGTAgcgc          >  2:3472286/1‑97 (MQ=255)
                                                                                                tataAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGt     >  2:430765/1‑100 (MQ=255)
                                                                                                 ataAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTc    >  1:2936361/1‑100 (MQ=255)
                                                                                                 ataAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTc    <  1:1797603/100‑1 (MQ=255)
                                                                                                  taAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCa   <  2:419240/100‑1 (MQ=255)
                                                                                                  taAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCa   >  1:1402033/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCAT  >  NC_000913/2286046‑2286244

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