Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,145 A→T 15.5% 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.5% 246.0 / 38.1 97intergenic (+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 (27/54);  new base T (8/7);  total (36/61)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.69e-01

Only 100 of 103 total aligned reads displayed.

CTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCG  >  NC_000913/2286046‑2286240
                                                                                                   |                                                                                               
cTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATtata                                                                                                 <  1:1659220/100‑1 (MQ=255)
cTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATtata                                                                                                 <  2:306509/100‑1 (MQ=255)
  tACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATaaa                                                                                               <  1:1158474/100‑1 (MQ=255)
  tACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATaaa                                                                                               <  1:1336777/100‑1 (MQ=255)
    cGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATc                                                                                             <  2:1447749/100‑1 (MQ=255)
     gtggtgAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATc                                                                                             <  2:1083213/99‑1 (MQ=255)
     gtggtgAAAGAGTGAAAAATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATtata                                                                                                 <  2:1561490/95‑1 (MQ=255)
      tggtgAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAg                                                                                           <  1:70598/100‑1 (MQ=255)
       ggtgAAAGAGTGAAAGATGAAGAACCACGGTGAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAg                                                                                           <  2:2548068/99‑1 (MQ=255)
        gtgaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGtt                                                                                         >  1:2562673/1‑100 (MQ=255)
         tgaaatagttaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGtt                                                                                         <  2:162254/99‑1 (MQ=255)
          gaaagggtgcaataCGGAAGACCACCTTTATGTTTGTTATTGAAAGTGTTGACAGTCGACTGGGGTTGTATCGGTAACTGATTAATTATAAATCAGTTAg                                                                                       <  2:1769400/100‑1 (MQ=255)
           aaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGc                                                                                      <  2:2513437/100‑1 (MQ=38)
           aaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGc                                                                                      <  1:2604007/100‑1 (MQ=38)
           aaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGc                                                                                      <  2:1923608/100‑1 (MQ=255)
            aagagtgaaagaTGAAAACCCACAGTAAACTTTGTTACTGCAAGTGCTCACAGACGAGAAGCGTTTTATCGCGAACTGAGTTATAATTAATCAGtcaccg                                                                                     <  2:2136181/100‑6 (MQ=255)
            aagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGCg                                                                                     <  1:1093131/100‑1 (MQ=38)
            aaaagtaaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGc                                                                                      <  2:2347578/96‑1 (MQ=37)
               agtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGCGa                                                                                    <  2:2386921/98‑1 (MQ=38)
                gtgaaagaTGAAAAACCACAGTTAAGGTTTTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTGTCGCTAACTGAGTAATTATAAATCAGTTAGCGAAat                                                                                 >  1:545254/1‑100 (MQ=255)
                 tgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAata                                                                                <  1:2524122/100‑1 (MQ=255)
                  gaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAatat                                                                               >  1:2445688/1‑100 (MQ=255)
                    aaggTGAAAAACCACAGTTACGTCTGTTATTGCAAGTGCCGACAGACGCGCAGCGTTTTATGGCCAACTGCTTAAATATCAGTCCGCTAGCGAActa                                                                                >  2:2457069/1‑97 (MQ=255)
                     agaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCtt                                                                            >  2:209040/1‑100 (MQ=255)
                     agaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCtt                                                                            >  2:1841336/1‑100 (MQ=255)
                     agaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCtt                                                                            >  2:196514/1‑100 (MQ=255)
                       aTGAGAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTAc                                                                          <  2:363603/100‑1 (MQ=255)
                        tGAAAAATCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACt                                                                         >  1:1445681/1‑100 (MQ=255)
                         gAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACtt                                                                        >  2:1901614/1‑100 (MQ=255)
                          acaaaCCTCAGTTACGTTTGATATTGCAAGTGCTGACAGACGACAAGCATCTTATTGCAAACTGATTAATTATAAATCAGTTAGCGAAATATCTTGCt                                                                         >  2:532817/3‑98 (MQ=255)
                          aaaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTg                                                                       <  2:1382309/100‑1 (MQ=255)
                          aaaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTg                                                                       >  1:1127550/1‑100 (MQ=255)
                          aaaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGATAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTg                                                                       <  2:1578473/100‑1 (MQ=255)
                           aaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGc                                                                      <  2:2315123/100‑1 (MQ=255)
                           aaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGc                                                                      <  2:1238390/100‑1 (MQ=255)
                            aaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTg                                                                       <  2:234111/98‑1 (MQ=255)
                               ccACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATc                                                                  <  1:1057545/100‑1 (MQ=255)
                                  caGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGgt                                                               <  1:1249760/100‑1 (MQ=255)
                                  caGTGAAGTGTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTGATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCgg                                                                <  2:823277/99‑1 (MQ=255)
                                   aGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGgtg                                                              <  1:2376021/100‑1 (MQ=255)
                                     ttAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGtgtg                                                            <  1:361977/100‑1 (MQ=255)
                                       aaGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGa                                                          <  1:2470622/100‑1 (MQ=255)
                                           ttGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAAc                                                      <  1:1555161/100‑1 (MQ=255)
                                             gTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACgg                                                    <  1:1073773/100‑1 (MQ=255)
                                               tATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgta                                                  <  2:1192962/100‑1 (MQ=255)
                                                aTTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgtag                                                 <  2:625239/100‑1 (MQ=255)
                                                  tGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgtagta                                               >  1:1475078/1‑100 (MQ=255)
                                                  tGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgtagta                                               <  1:932248/100‑1 (MQ=255)
                                                   gCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTAt                                              <  1:713297/100‑1 (MQ=255)
                                                    cAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTAtt                                             <  2:191233/100‑1 (MQ=255)
                                                     aaGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTa                                            <  1:2693338/100‑1 (MQ=255)
                                                     aaGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTa                                            <  1:2303357/100‑1 (MQ=255)
                                                       gTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTagc                                          <  2:80008/100‑1 (MQ=255)
                                                       gTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTagc                                          <  2:357486/100‑1 (MQ=255)
                                                       gTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTagc                                          <  2:1768529/100‑1 (MQ=255)
                                                        tGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTagca                                         <  2:530777/100‑1 (MQ=255)
                                                           tGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGcc                                      >  1:618458/1‑100 (MQ=255)
                                                            gACAGACGAGAAGCGTTTTATCGCTAACTGATTAATGATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCa                                     >  1:1939974/1‑100 (MQ=255)
                                                              cAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGTAAAACGGTAGTATTAGCAGGCACg                                   <  2:1329782/100‑1 (MQ=255)
                                                              cAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACg                                   <  1:1217734/100‑1 (MQ=255)
                                                                gACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAg                                 <  2:269810/100‑1 (MQ=255)
                                                                gACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAg                                 >  1:696889/1‑100 (MQ=255)
                                                                 aCGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCTGTGTGGAAAACGGTAGTATTAGCAGCCACGAg                                 <  1:2609171/99‑1 (MQ=255)
                                                                 aCGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGt                                >  2:2178763/1‑100 (MQ=255)
                                                                  cGAGAAGCGTTTTATCGCTACCTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTc                               >  2:1675726/1‑100 (MQ=255)
                                                                  cGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTc                               >  1:1873418/1‑100 (MQ=255)
                                                                   gagaAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGTTGTCGAAAACGGTAGTATTAGCAGCCACGAGTCg                              <  2:953996/100‑1 (MQ=255)
                                                                     gaAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGc                            <  1:1018630/100‑1 (MQ=255)
                                                                     gaAGCGTTTTATCACTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGc                            >  1:1861503/1‑100 (MQ=255)
                                                                     gaAGCGTCTTACCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGACACGAGTCgg                             >  2:1444514/1‑99 (MQ=255)
                                                                      aaGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCa                           <  1:2409245/100‑1 (MQ=255)
                                                                        gCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGCCGGCACg                         >  1:814059/1‑100 (MQ=255)
                                                                         cGTTTTATCGCTAACTGATTAATTAGAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACTGTAGTATTAGCAGCCACGAGTCGGCACg                         <  2:980416/99‑1 (MQ=255)
                                                                          gTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTa                       >  1:898017/1‑100 (MQ=255)
                                                                           ttttATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAg                      <  2:2562655/100‑1 (MQ=255)
                                                                            tttATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAgc                     <  2:632097/100‑1 (MQ=255)
                                                                            tttATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAgc                     <  2:2101178/100‑1 (MQ=255)
                                                                             ttATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAgcg                    >  1:1549/1‑100 (MQ=255)
                                                                             ttATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAgcg                    <  1:2590261/100‑1 (MQ=255)
                                                                              tttCGCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAgcgc                   >  2:313673/3‑100 (MQ=38)
                                                                               aTCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCa                  <  2:2232682/100‑1 (MQ=255)
                                                                                tCGCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAg                 >  1:6854/1‑100 (MQ=38)
                                                                                tCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAg                 <  2:696887/100‑1 (MQ=255)
                                                                                  gCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCTGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAg                 >  2:2268651/1‑98 (MQ=38)
                                                                                  gCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGcc               >  1:655736/1‑100 (MQ=38)
                                                                                  gCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGcc               >  1:1722545/1‑100 (MQ=38)
                                                                                  gCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGAAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGcc               >  2:860768/1‑100 (MQ=38)
                                                                                  gCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGcc               >  1:656977/1‑100 (MQ=255)
                                                                                   cTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCt              >  1:309397/1‑100 (MQ=38)
                                                                                   cTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCt              >  1:479937/1‑100 (MQ=38)
                                                                                     aaCTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTgg            >  2:2351754/1‑100 (MQ=255)
                                                                                      aCTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGt           <  1:509600/100‑1 (MQ=255)
                                                                                       cTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTa          <  2:2218346/100‑1 (MQ=255)
                                                                                         gATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAgc        >  1:2213763/1‑100 (MQ=255)
                                                                                          attaattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAgcg       <  1:836718/100‑1 (MQ=255)
                                                                                             aattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCAc    <  2:1844738/100‑1 (MQ=255)
                                                                                             aattaCAAATCAGTTCGCGACATATCCCACTTGCACTCGGTGTGGAAAACCGTAGCATTCGCAGCCACGAGTCGGCAAGTAGCGCAGCCCGGTCgcgccc    <  2:2614524/100‑3 (MQ=255)
                                                                                              actataAATCAGTAAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCAcc   >  1:1866048/3‑100 (MQ=255)
                                                                                               ttaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCg  <  1:1233647/100‑1 (MQ=255)
                                                                                               ttaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCg  <  2:861720/100‑1 (MQ=255)
                                                                                                   |                                                                                               
CTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCG  >  NC_000913/2286046‑2286240

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: