Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,500,814 Δ1 bp 100% intergenic (‑23/‑446) yjgX ← / → yjgZ pseudogene fragment/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,500,8110T.94.2% 350.7 / 13.9 86intergenic (‑20/‑449)yjgX/yjgZpseudogene fragment/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base . (38/43);  total (40/46)
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

TCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTGCAAGAGGGGG  >  NC_000913/4500713‑4500906
                                                                                                  |                                                                                               
tCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGtt                                                                                                <  1:978838/100‑1 (MQ=255)
tCCGCCAGAAAATCACCACGAATACCAGAAACGCGGCACTGACCAGAAAATGATCAACAGGAACCACGGGGCGGTTCACTTGGCTTTCTCCTGTAGTGtt                                                                                                <  2:2023373/100‑1 (MQ=255)
 ccGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGttt                                                                                               >  2:2177616/1‑100 (MQ=255)
  cGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGtttt                                                                                              >  2:1850027/1‑100 (MQ=255)
  cGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGtttc                                                                                              <  1:1914356/100‑2 (MQ=255)
  cGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGtttc                                                                                              <  1:1824221/100‑2 (MQ=255)
  cGCCAGAAAATCACCACGAATACCAGAATCGCGCCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCAGTAGTGtttc                                                                                              >  2:1685820/1‑99 (MQ=255)
   gCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGtttct                                                                                             >  1:1813271/1‑98 (MQ=255)
    ccAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGtttctg                                                                                            <  2:1903452/100‑4 (MQ=255)
    ccAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGtttctg                                                                                            >  2:2554669/1‑97 (MQ=255)
     cAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGtttctgt                                                                                           >  2:791289/1‑96 (MQ=255)
     cAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGtttctgt                                                                                           >  1:1960431/1‑96 (MQ=255)
       gAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAg                                                                                        <  2:1930142/100‑1 (MQ=255)
        aaaaTCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAgg                                                                                       <  1:888942/100‑1 (MQ=255)
        aaaaTCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAgg                                                                                       <  1:1949040/100‑1 (MQ=255)
         aaaTCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGt                                                                                      <  1:50728/100‑1 (MQ=255)
            tCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCt                                                                                   <  1:1117286/100‑1 (MQ=255)
            tCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCt                                                                                   <  1:1255829/100‑1 (MQ=255)
             caccacGAATACCAGAATCGCGGCACTTACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCt                                                                                   <  1:1079771/99‑1 (MQ=255)
             caccacGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTg                                                                                  >  1:2182102/1‑100 (MQ=255)
              accacGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGa                                                                                 >  2:1324096/1‑100 (MQ=255)
              accacGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGa                                                                                 >  1:1167069/1‑100 (MQ=255)
              accacGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGa                                                                                 >  2:386336/1‑100 (MQ=255)
               ccacGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGaa                                                                                >  1:914408/1‑100 (MQ=255)
                 acGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAAcc                                                                              >  1:2276326/1‑100 (MQ=255)
                        ccAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTTCTGTAGGTTCTGAACCGGTTCt                                                                        <  1:1850041/100‑1 (MQ=255)
                        ccAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTTCTGTAGGTTCTGAACCGGTTCt                                                                        >  1:1760890/1‑100 (MQ=255)
                        ccAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTa                                                                       >  1:1286109/1‑100 (MQ=255)
                        ccAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTa                                                                       >  1:2465088/1‑100 (MQ=255)
                            aaTCGCGGCACTGACCAGAAAATGATAATCAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCgg                                                                   <  2:1505701/100‑1 (MQ=255)
                            aaTCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCgg                                                                   <  1:970046/100‑1 (MQ=255)
                             aTCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCggg                                                                  <  2:1286101/100‑1 (MQ=255)
                                gcgGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGtt                                                               >  2:177805/1‑100 (MQ=255)
                                gcgGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGtt                                                               >  2:870628/1‑100 (MQ=255)
                                 cgGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTtg                                                              <  2:87513/100‑1 (MQ=255)
                                 cgGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACCTTGCTTTCTCCTGTAGTGCTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTtg                                                              >  1:2232019/1‑100 (MQ=255)
                                   gCACTGACCAGAAAATGATAAACAGGAAACACGGAGCGGTTCACTATGCTTTCTCCTGTAGTGTTTCTGTAGGTACTCAACCGGTTATAGCGGGGTtg                                                              >  2:203988/1‑98 (MQ=255)
                                     aCTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACGGGTTCTAGCGGGGTTGTGAg                                                          <  2:1569102/100‑1 (MQ=255)
                                     aCTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAg                                                          >  2:2278972/1‑100 (MQ=255)
                                     aCTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAg                                                          <  1:395587/100‑1 (MQ=255)
                                     aCTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACAGGTTCTAGCGGGGTTGTGAg                                                          <  2:1569118/100‑1 (MQ=255)
                                      cTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTTCTGTAGGTTCTGAGCCGGTTCTAGCGGGGTTGTGAg                                                          <  1:289076/100‑1 (MQ=255)
                                      cTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAg                                                          <  2:1686449/99‑1 (MQ=255)
                                       tGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCGGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGc                                                        <  2:225235/100‑1 (MQ=255)
                                           cAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGAt                                                    <  1:547333/100‑1 (MQ=255)
                                           cAGAAAAGGATAAACAGGAGACACGGGGCGGTTCTCTTTGCTTTCTCCTGAAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGCTTTGAGGCGGAt                                                    >  2:2343263/1‑100 (MQ=255)
                                            aGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATg                                                   >  1:2332758/1‑100 (MQ=255)
                                            aGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATg                                                   <  2:1673380/100‑1 (MQ=255)
                                            aGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATg                                                   <  1:716616/100‑1 (MQ=255)
                                            aGAAAATGATAAACAGGAAACACGGGGCGGTTCAATTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATg                                                   >  2:1935614/1‑100 (MQ=255)
                                             gAAAATGATGAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATgg                                                  <  1:907829/100‑1 (MQ=255)
                                             gAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATgg                                                  <  1:928198/100‑1 (MQ=255)
                                              aaaaTGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGc                                                 <  1:849683/100‑1 (MQ=255)
                                               aaaTGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGc                                                 <  2:1760877/100‑1 (MQ=255)
                                                 aTGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAAt                                              >  1:440418/1‑100 (MQ=255)
                                                    aTAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGc                                           <  2:915939/100‑1 (MQ=255)
                                                    aTAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGc                                           >  2:1003324/1‑100 (MQ=255)
                                                     tAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGcc                                          >  1:1484165/1‑100 (MQ=255)
                                                       atctagAAACACGGGGCGGTTCAGTTTGCTTTGTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCgg                                        <  1:803066/95‑1 (MQ=255)
                                                       aaCAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCgg                                        <  2:2543057/100‑1 (MQ=255)
                                                       aaCAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCgg                                        >  2:1635750/1‑100 (MQ=255)
                                                        aCAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCggg                                       <  1:1204827/100‑1 (MQ=255)
                                                         cAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGc                                      >  1:369328/1‑100 (MQ=255)
                                                         cAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGc                                      <  1:496472/100‑1 (MQ=255)
                                                            gAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGctct                                   <  2:1223470/100‑1 (MQ=255)
                                                             aaaCACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCtt                                  >  1:1871961/1‑100 (MQ=255)
                                                              aaCACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTc                                 <  2:589668/100‑1 (MQ=255)
                                                                cacGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCt                               >  2:2097199/1‑100 (MQ=255)
                                                                cacGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCt                               >  1:1198638/1‑100 (MQ=255)
                                                                cacGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCt                               >  2:1384965/1‑100 (MQ=255)
                                                                  cGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGc                             >  2:1451785/1‑100 (MQ=255)
                                                                   ggggCGGTTCACTTTGCTTTCTCCTGTAGTGTTTGTGTAGGTTCTGAACCGGTTCTAGCGGTGTTGTGAGGCTGATGGCAATCGCGGGGCTCTTCCTGCa                            >  2:142486/1‑100 (MQ=255)
                                                                    gggCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCAt                           >  2:1499527/1‑100 (MQ=255)
                                                                     ggCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATg                          <  2:1038257/100‑1 (MQ=255)
                                                                     ggCGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATg                          <  1:671788/100‑1 (MQ=255)
                                                                       cGGTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCCTGtt                        >  1:1567705/1‑100 (MQ=255)
                                                                        ggTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTa                       >  2:783852/1‑100 (MQ=255)
                                                                        ggTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTa                       >  1:272272/1‑100 (MQ=255)
                                                                        ggTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTa                       <  2:1532995/100‑1 (MQ=255)
                                                                        ggTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTa                       <  2:1775446/100‑1 (MQ=255)
                                                                         gatctcTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTac                      <  1:2027299/95‑1 (MQ=255)
                                                                         gTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTTTTCCTGCATGTTac                      <  2:1310077/100‑1 (MQ=255)
                                                                         gTTCACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTac                      <  1:663181/100‑1 (MQ=255)
                                                                           tCACTTTGCTTTCTCCTGTGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACt                   <  1:1994258/100‑1 (MQ=255)
                                                                            cACTTTGCTTTCTCCTGTCGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTacac                    >  2:344556/1‑99 (MQ=255)
                                                                            cACTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACt                   >  2:624742/1‑100 (MQ=255)
                                                                             aCTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTg                  >  1:1173911/1‑100 (MQ=255)
                                                                              cTTTGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTacagatc                 >  2:2027285/1‑96 (MQ=255)
                                                                               tttGCTTTCTCCTGTAGTGTTTCTGTGGGTTCTGAACCGGTTCTATCGGGGTTGTGAGGCTGATGGCAATCGCCGGACTCTTCCTGCATGTTACACTGGt                <  2:1377237/100‑1 (MQ=255)
                                                                                ttGCTTTTTCCTGTAGTGTTTCTGTAGGTTCTGAAACGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTa               <  2:29061/100‑1 (MQ=255)
                                                                                 tGCTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTAc              <  1:581294/100‑1 (MQ=255)
                                                                                   cTTTCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTg            >  1:2466709/1‑100 (MQ=255)
                                                                                    tttCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTGc           <  2:1762317/100‑1 (MQ=255)
                                                                                    tttCTCCTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTGc           >  1:2491002/1‑100 (MQ=255)
                                                                                         ccTGTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTGCAagag      >  2:783899/1‑100 (MQ=255)
                                                                                            gTAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTGCAAGAgggg   <  1:199608/100‑1 (MQ=255)
                                                                                             tAGTGTTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTGCAAGAgggg   >  1:710840/1‑100 (MQ=255)
                                                                                             tAGTGTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTGCAAGAggggg  <  2:1173903/100‑1 (MQ=255)
                                                                                                  |                                                                                               
TCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTGCAAGAGGGGG  >  NC_000913/4500713‑4500906

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: