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)
       gAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAg                                                                                        <  2:1930142/100‑1 (MQ=255)
        aaaaTCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAgg                                                                                       <  1:888942/100‑1 (MQ=255)
        aaaaTCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAgg                                                                                       <  1:1949040/100‑1 (MQ=255)
         aaaTCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGt                                                                                      <  1:50728/100‑1 (MQ=255)
            tCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCt                                                                                   <  1:1117286/100‑1 (MQ=255)
            tCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCt                                                                                   <  1:1255829/100‑1 (MQ=255)
             caccacGAATACCAGAATCGCGGCACTTACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCt                                                                                   <  1:1079771/99‑1 (MQ=255)
             caccacGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTg                                                                                  >  1:2182102/1‑100 (MQ=255)
              accacGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGa                                                                                 >  2:1324096/1‑100 (MQ=255)
              accacGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGa                                                                                 >  1:1167069/1‑100 (MQ=255)
              accacGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGa                                                                                 >  2:386336/1‑100 (MQ=255)
               ccacGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGaa                                                                                >  1:914408/1‑100 (MQ=255)
                 acGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAAcc                                                                              >  1:2276326/1‑100 (MQ=255)
                        ccAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTTCTGTAGGTTCTGAACCGGTTCt                                                                        <  1:1850041/100‑1 (MQ=255)
                        ccAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTTCTGTAGGTTCTGAACCGGTTCt                                                                        >  1:1760890/1‑100 (MQ=255)
                        ccAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTa                                                                       >  1:1286109/1‑100 (MQ=255)
                        ccAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTa                                                                       >  1:2465088/1‑100 (MQ=255)
                            aaTCGCGGCACTGACCAGAAAATGATAATCAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCgg                                                                   <  2:1505701/100‑1 (MQ=255)
                            aaTCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCgg                                                                   <  1:970046/100‑1 (MQ=255)
                             aTCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCggg                                                                  <  2:1286101/100‑1 (MQ=255)
                                gcgGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGtt                                                               >  2:177805/1‑100 (MQ=255)
                                gcgGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGtt                                                               >  2:870628/1‑100 (MQ=255)
                                 cgGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTtg                                                              <  2:87513/100‑1 (MQ=255)
                                 cgGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACCTTGCTTTCTCCTGTAGTG‑CTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTtg                                                              >  1:2232019/1‑100 (MQ=255)
                                   gCACTGACCAGAAAATGATAAACAGGAAACACGGAGCGGTTCACTATGCTTTCTCCTGTAGTG‑TTTCTGTAGGTACTCAACCGGTTATAGCGGGGTtg                                                              >  2:203988/1‑98 (MQ=255)
                                     aCTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACGGGTTCTAGCGGGGTTGTGAg                                                          <  2:1569102/100‑1 (MQ=255)
                                     aCTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAg                                                          >  2:2278972/1‑100 (MQ=255)
                                     aCTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAg                                                          <  1:395587/100‑1 (MQ=255)
                                     aCTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACAGGTTCTAGCGGGGTTGTGAg                                                          <  2:1569118/100‑1 (MQ=255)
                                      cTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTTCTGTAGGTTCTGAGCCGGTTCTAGCGGGGTTGTGAg                                                          <  1:289076/100‑1 (MQ=255)
                                      cTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAg                                                          <  2:1686449/99‑1 (MQ=255)
                                       tGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCGGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGc                                                        <  2:225235/100‑1 (MQ=255)
                                           cAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGAt                                                    <  1:547333/100‑1 (MQ=255)
                                           cAGAAAAGGATAAACAGGAGACACGGGGCGGTTCTCTTTGCTTTCTCCTGAAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGCTTTGAGGCGGAt                                                    >  2:2343263/1‑100 (MQ=255)
                                            aGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATg                                                   >  1:2332758/1‑100 (MQ=255)
                                            aGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATg                                                   <  2:1673380/100‑1 (MQ=255)
                                            aGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATg                                                   <  1:716616/100‑1 (MQ=255)
                                            aGAAAATGATAAACAGGAAACACGGGGCGGTTCAATTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATg                                                   >  2:1935614/1‑100 (MQ=255)
                                             gAAAATGATGAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATgg                                                  <  1:907829/100‑1 (MQ=255)
                                             gAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATgg                                                  <  1:928198/100‑1 (MQ=255)
                                              aaaaTGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGc                                                 <  1:849683/100‑1 (MQ=255)
                                               aaaTGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGc                                                 <  2:1760877/100‑1 (MQ=255)
                                                 aTGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAAt                                              >  1:440418/1‑100 (MQ=255)
                                                    aTAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGc                                           <  2:915939/100‑1 (MQ=255)
                                                    aTAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGc                                           >  2:1003324/1‑100 (MQ=255)
                                                     tAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGcc                                          >  1:1484165/1‑100 (MQ=255)
                                                       atctagAAACACGGGGCGGTTCAGTTTGCTTTGTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCgg                                        <  1:803066/95‑1 (MQ=255)
                                                       aaCAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCgg                                        <  2:2543057/100‑1 (MQ=255)
                                                       aaCAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCgg                                        >  2:1635750/1‑100 (MQ=255)
                                                        aCAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCggg                                       <  1:1204827/100‑1 (MQ=255)
                                                         cAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGc                                      >  1:369328/1‑100 (MQ=255)
                                                         cAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGc                                      <  1:496472/100‑1 (MQ=255)
                                                            gAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGctct                                   <  2:1223470/100‑1 (MQ=255)
                                                             aaaCACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCtt                                  >  1:1871961/1‑100 (MQ=255)
                                                              aaCACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTc                                 <  2:589668/100‑1 (MQ=255)
                                                                cacGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCt                               >  2:2097199/1‑100 (MQ=255)
                                                                cacGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCt                               >  1:1198638/1‑100 (MQ=255)
                                                                cacGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCt                               >  2:1384965/1‑100 (MQ=255)
                                                                  cGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGc                             >  2:1451785/1‑100 (MQ=255)
                                                                   ggggCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTGTGTAGGTTCTGAACCGGTTCTAGCGGTGTTGTGAGGCTGATGGCAATCGCGGGGCTCTTCCTGCa                            >  2:142486/1‑100 (MQ=255)
                                                                    gggCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCAt                           >  2:1499527/1‑100 (MQ=255)
                                                                     ggCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATg                          <  2:1038257/100‑1 (MQ=255)
                                                                     ggCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATg                          <  1:671788/100‑1 (MQ=255)
                                                                       cGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCCTGtt                        >  1:1567705/1‑100 (MQ=255)
                                                                        ggTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTa                       >  2:783852/1‑100 (MQ=255)
                                                                        ggTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTa                       >  1:272272/1‑100 (MQ=255)
                                                                        ggTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTa                       <  2:1532995/100‑1 (MQ=255)
                                                                        ggTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTa                       <  2:1775446/100‑1 (MQ=255)
                                                                         gatctcTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTac                      <  1:2027299/95‑1 (MQ=255)
                                                                         gTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTTTTCCTGCATGTTac                      <  2:1310077/100‑1 (MQ=255)
                                                                         gTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTac                      <  1:663181/100‑1 (MQ=255)
                                                                           tCACTTTGCTTTCTCCTGT‑GTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACt                   <  1:1994258/100‑1 (MQ=255)
                                                                            cACTTTGCTTTCTCCTGTCGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTacac                    >  2:344556/1‑99 (MQ=255)
                                                                            cACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACt                   >  2:624742/1‑100 (MQ=255)
                                                                             aCTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTg                  >  1:1173911/1‑100 (MQ=255)
                                                                              cTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTacagatc                 >  2:2027285/1‑96 (MQ=255)
                                                                               tttGCTTTCTCCTGTAGTG‑TTTCTGTGGGTTCTGAACCGGTTCTATCGGGGTTGTGAGGCTGATGGCAATCGCCGGACTCTTCCTGCATGTTACACTGGt                <  2:1377237/100‑1 (MQ=255)
                                                                                ttGCTTTTTCCTGTAGTG‑TTTCTGTAGGTTCTGAAACGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTa               <  2:29061/100‑1 (MQ=255)
                                                                                 tGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTAc              <  1:581294/100‑1 (MQ=255)
                                                                                   cTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTg            >  1:2466709/1‑100 (MQ=255)
                                                                                    tttCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTGc           <  2:1762317/100‑1 (MQ=255)
                                                                                    tttCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTGc           >  1:2491002/1‑100 (MQ=255)
                                                                                         ccTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTGCAagag      >  2:783899/1‑100 (MQ=255)
                                                                                            gTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTGCAAGAgggg   <  1:199608/100‑1 (MQ=255)
                                                                                             tAGTGTTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTGCAAGAgggg   >  1:710840/1‑100 (MQ=255)
                                                                                             tAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTGCAAGAggggg  <  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: ‑