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)
      gAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCg                                                                                        <  1:176990/100‑1 (MQ=255)
       aaaaTGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGcccgctc                                                                                          >  2:662488/1‑94 (MQ=255)
       aaaaCGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGt                                                                                       <  2:369020/100‑1 (MQ=255)
         aaTGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTcc                                                                                     <  2:1898594/100‑1 (MQ=255)
          aTGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTcc                                                                                     <  2:2237739/99‑1 (MQ=255)
            gACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCctgc                                                                                  <  2:197694/100‑1 (MQ=255)
            gACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCctgc                                                                                  <  1:782711/100‑1 (MQ=255)
            atctaTCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCctgc                                                                                  <  1:847826/96‑1 (MQ=255)
             aCGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCctgct                                                                                 <  1:1155729/100‑1 (MQ=255)
              cGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTGTGCTGTACCTAACACGGGAGCTTAACCCCGCTCGTc                                                                                      >  2:2427180/1‑94 (MQ=255)
               gATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGa                                                                               >  1:1218281/1‑100 (MQ=255)
               gATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGa                                                                               >  2:134143/1‑100 (MQ=255)
                aTCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCctgcagat                                                                              >  2:847819/1‑96 (MQ=255)
                aTCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAg                                                                              >  2:1275504/1‑100 (MQ=255)
                aTCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAg                                                                              >  2:987649/1‑100 (MQ=255)
                aTCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAg                                                                              >  1:1550472/1‑100 (MQ=255)
                 tCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGt                                                                             >  1:2207850/1‑100 (MQ=255)
                 tCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGt                                                                             >  2:1099248/1‑100 (MQ=255)
                  cTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTa                                                                            >  2:443945/1‑100 (MQ=255)
                   tGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTAt                                                                           >  2:996764/1‑100 (MQ=255)
                    gagaAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATc                                                                          >  1:176455/1‑100 (MQ=255)
                      gaAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCaa                                                                        <  1:434488/100‑1 (MQ=255)
                      gaAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCaa                                                                        <  1:1661730/100‑1 (MQ=255)
                      gaAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCaa                                                                        >  1:677441/1‑100 (MQ=255)
                         gCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAAt                                                                     >  1:925095/1‑100 (MQ=255)
                          cGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATg                                                                    >  2:823765/1‑100 (MQ=255)
                          cGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATg                                                                    >  1:1474570/1‑100 (MQ=255)
                           gATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATgg                                                                   <  1:749570/100‑1 (MQ=255)
                            aTAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGc                                                                  <  1:1050215/100‑1 (MQ=255)
                            aTAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGc                                                                  <  1:470653/100‑1 (MQ=255)
                            aTAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGc                                                                  <  1:1868153/100‑1 (MQ=255)
                             tAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCt                                                                 >  2:1105404/1‑100 (MQ=255)
                             tAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAAAGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCt                                                                 >  2:1797170/1‑100 (MQ=255)
                              aaGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTa                                                                <  1:277770/100‑1 (MQ=255)
                               aGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTAc                                                               >  2:2410895/1‑100 (MQ=255)
                                gcggcgAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACt                                                              >  2:774582/1‑100 (MQ=255)
                                gcggcgAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACt                                                              >  1:434138/1‑100 (MQ=255)
                                 cggcgAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTg                                                             >  1:331033/1‑100 (MQ=255)
                                       gcgaAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCa                                                        >  2:2031582/1‑99 (MQ=255)
                                       gcgaAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAg                                                       >  2:1191182/1‑100 (MQ=255)
                                         gaAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGt                                                     >  1:2365665/1‑100 (MQ=255)
                                         gaAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGt                                                     >  1:1406886/1‑100 (MQ=255)
                                          aaaTCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGt                                                     <  2:751374/99‑1 (MQ=255)
                                             tCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACa                                                 <  2:1640548/100‑1 (MQ=255)
                                             tCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACa                                                 <  2:1530207/100‑1 (MQ=255)
                                               gTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATc                                               <  1:507830/100‑1 (MQ=255)
                                               gTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATc                                               <  2:469062/100‑1 (MQ=255)
                                                tGCTGTTGTTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACa                                                 <  2:952822/98‑1 (MQ=255)
                                                 gCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCaa                                             <  1:1199724/100‑1 (MQ=255)
                                                  cTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCaaa                                            >  1:2072638/1‑100 (MQ=255)
                                                   tGTTGGTGGGATCCACTTTTTTGCTGTAATTAAAACGGAAGCTTAAGCCCGCTCGTCATGCTGAGTATCAAAATGGCTACTTGCCCAGGTGACATCAAAg                                           <  2:634018/100‑1 (MQ=255)
                                                      tggtggCATTTACTTGTTTGCTCTAATTAAAACGGTAGCTTAAGCCCGCTCCTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGAc                                        <  2:922109/100‑1 (MQ=255)
                                                      tggtggCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGAc                                        <  2:419072/100‑1 (MQ=255)
                                                      tggtggCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTAGCAGAATGGCGACTGGCCCAGGTGACATCAAAGGAc                                        <  1:2427194/100‑1 (MQ=255)
                                                       ggtggCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACa                                       >  1:25582/1‑100 (MQ=255)
                                                       ggtggCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACa                                       >  2:1944982/1‑100 (MQ=255)
                                                          ggCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGt                                     >  2:2053934/1‑99 (MQ=255)
                                                              tCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGcc                                <  1:2471824/100‑1 (MQ=255)
                                                               ccACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCa                               >  1:666939/1‑100 (MQ=255)
                                                                cACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCaa                              >  2:2367088/1‑100 (MQ=255)
                                                                cACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCaa                              >  2:1292993/1‑100 (MQ=255)
                                                                   ttGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGcc                           <  1:1835845/100‑1 (MQ=255)
                                                                   ttGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGcc                           <  1:1835812/100‑1 (MQ=255)
                                                                    tGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCca                          >  1:1408764/1‑100 (MQ=255)
                                                                      tttGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCcaca                        <  1:1159547/100‑1 (MQ=255)
                                                                       ttGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAg                       <  2:1576764/100‑1 (MQ=255)
                                                                       ttGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAg                       >  1:1450895/1‑100 (MQ=255)
                                                                        tGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGc                      <  2:752650/100‑1 (MQ=255)
                                                                          cTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCat                    <  1:42033/100‑1 (MQ=255)
                                                                           tGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCata                   >  2:2363483/1‑100 (MQ=255)
                                                                            gTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCatat                  >  1:772015/1‑100 (MQ=255)
                                                                             tAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATAtt                 <  2:1589356/100‑1 (MQ=255)
                                                                             tAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATAtt                 >  1:573440/1‑100 (MQ=255)
                                                                              aaTTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATAttt                <  1:630869/100‑1 (MQ=255)
                                                                               aTTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAACGGACAGTGCGCCAAGCCACAGCATAtttt               <  2:2365650/100‑1 (MQ=255)
                                                                               aTTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATAtttt               >  2:764008/1‑100 (MQ=255)
                                                                               aTTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATAtttt               >  1:871910/1‑100 (MQ=255)
                                                                               aTTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATAtttt               >  1:161314/1‑100 (MQ=255)
                                                                                 tAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGc             <  1:1843597/100‑1 (MQ=255)
                                                                                  aaaaCGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGcc            <  2:682268/100‑1 (MQ=255)
                                                                                    aaCGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCg          <  2:2258271/100‑1 (MQ=255)
                                                                                        gTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATAc      <  1:1511266/100‑1 (MQ=255)
                                                                                           gCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCa   >  2:914069/1‑100 (MQ=255)
                                                                                           gCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCa   >  1:2457121/1‑100 (MQ=255)
                                                                                            cTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCAt  >  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: