Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,596,879:1 +G 16.1% intergenic (‑142/+26) ypfJ ← / ← purC putative neutral zinc metallopeptidase/phosphoribosylaminoimidazole‑succinocarboxamide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,596,8781.G16.1% 450.1 / 91.8 137intergenic (‑141/+27)ypfJ/purCputative neutral zinc metallopeptidase/phosphoribosylaminoimidazole‑succinocarboxamide synthetase
Reads supporting (aligned to +/- strand):  ref base . (52/63);  new base G (11/11);  total (63/74)
Fisher's exact test for biased strand distribution p-value = 8.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.95e-01

Only 100 of 142 total aligned reads displayed.

GGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGG  >  NC_000913/2596780‑2596977
                                                                                                   |                                                                                                   
ggtTCAGAGGAAGATTACCACAGGAAAGCGCGACGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTg                                                                                                    >  1:1917050/1‑100 (MQ=255)
 gtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGc                                                                                                   <  2:1637175/100‑1 (MQ=255)
  ttCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCa                                                                                                  <  1:238561/100‑1 (MQ=255)
  ttCAGAGGAAGATTACCACACGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAGTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCa                                                                                                  >  1:2041033/1‑100 (MQ=255)
   tCAGAGGCAGGTTACCACAGGAAAGCGGGATGAAATCGCGGAGGAATGTGTGTGCAATGCATATATTCCGGGAACGGTTGTCTTGCCGTGCCAGCTGCAc                                                                                                 <  1:441863/100‑1 (MQ=255)
    cAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACg                                                                                                >  2:2334680/1‑100 (MQ=255)
      gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGc                                                                                              >  2:1269666/1‑100 (MQ=255)
      gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGc                                                                                              >  1:949943/1‑100 (MQ=255)
      gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGc                                                                                              >  1:3908/1‑100 (MQ=255)
      gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGc                                                                                              <  1:2233936/100‑1 (MQ=255)
       agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCa                                                                                             <  2:929689/100‑1 (MQ=255)
       agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCa                                                                                             <  2:2030778/100‑1 (MQ=255)
       agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCa                                                                                             <  1:2327472/100‑1 (MQ=255)
        ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCaa                                                                                            <  2:2531787/100‑1 (MQ=255)
        ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCaa                                                                                            <  2:1088142/100‑1 (MQ=255)
        ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCaa                                                                                            <  2:2357695/100‑1 (MQ=255)
        ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCaa                                                                                            <  2:717183/100‑1 (MQ=255)
        ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGAGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCGCGGCaa                                                                                            >  1:558423/1‑100 (MQ=255)
         gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCATGGCAAg                                                                                           <  1:1670128/100‑1 (MQ=255)
          aaGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGCGCGGTTGTCTTGCCGTGCCAGCTGCAc                                                                                                 >  1:2722339/1‑93 (MQ=255)
           aGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCAAgac                                                                                         <  2:767666/100‑2 (MQ=39)
            gATTACCACAGGAAAGCGTGATGAAATGGCTGTGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTTTCTTGCAGCTGGCACGGCAAg                                                                                           <  1:497002/97‑1 (MQ=255)
             aTTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgatga                                                                                       <  1:1614946/100‑1 (MQ=255)
             aTTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgatga                                                                                       <  2:1163401/100‑1 (MQ=255)
             aTTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgatga                                                                                       <  2:1498547/100‑1 (MQ=255)
                aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGAtgatga                                                                                    <  1:182786/100‑1 (MQ=255)
                aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGAtgatga                                                                                    >  2:516565/1‑100 (MQ=255)
                aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGAtgatga                                                                                    >  2:1723403/1‑100 (MQ=255)
                 ccACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGaa                                                                                   >  2:1875401/1‑100 (MQ=255)
                  cacaGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAAc                                                                                  >  1:238373/1‑100 (MQ=255)
                   acaGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACa                                                                                 <  2:450234/100‑1 (MQ=255)
                   acaGGAAAGCGGGATGAAATGGCGGAGCAATGAGCGTTCAATGCATATATTCCGAGCGCGGCTGTCTTGCCGTGCCAGCTGCAGGGCAAGATGATGGAc                                                                                  >  1:134552/1‑99 (MQ=255)
                    caGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAg                                                                                >  1:829676/1‑100 (MQ=255)
                     aGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGa                                                                               >  1:1471946/1‑100 (MQ=255)
                          aGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAAt                                                                          >  2:79057/1‑100 (MQ=255)
                           gCGGTATGAGATGGGGGAGGAATGAGCGTTGAATTCAGATGTGGCGAGGGCGGTTGTCTTGCCGTGCCAGCTGTGCGGCAAGATGATGAACAGAAAAATc                                                                         <  1:2326839/100‑1 (MQ=255)
                           gCGGGATGAAATGGCGGAGGAATGAGCGTTCAATTCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGTACAGAAAAATc                                                                         <  1:1529476/100‑1 (MQ=255)
                            cGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCa                                                                        >  2:1053599/1‑100 (MQ=255)
                             cgaaTGAAATGGCGGATGAATGAGCGTTCAACCCATCTCTTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAACTc                                                                         <  1:1357579/95‑1 (MQ=255)
                                aTGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTcc                                                                    <  2:110001/100‑1 (MQ=255)
                                  gAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAg                                                                  >  1:2503617/1‑100 (MQ=255)
                                  gAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAg                                                                  <  2:631777/100‑1 (MQ=255)
                                    aaTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCt                                                                <  2:1471952/100‑1 (MQ=255)
                                       ggCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTa                                                             <  2:1310088/100‑1 (MQ=255)
                                       ggCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCACTCCAGCTGTa                                                             >  1:2416987/1‑100 (MQ=255)
                                        gCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTac                                                            <  2:610398/100‑1 (MQ=255)
                                        gCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTac                                                            <  2:613965/100‑1 (MQ=255)
                                         cGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTaca                                                           >  2:1044559/1‑100 (MQ=255)
                                         cGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTaca                                                           <  1:2268582/100‑1 (MQ=255)
                                          ggaggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCCGCTGTacac                                                          >  1:2601302/1‑100 (MQ=255)
                                             ggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCa                                                       >  2:256465/1‑100 (MQ=255)
                                               aaTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAATAGAAAAATCAGGCCAGCTGTACACCCAgg                                                     <  1:1195244/100‑1 (MQ=255)
                                                aTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggc                                                    >  2:635074/1‑100 (MQ=255)
                                                 tGAGCGTTCAATGCCTATATTCCGAGAGCGGTTGTTTTGCCGTCCCAGCTGCACGGCAAGAGGATGAACAGAAAAATCAGTCCCGCTGTACACCCAgtc                                                    >  1:1703126/1‑97 (MQ=255)
                                                   aGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcggc                                                 >  2:1633024/1‑100 (MQ=255)
                                                   aGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggc                                                    >  2:875820/1‑97 (MQ=255)
                                                    gCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCCCGGCAATATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcg                                                   >  1:628309/1‑97 (MQ=255)
                                                       ttCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcgggc                                             <  2:2673905/100‑1 (MQ=255)
                                                         cAATACATATATTCCGAGAGCGGTTGACTTGCCGTGCCAGCTGCACGGCAAGAAAAAAAAAAAAAAAATCAGTCCAGCTGTACACCCAGGCggcgggcg                                            <  1:2868482/99‑1 (MQ=255)
                                                            tGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACgg                                        <  2:167127/100‑1 (MQ=255)
                                                            tGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACgg                                        >  2:122278/1‑100 (MQ=255)
                                                               atatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCtt                                     >  2:2687622/1‑100 (MQ=255)
                                                               atatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCtt                                     >  1:2059699/1‑100 (MQ=255)
                                                                tatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTc                                    <  2:1426113/100‑1 (MQ=255)
                                                                tatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTc                                    >  2:1417897/1‑100 (MQ=255)
                                                                tatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTc                                    <  1:927508/100‑1 (MQ=255)
                                                                 atatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCa                                   <  2:2410477/100‑1 (MQ=255)
                                                                   atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                                 <  1:2851937/100‑1 (MQ=255)
                                                                    ttCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAg                                >  1:874699/1‑100 (MQ=255)
                                                                           agCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGa                         >  2:1826605/1‑100 (MQ=255)
                                                                            gCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGAt                        >  2:1655084/1‑100 (MQ=255)
                                                                             cGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATc                       >  2:2832086/1‑100 (MQ=255)
                                                                             cGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGAGGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATc                       <  1:1052972/100‑1 (MQ=255)
                                                                               gTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAg                     <  1:1713357/100‑1 (MQ=255)
                                                                               gTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAg                     <  2:1932583/100‑1 (MQ=255)
                                                                               gTTGGCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAg                     <  1:1940696/100‑1 (MQ=255)
                                                                                ttGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAgg                    <  2:2794309/100‑1 (MQ=255)
                                                                                ttGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAgg                    <  1:919118/100‑1 (MQ=255)
                                                                                ttGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAgg                    <  2:1685057/100‑1 (MQ=255)
                                                                                 catCTTGCCGTGCAGCTGGCACGGCAAGATGATGAACAGTAAAATCAGTCCAGCTGTACACTCTGGCGGCGGTCGACGGCTTCATATGCTTCGagctgg                    >  1:2621103/3‑94 (MQ=25)
                                                                                 catCTTGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGc                   >  1:2398215/3‑100 (MQ=38)
                                                                                  gTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGcc                  <  2:1632600/100‑1 (MQ=255)
                                                                                  gTCTTGCCGTGCCAGCTGCACGGCAAGAGGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGc                   <  1:1895519/99‑1 (MQ=255)
                                                                                   tCTTGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCa                 >  2:2347423/1‑100 (MQ=39)
                                                                                     ttGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAcc               <  2:1121068/100‑1 (MQ=255)
                                                                                     ttGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAcc               >  2:1002542/1‑100 (MQ=255)
                                                                                      tGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCg              >  1:2094448/1‑100 (MQ=255)
                                                                                      tGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCg              >  1:396817/1‑100 (MQ=255)
                                                                                       gCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGa             >  2:220136/1‑100 (MQ=255)
                                                                                       gCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGa             <  2:829679/100‑1 (MQ=255)
                                                                                       gCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGa             >  2:975622/1‑100 (MQ=255)
                                                                                          gTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGc          <  2:2090040/100‑1 (MQ=255)
                                                                                          gTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGc          <  2:460814/100‑1 (MQ=255)
                                                                                          gTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGc          <  1:516562/100‑1 (MQ=255)
                                                                                            gCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGctc        >  2:1693013/1‑100 (MQ=255)
                                                                                             ccAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGctct       <  2:238374/100‑1 (MQ=255)
                                                                                             ccAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGctct       <  2:2519760/100‑1 (MQ=255)
                                                                                                 cTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCACACTAGTCTTCGAACCGGCCACCGAGGCTCTGGc    >  1:119897/1‑99 (MQ=255)
                                                                                                  tGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATTAGGCCACCGCGGATCTcgcg   >  1:1896198/1‑99 (MQ=255)
                                                                                                  tGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCg   >  1:1689152/1‑99 (MQ=255)
                                                                                                   |                                                                                                   
GGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGG  >  NC_000913/2596780‑2596977

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: