Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 750,896 Δ1 bp 100% coding (1274/2448 nt) ybgQ ← putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913750,8930C.97.6% 372.5 / 0.7 84coding (1277/2448 nt)ybgQputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (43/39);  total (44/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCAT  >  NC_000913/750795‑750987
                                                                                                  |                                                                                              
aGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGcc                                                                                               <  1:1326973/100‑1 (MQ=255)
 gCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGccc                                                                                              >  1:2358381/1‑100 (MQ=255)
 gCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGccc                                                                                              >  2:1159539/1‑100 (MQ=255)
   tgtgAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGcccgc                                                                                            <  1:2231576/100‑3 (MQ=255)
   tgtgAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGcccgc                                                                                            >  1:1906137/1‑98 (MQ=255)
   tgtgAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGcccgc                                                                                            <  2:2390462/100‑3 (MQ=255)
    gtgAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGcccgct                                                                                           >  2:1843583/1‑97 (MQ=255)
    gtgAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGcccgct                                                                                           >  2:2471810/1‑97 (MQ=255)
     tgAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGcccgctc                                                                                          >  1:448405/1‑96 (MQ=255)
     tgAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGcccgctc                                                                                          <  1:955030/100‑5 (MQ=255)
      gAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCg                                                                                        <  1:176990/100‑1 (MQ=255)
       aaaaTGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGcccgctc                                                                                          >  2:662488/1‑94 (MQ=255)
       aaaaCGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGt                                                                                       <  2:369020/100‑1 (MQ=255)
         aaTGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTcc                                                                                     <  2:1898594/100‑1 (MQ=255)
          aTGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTcc                                                                                     <  2:2237739/99‑1 (MQ=255)
            gACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCctgc                                                                                  <  2:197694/100‑1 (MQ=255)
            gACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCctgc                                                                                  <  1:782711/100‑1 (MQ=255)
            atctaTCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCctgc                                                                                  <  1:847826/96‑1 (MQ=255)
             aCGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCctgct                                                                                 <  1:1155729/100‑1 (MQ=255)
              cGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTGTGCTGTACCTAACACGGGAGCTTAA‑CCCCGCTCGTc                                                                                      >  2:2427180/1‑94 (MQ=255)
               gATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGa                                                                               >  1:1218281/1‑100 (MQ=255)
               gATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGa                                                                               >  2:134143/1‑100 (MQ=255)
                aTCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCctgcagat                                                                              >  2:847819/1‑96 (MQ=255)
                aTCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAg                                                                              >  2:1275504/1‑100 (MQ=255)
                aTCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAg                                                                              >  2:987649/1‑100 (MQ=255)
                aTCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAg                                                                              >  1:1550472/1‑100 (MQ=255)
                 tCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGt                                                                             >  1:2207850/1‑100 (MQ=255)
                 tCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGt                                                                             >  2:1099248/1‑100 (MQ=255)
                  cTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTa                                                                            >  2:443945/1‑100 (MQ=255)
                   tGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTAt                                                                           >  2:996764/1‑100 (MQ=255)
                    gagaAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATc                                                                          >  1:176455/1‑100 (MQ=255)
                      gaAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCaa                                                                        <  1:434488/100‑1 (MQ=255)
                      gaAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCaa                                                                        <  1:1661730/100‑1 (MQ=255)
                      gaAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCaa                                                                        >  1:677441/1‑100 (MQ=255)
                         gCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAAt                                                                     >  1:925095/1‑100 (MQ=255)
                          cGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATg                                                                    >  2:823765/1‑100 (MQ=255)
                          cGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATg                                                                    >  1:1474570/1‑100 (MQ=255)
                           gATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATgg                                                                   <  1:749570/100‑1 (MQ=255)
                            aTAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGc                                                                  <  1:1050215/100‑1 (MQ=255)
                            aTAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGc                                                                  <  1:470653/100‑1 (MQ=255)
                            aTAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGc                                                                  <  1:1868153/100‑1 (MQ=255)
                             tAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCt                                                                 >  2:1105404/1‑100 (MQ=255)
                             tAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAAAGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCt                                                                 >  2:1797170/1‑100 (MQ=255)
                              aaGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTa                                                                <  1:277770/100‑1 (MQ=255)
                               aGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTAc                                                               >  2:2410895/1‑100 (MQ=255)
                                gcggcgAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACt                                                              >  2:774582/1‑100 (MQ=255)
                                gcggcgAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACt                                                              >  1:434138/1‑100 (MQ=255)
                                 cggcgAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTg                                                             >  1:331033/1‑100 (MQ=255)
                                       gcgaAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCa                                                        >  2:2031582/1‑99 (MQ=255)
                                       gcgaAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAg                                                       >  2:1191182/1‑100 (MQ=255)
                                         gaAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGt                                                     >  1:2365665/1‑100 (MQ=255)
                                         gaAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGt                                                     >  1:1406886/1‑100 (MQ=255)
                                          aaaTCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGt                                                     <  2:751374/99‑1 (MQ=255)
                                             tCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACa                                                 <  2:1640548/100‑1 (MQ=255)
                                             tCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACa                                                 <  2:1530207/100‑1 (MQ=255)
                                               gTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATc                                               <  1:507830/100‑1 (MQ=255)
                                               gTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATc                                               <  2:469062/100‑1 (MQ=255)
                                                tGCTGTTGTTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACa                                                 <  2:952822/98‑1 (MQ=255)
                                                 gCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCaa                                             <  1:1199724/100‑1 (MQ=255)
                                                  cTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCaaa                                            >  1:2072638/1‑100 (MQ=255)
                                                   tGTTGGTGGGATCCACTTTTTTGCTGTAATTAAAACGGAAGCTTAAG‑CCCGCTCGTCATGCTGAGTATCAAAATGGCTACTTGCCCAGGTGACATCAAAg                                           <  2:634018/100‑1 (MQ=255)
                                                      tggtggCATTTACTTGTTTGCTCTAATTAAAACGGTAGCTTAAG‑CCCGCTCCTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGAc                                        <  2:922109/100‑1 (MQ=255)
                                                      tggtggCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGAc                                        <  2:419072/100‑1 (MQ=255)
                                                      tggtggCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTAGCAGAATGGCGACTGGCCCAGGTGACATCAAAGGAc                                        <  1:2427194/100‑1 (MQ=255)
                                                       ggtggCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACa                                       >  1:25582/1‑100 (MQ=255)
                                                       ggtggCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACa                                       >  2:1944982/1‑100 (MQ=255)
                                                          ggCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGt                                     >  2:2053934/1‑99 (MQ=255)
                                                              tCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGcc                                <  1:2471824/100‑1 (MQ=255)
                                                               ccACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCa                               >  1:666939/1‑100 (MQ=255)
                                                                cACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCaa                              >  2:2367088/1‑100 (MQ=255)
                                                                cACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCaa                              >  2:1292993/1‑100 (MQ=255)
                                                                   ttGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGcc                           <  1:1835845/100‑1 (MQ=255)
                                                                   ttGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGcc                           <  1:1835812/100‑1 (MQ=255)
                                                                    tGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCca                          >  1:1408764/1‑100 (MQ=255)
                                                                      tttGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCcaca                        <  1:1159547/100‑1 (MQ=255)
                                                                       ttGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAg                       <  2:1576764/100‑1 (MQ=255)
                                                                       ttGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAg                       >  1:1450895/1‑100 (MQ=255)
                                                                        tGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGc                      <  2:752650/100‑1 (MQ=255)
                                                                          cTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCat                    <  1:42033/100‑1 (MQ=255)
                                                                           tGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCata                   >  2:2363483/1‑100 (MQ=255)
                                                                            gTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCatat                  >  1:772015/1‑100 (MQ=255)
                                                                             tAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATAtt                 <  2:1589356/100‑1 (MQ=255)
                                                                             tAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATAtt                 >  1:573440/1‑100 (MQ=255)
                                                                              aaTTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATAttt                <  1:630869/100‑1 (MQ=255)
                                                                               aTTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAACGGACAGTGCGCCAAGCCACAGCATAtttt               <  2:2365650/100‑1 (MQ=255)
                                                                               aTTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATAtttt               >  2:764008/1‑100 (MQ=255)
                                                                               aTTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATAtttt               >  1:871910/1‑100 (MQ=255)
                                                                               aTTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATAtttt               >  1:161314/1‑100 (MQ=255)
                                                                                 tAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGc             <  1:1843597/100‑1 (MQ=255)
                                                                                  aaaaCGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGcc            <  2:682268/100‑1 (MQ=255)
                                                                                    aaCGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCg          <  2:2258271/100‑1 (MQ=255)
                                                                                        gTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATAc      <  1:1511266/100‑1 (MQ=255)
                                                                                           gCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCa   >  2:914069/1‑100 (MQ=255)
                                                                                           gCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCa   >  1:2457121/1‑100 (MQ=255)
                                                                                            cTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCAt  >  2:877125/1‑100 (MQ=255)
                                                                                                  |                                                                                              
AGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCAT  >  NC_000913/750795‑750987

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