Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,145 A→T 15.2% 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.2% 374.5 / 55.5 145intergenic (+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 (60/63);  new base T (7/15);  total (67/78)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.18e-01

Only 100 of 157 total aligned reads displayed.

CTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCAT  >  NC_000913/2286046‑2286244
                                                                                                   |                                                                                                   
cTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATtata                                                                                                     <  2:419059/100‑1 (MQ=255)
 ttACGTGGTGAAAGAGTGAAAGATGAAGAACCACAGTTAAGTTTGTTATTGCAAGTGTTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATtata                                                                                                     <  2:1879215/99‑1 (MQ=255)
  tACGTCGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGATAAGCGTTTTATCGCTAACTGATTAATTATaaa                                                                                                   <  2:10171/100‑1 (MQ=255)
    cGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATc                                                                                                 >  1:709466/1‑100 (MQ=255)
    cGTGGTGAAAGAGTGAAAGAAAAAAAACCACAGTTAAGTTTGTTATTGCAAGTACTAACAGACGAGAAACGTTTTATCGCTAACTGATTAATTATAAATc                                                                                                 >  2:588069/1‑100 (MQ=255)
     gtggtgAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATc                                                                                                 <  1:3765244/99‑1 (MQ=255)
     gtggtgAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCa                                                                                                >  2:2119938/1‑100 (MQ=255)
      tggtgAAAGAGTGAAAGATGGAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAg                                                                                               >  1:3508549/1‑100 (MQ=255)
       ggtgAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGGAAGTGTTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAg                                                                                               <  2:3911130/99‑1 (MQ=255)
       ggtgAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGt                                                                                              <  2:3746433/100‑1 (MQ=255)
        gtgaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGtt                                                                                             >  2:369545/1‑100 (MQ=255)
         tgaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATGAATTATAAATAAGTTa                                                                                            >  2:2419586/1‑100 (MQ=255)
          gaaagagtgaaagaTGAAAAACCACGGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAg                                                                                           <  2:4004307/100‑1 (MQ=38)
          gaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAg                                                                                           <  2:3249150/100‑1 (MQ=38)
          gaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAg                                                                                           >  1:3942401/1‑100 (MQ=255)
          gaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAg                                                                                           <  1:1347498/100‑1 (MQ=255)
           aaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGc                                                                                          <  1:2571957/100‑1 (MQ=38)
           aaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGc                                                                                          <  1:3985085/100‑1 (MQ=38)
           aaagagtgaaagaTGAAAAAACACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGc                                                                                          <  1:106508/100‑1 (MQ=38)
           aaagaaaaaaagaAAAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAg                                                                                           <  2:3943174/91‑1 (MQ=255)
           aaaaagtgaaaaaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGCAGCGTTTTATCGCTAACTGATTAATTATAAATc                                                                                                 <  2:2729086/93‑1 (MQ=255)
            aagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGCg                                                                                         <  1:655287/100‑1 (MQ=38)
            aagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGCg                                                                                         <  2:823996/100‑1 (MQ=38)
            aagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCGAACTGATTTATAATTAATCAGTTAGc                                                                                          <  2:682403/99‑1 (MQ=38)
             agagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGCGa                                                                                        <  1:388366/100‑1 (MQ=38)
             agagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGAAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGa                                                                                        <  2:2616177/100‑1 (MQ=255)
              gagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGaa                                                                                       >  2:658384/1‑100 (MQ=255)
              gagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGaa                                                                                       >  2:662364/1‑100 (MQ=255)
                 tgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAata                                                                                    <  1:2357372/100‑1 (MQ=255)
                 tgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAata                                                                                    <  2:2892876/100‑1 (MQ=255)
                    tctaTGAAAACCCACAGTTAAGTTTGTTATTGCAAGTGTTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCTAAATATCt                                                                                 <  2:2828922/97‑1 (MQ=255)
                    aagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCt                                                                                 <  2:3418728/100‑1 (MQ=255)
                    aagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCt                                                                                 >  1:984465/1‑100 (MQ=255)
                      gaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTa                                                                               <  1:3221538/100‑1 (MQ=255)
                       aTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTAc                                                                              <  1:1235402/100‑1 (MQ=255)
                        tGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACt                                                                             >  2:3119598/1‑100 (MQ=255)
                        tGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACt                                                                             >  1:1335607/1‑100 (MQ=255)
                           aaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGc                                                                          <  1:2936693/100‑1 (MQ=255)
                           aaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGc                                                                          <  2:3104027/100‑1 (MQ=255)
                           aaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGc                                                                          <  1:2012104/100‑1 (MQ=255)
                              aCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAAt                                                                       >  1:3471825/1‑100 (MQ=255)
                                  caGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCgat                                                                   >  2:3792807/1‑98 (MQ=255)
                                  caGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGgt                                                                   <  1:3967289/100‑1 (MQ=255)
                                   aGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGgtg                                                                  <  1:1875419/100‑1 (MQ=255)
                                   aGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGgtg                                                                  <  2:1880344/100‑1 (MQ=255)
                                   aGTTAAGTTTGTGATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGgtg                                                                  <  2:1076718/100‑1 (MQ=255)
                                    gTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGgtgt                                                                 <  1:2683575/100‑1 (MQ=255)
                                     ttAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGtgtg                                                                <  1:1103571/100‑1 (MQ=255)
                                      tAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTgg                                                               <  1:1759044/100‑1 (MQ=255)
                                          tttGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCTCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGaaaa                                                           >  2:3392828/1‑100 (MQ=255)
                                           ttGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAAc                                                          <  2:1879108/100‑1 (MQ=255)
                                             gTTATTGCAAGTGCTGACAGACGAGCAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACgg                                                        >  2:2661452/1‑100 (MQ=255)
                                              ttATTGGAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACgg                                                        <  2:2263181/99‑1 (MQ=255)
                                                aTTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgtag                                                     <  1:106285/100‑1 (MQ=255)
                                                 ttGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgtagt                                                    <  1:1970727/100‑1 (MQ=255)
                                                  tGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgtagta                                                   >  1:1070093/1‑100 (MQ=255)
                                                   gCAAGTGCTGACAGACGAGAAGCGCTTTATCGCTAACTGAATAATTATACATCAGTTAGCGAAATATCTTACTTGCAACCGGTATGGAAAACGgtagta                                                   >  2:2386650/1‑99 (MQ=255)
                                                    cAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTAtt                                                 <  1:2659503/100‑1 (MQ=255)
                                                      aGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTag                                               >  1:168651/1‑100 (MQ=255)
                                                        tGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTagca                                             >  1:2113935/1‑100 (MQ=255)
                                                         gCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTagcag                                            >  1:421541/1‑100 (MQ=255)
                                                          cTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTagcagc                                           >  2:3566560/1‑100 (MQ=255)
                                                           taaCAGACGAGAAACGTTTTATCGCTAACTGATTAATTATAACTCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGcc                                          >  1:2786629/3‑100 (MQ=255)
                                                           tGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGcc                                          >  2:1568002/1‑100 (MQ=255)
                                                            gACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCa                                         >  2:3375516/1‑100 (MQ=255)
                                                              cAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACg                                       <  2:281111/100‑1 (MQ=255)
                                                              cAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACg                                       <  2:1315265/100‑1 (MQ=255)
                                                              cAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACg                                       <  2:1578389/100‑1 (MQ=255)
                                                              cAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACg                                       <  1:2576874/100‑1 (MQ=255)
                                                                gACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAg                                     >  1:2821984/1‑100 (MQ=255)
                                                                gACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAg                                     >  1:3298527/1‑100 (MQ=255)
                                                                gACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAg                                     >  2:2588404/1‑100 (MQ=255)
                                                                  cGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTc                                   >  2:1154496/1‑100 (MQ=255)
                                                                   gagaAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCg                                  >  2:1222289/1‑100 (MQ=255)
                                                                     gaAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGc                                >  1:713799/1‑100 (MQ=255)
                                                                       aGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCAc                              >  1:3144143/1‑100 (MQ=255)
                                                                             atttCGCTAACTGATTTATAATTAATCAGTTAGCTAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAgcg                        >  2:667355/4‑100 (MQ=37)
                                                                              tATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAgcgc                       >  1:2129810/1‑100 (MQ=255)
                                                                              tATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTAGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAgcgc                       >  1:1970099/1‑100 (MQ=255)
                                                                                 cGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAg                     >  2:4040668/1‑99 (MQ=255)
                                                                                  gCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGcc                   >  2:1507076/1‑100 (MQ=38)
                                                                                  gCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGcc                   >  1:3843751/1‑100 (MQ=38)
                                                                                  gCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGcc                   >  1:443271/1‑100 (MQ=38)
                                                                                  gCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGcc                   >  1:1020812/1‑100 (MQ=255)
                                                                                  gCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGcc                   >  2:514093/1‑100 (MQ=255)
                                                                                   cTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCt                  >  2:3206767/1‑100 (MQ=255)
                                                                                    tAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTg                 <  1:2097910/100‑1 (MQ=255)
                                                                                       cTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCAACGAGTCGGCACGTAGCGCAGCCTGGt               >  1:2231898/1‑99 (MQ=255)
                                                                                        tGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAg             >  2:3528620/1‑100 (MQ=255)
                                                                                         gATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAgc            <  2:2915265/100‑1 (MQ=255)
                                                                                          attaattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAgcg           >  1:2435915/1‑100 (MQ=255)
                                                                                           ttaattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAgcgc          <  1:658390/100‑1 (MQ=255)
                                                                                           ttaattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAgcgc          >  1:1928422/1‑100 (MQ=255)
                                                                                           ctaattaTAAATCAGTTAGCGAAATATATTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCAAGAGTCGGCACGTAGCGCAGCCTGGTAgcgc          >  2:1896985/2‑100 (MQ=255)
                                                                                             aattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCAc        >  1:1442870/1‑100 (MQ=255)
                                                                                             aattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCAc        <  1:3044801/100‑1 (MQ=255)
                                                                                               ttaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCg      <  1:2129783/100‑1 (MQ=255)
                                                                                                 ataAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTc    >  1:3754759/1‑100 (MQ=255)
                                                                                                  taAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCa   <  1:2849736/100‑1 (MQ=255)
                                                                                                   aaaTCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCAt  >  1:780265/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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: