Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,139 A→T 15.6% intergenic (+3/‑72) 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,1390AT15.6% 249.4 / 37.7 96intergenic (+3/‑72)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 (35/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.90e-01

Only 100 of 113 total aligned reads displayed.

TACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCAC  >  NC_000913/2286041‑2286238
                                                                                                  |                                                                                                   
tACATCTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGattaa                                                                                                    <  2:411607/100‑1 (MQ=255)
tACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGattaa                                                                                                    >  2:1726617/1‑100 (MQ=255)
tACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGattaa                                                                                                    >  2:2273127/1‑100 (MQ=255)
tACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGattaa                                                                                                    <  2:1376287/100‑1 (MQ=255)
tACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGattaa                                                                                                    >  2:2326665/1‑100 (MQ=255)
tACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTATCGCTAACTGattaat                                                                                                   >  1:2364782/1‑100 (MQ=255)
 aCAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGattaat                                                                                                   <  1:2008146/100‑1 (MQ=255)
  cAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGattaatt                                                                                                  >  1:357486/1‑100 (MQ=255)
  cAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGattaatt                                                                                                  >  1:1192967/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                                                                                               <  1:1659220/100‑1 (MQ=255)
     cTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATtata                                                                                               <  2:306509/100‑1 (MQ=255)
       tACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATaaa                                                                                             <  1:1336777/100‑1 (MQ=255)
       tACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATaaa                                                                                             <  1:1158474/100‑1 (MQ=255)
         cGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATc                                                                                           <  2:1447749/100‑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                                                                                    <  1:2604007/100‑1 (MQ=38)
                aaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGc                                                                                    <  2:2513437/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)
                         aaggTGAAAAACCACAGTTACGTCTGTTATTGCAAGTGCCGACAGACGCGCAGCGTTTTATGGCCAACTGCTTAAATATCAGTCCGCTAGCGAActa                                                                              >  2:2457069/1‑97 (MQ=255)
                          agaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCtt                                                                          >  2:196514/1‑100 (MQ=255)
                          agaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCtt                                                                          >  2:1841336/1‑100 (MQ=255)
                          agaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCtt                                                                          >  2:209040/1‑100 (MQ=255)
                            aTGAGAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTAc                                                                        <  2:363603/100‑1 (MQ=255)
                             tGAAAAATCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACt                                                                       >  1:1445681/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)
                               aaaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGATAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTg                                                                     <  2:1578473/100‑1 (MQ=255)
                                aaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGc                                                                    <  2:1238390/100‑1 (MQ=255)
                                aaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGc                                                                    <  2:2315123/100‑1 (MQ=255)
                                    ccACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATc                                                                <  1:1057545/100‑1 (MQ=255)
                                     cacaGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCg                                                               <  2:2037388/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: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:1768529/100‑1 (MQ=255)
                                                            gTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTagc                                        <  2:80008/100‑1 (MQ=255)
                                                            gTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTagc                                        <  2:357486/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                               >  1:696889/1‑100 (MQ=255)
                                                                      aCGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCTGTGTGGAAAACGGTAGTATTAGCAGCCACGAg                               <  1:2609171/99‑1 (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)
                                                                          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:2101178/100‑1 (MQ=255)
                                                                                  ttATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAgcg                  <  1:2590261/100‑1 (MQ=255)
                                                                                  ttATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAgcg                  >  1:1549/1‑100 (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)
                                                                                        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)
                                                                                                  |                                                                                                   
TACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCAC  >  NC_000913/2286041‑2286238

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: