Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,142 T→A 16.1% intergenic (+6/‑69) 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,1420TA16.1% 240.3 / 33.7 93intergenic (+6/‑69)yejM/proLessential inner membrane DUF3413 domain‑containing protein; lipid A production and membrane permeability factor/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base T (26/51);  new base A (8/7);  total (35/58)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

CAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCG  >  NC_000913/2286043‑2286240
                                                                                                   |                                                                                                  
cAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGattaatt                                                                                                    >  1:1192967/1‑100 (MQ=255)
cAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGattaatt                                                                                                    >  1:357486/1‑100 (MQ=255)
cAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGattaatt                                                                                                    >  1:1342978/1‑100 (MQ=255)
 aaCTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGattaatta                                                                                                   <  2:824652/100‑1 (MQ=255)
 aaCTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGattaatta                                                                                                   <  2:2015590/100‑1 (MQ=255)
 aaCTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGattaatta                                                                                                   <  2:2203561/100‑1 (MQ=255)
   cTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATtata                                                                                                 <  2:306509/100‑1 (MQ=255)
   cTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATtata                                                                                                 <  1:1659220/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:1841336/1‑100 (MQ=255)
                        agaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCtt                                                                            >  2:196514/1‑100 (MQ=255)
                        agaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCtt                                                                            >  2:209040/1‑100 (MQ=255)
                          aTGAGAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTAc                                                                          <  2:363603/100‑1 (MQ=255)
                            gAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACtt                                                                        >  2:1901614/1‑100 (MQ=255)
                             acaaaCCTCAGTTACGTTTGATATTGCAAGTGCTGACAGACGACAAGCATCTTATTGCAAACTGATTAATTATAAATCAGTTAGCGAAATATCTTGCt                                                                         >  2:532817/3‑98 (MQ=255)
                             aaaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTg                                                                       >  1:1127550/1‑100 (MQ=255)
                             aaaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTg                                                                       <  2:1382309/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)
                                          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:932248/100‑1 (MQ=255)
                                                     tGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgtagta                                               >  1:1475078/1‑100 (MQ=255)
                                                      gCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTAt                                              <  1:713297/100‑1 (MQ=255)
                                                      gCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTAt                                              >  1:415212/1‑100 (MQ=255)
                                                       cAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTAtt                                             <  2:191233/100‑1 (MQ=255)
                                                        aaGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTa                                            <  1:2303357/100‑1 (MQ=255)
                                                        aaGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTa                                            <  1:2693338/100‑1 (MQ=255)
                                                          gTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTagc                                          <  2:357486/100‑1 (MQ=255)
                                                          gTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTagc                                          <  2:1768529/100‑1 (MQ=255)
                                                          gTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTagc                                          <  2:80008/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)
                                                                            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:1722545/1‑100 (MQ=38)
                                                                                     gCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGcc               >  1:655736/1‑100 (MQ=38)
                                                                                     gCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGAAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGcc               >  2:860768/1‑100 (MQ=38)
                                                                                     gCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGcc               >  1:656977/1‑100 (MQ=255)
                                                                                      cTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCt              >  1:479937/1‑100 (MQ=38)
                                                                                      cTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCt              >  1:309397/1‑100 (MQ=38)
                                                                                         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)
                                                                                                   |                                                                                                  
CAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCG  >  NC_000913/2286043‑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: