Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,596,879:1 +G 19.9% 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.G19.9% 421.0 / 124.6 146intergenic (‑141/+27)ypfJ/purCputative neutral zinc metallopeptidase/phosphoribosylaminoimidazole‑succinocarboxamide synthetase
Reads supporting (aligned to +/- strand):  ref base . (51/66);  new base G (12/17);  total (63/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.41e-01

Only 100 of 152 total aligned reads displayed.

TGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGG  >  NC_000913/2596779‑2596977
                                                                                                    |                                                                                                   
tggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCt                                                                                                      >  1:2747059/1‑100 (MQ=255)
 ggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTg                                                                                                    >  1:1504568/1‑100 (MQ=255)
  gtTCAGCGGAAGATTACCACAGGAAACCGGGATGAAATGTCGGAGGAATGACCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTACCGTGCCAGCTg                                                                                                    >  2:2408916/1‑99 (MQ=255)
  gtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGc                                                                                                   >  2:1333348/1‑100 (MQ=255)
  gtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGc                                                                                                   >  2:494723/1‑100 (MQ=255)
    tCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCAc                                                                                                 >  1:211375/1‑100 (MQ=255)
     cAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACg                                                                                                <  2:3273092/100‑1 (MQ=255)
     cAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATAGATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACg                                                                                                <  1:2356099/100‑1 (MQ=255)
      agagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACgg                                                                                               <  2:675487/100‑1 (MQ=255)
      agagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACgg                                                                                               >  2:253686/1‑100 (MQ=255)
      agagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACgg                                                                                               <  1:4060052/100‑1 (MQ=255)
      agagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACg                                                                                                >  1:173569/1‑99 (MQ=255)
      agagGAAGATTACCACAGGAAAGCGGGATGAAATGGCAGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACg                                                                                                <  1:53460/99‑1 (MQ=255)
       gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGc                                                                                              >  2:3889558/1‑100 (MQ=255)
       gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGc                                                                                              <  1:3527836/100‑1 (MQ=255)
        agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCa                                                                                             <  2:2059640/100‑1 (MQ=255)
         ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTTAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCa                                                                                             >  1:4368678/1‑99 (MQ=255)
         ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCaa                                                                                            <  1:2856594/100‑1 (MQ=255)
         ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCaa                                                                                            <  2:2872847/100‑1 (MQ=255)
          gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCAAg                                                                                           <  2:2954893/100‑1 (MQ=39)
          gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCAAg                                                                                           <  1:1158992/100‑1 (MQ=39)
           aacaTTACCCCAACAAACCGGTATGAAATGGCGGAGGAATCAGCGTTCAACGCATATATTCCGATAGCGTTTGTCTTGCCTTGCAGCTGGCACGGCa                                                                                             <  1:3563226/94‑1 (MQ=255)
           aaGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAga                                                                                          >  2:3251921/1‑100 (MQ=255)
           aaGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAga                                                                                          <  2:2273214/100‑1 (MQ=255)
           aaGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCAAga                                                                                          <  2:2058543/100‑1 (MQ=39)
           aaGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCAAga                                                                                          <  1:2132609/100‑1 (MQ=39)
           aaGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCAAga                                                                                          <  1:3652805/100‑1 (MQ=39)
            aGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgat                                                                                         >  2:1592512/1‑100 (MQ=255)
              aTTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgatga                                                                                       <  2:1404729/100‑1 (MQ=255)
                tACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGatgatg                                                                                     <  1:2278360/100‑1 (MQ=255)
                 aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGAtgatga                                                                                    <  2:3922605/100‑1 (MQ=255)
                 aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCAAgacaaccg                                                                                    <  1:2759877/100‑7 (MQ=37)
                 aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCAAgacaac                                                                                      <  1:3893976/98‑5 (MQ=38)
                     caGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAg                                                                                >  2:2950401/1‑100 (MQ=255)
                       ggAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGaa                                                                              >  2:2893845/1‑100 (MQ=255)
                          aaGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGaaaaa                                                                           <  2:448966/100‑1 (MQ=255)
                          aaGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGaaaaa                                                                           <  2:4214741/100‑1 (MQ=255)
                           aGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCTAGCTGCACGGCAAGATGCTGAACAGAAAAAt                                                                          >  2:1792833/1‑100 (MQ=255)
                            gCGGGATGAAATGGCGGAGGAATGGGCGTTCAATTCATATATCCCGAGCGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGACGAACAGAAAAATc                                                                         <  1:26380/100‑1 (MQ=255)
                            gCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATc                                                                         <  2:1184713/100‑1 (MQ=255)
                            gCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATc                                                                         >  2:642203/1‑100 (MQ=255)
                                    aaaTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGc                                                                 >  2:1764450/1‑100 (MQ=255)
                                      aTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTg                                                               <  1:1082956/100‑1 (MQ=255)
                                        ggCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGtt                                                             <  2:1980990/100‑2 (MQ=255)
                                        ggCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTa                                                             <  2:1495367/100‑1 (MQ=255)
                                         gCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTac                                                            >  2:2019309/1‑100 (MQ=255)
                                          cGGAGGAATGAGCGTTCAATGCAAATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTaca                                                           >  1:18069/1‑100 (MQ=255)
                                           ggaggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTacac                                                          <  2:4404960/100‑1 (MQ=255)
                                            gaggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACAcc                                                         <  2:2748974/100‑1 (MQ=255)
                                             aggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACAccc                                                        <  2:2735499/100‑1 (MQ=255)
                                              ggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCa                                                       <  2:2796951/100‑1 (MQ=255)
                                                aaTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAgg                                                     <  2:2122328/100‑1 (MQ=255)
                                                  tGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcg                                                   <  2:3371583/100‑1 (MQ=255)
                                                    aGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcggc                                                 <  1:2228579/100‑1 (MQ=255)
                                                     gCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCGGgcggcg                                                <  2:2315511/100‑1 (MQ=255)
                                                        ttCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcgggc                                             <  2:3962444/100‑1 (MQ=255)
                                                          cAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcgggc                                             >  2:4175895/1‑98 (MQ=255)
                                                           aaTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGAc                                          >  1:3057666/1‑100 (MQ=255)
                                                           aaTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGAc                                          >  2:3802741/1‑100 (MQ=255)
                                                            aTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACg                                         >  1:3350926/1‑100 (MQ=255)
                                                             ggCATATGTTCCGAGAGGGGATGTCGTGCCGTGACAGCGGCACGGCAATATGATGACCAGAAAAATGAGTCCAGCTGTACACCCAGGCGGCGGGCGACgg                                        <  1:3211163/99‑1 (MQ=255)
                                                              gatcttatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGc                                       <  1:840085/95‑1 (MQ=255)
                                                              gCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGc                                       <  2:492286/100‑1 (MQ=255)
                                                               cATATATTCCGAGATCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCt                                      <  1:679548/100‑1 (MQ=255)
                                                                atatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCtt                                     >  2:3358550/1‑100 (MQ=255)
                                                                atatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCtn                                     <  1:941152/100‑2 (MQ=255)
                                                                 tatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTc                                    <  1:974249/100‑1 (MQ=255)
                                                                 tatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTc                                    <  2:3908581/100‑1 (MQ=255)
                                                                  atatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAGGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCa                                   >  1:1614059/1‑100 (MQ=255)
                                                                   tatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGcagatc                                  >  2:840088/1‑95 (MQ=255)
                                                                    atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                                 >  2:2714938/1‑100 (MQ=255)
                                                                     ttCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAg                                <  1:4135296/100‑1 (MQ=255)
                                                                      tCCGAGGGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAg                                <  2:469214/99‑1 (MQ=255)
                                                                      tCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAgg                               >  1:665193/1‑100 (MQ=255)
                                                                        cGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCt                             >  1:3241114/1‑100 (MQ=255)
                                                                        cGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCt                             <  2:3582357/100‑1 (MQ=255)
                                                                          agagCGGTTGTCTTGCCGTTCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTc                           >  1:1106317/1‑100 (MQ=255)
                                                                           gagCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCg                          <  2:62488/100‑1 (MQ=255)
                                                                            agCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGa                         <  2:1387086/100‑1 (MQ=255)
                                                                             gCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGAt                        <  1:3906854/100‑1 (MQ=255)
                                                                              tcatcatCTTGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATc                       >  2:531289/7‑100 (MQ=37)
                                                                              cGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATc                       >  1:3613167/1‑100 (MQ=255)
                                                                               ggTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCa                      <  2:691063/100‑1 (MQ=255)
                                                                               ggTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCa                      <  2:525540/100‑1 (MQ=255)
                                                                               ggTTGCCTTGCCGTGCCAGCTGCACGTCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCgcac                       <  1:1302912/99‑4 (MQ=255)
                                                                                 ttGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAgg                    <  2:3040771/100‑1 (MQ=255)
                                                                                   atCTTGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGcc                  >  2:3119400/2‑100 (MQ=39)
                                                                                    tCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCa                 >  2:4384593/1‑100 (MQ=255)
                                                                                    tCTTGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCa                 >  2:2143205/1‑100 (MQ=39)
                                                                                     cTTGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTGCAGCTTTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAc                >  1:2598471/1‑100 (MQ=38)
                                                                                      ttGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAcc               >  1:742127/1‑100 (MQ=255)
                                                                                      ttGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAcc               >  1:1883874/1‑100 (MQ=255)
                                                                                       tGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCg              >  2:4133253/1‑100 (MQ=255)
                                                                                       tGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCg              >  2:4111783/1‑100 (MQ=255)
                                                                                         ccGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACTGCTTCATAGGCTTCGATCAGGCCACCGAg            <  1:2466999/100‑1 (MQ=255)
                                                                                           gTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGc          <  1:1528005/100‑1 (MQ=255)
                                                                                             gCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGctc        <  1:2208668/100‑1 (MQ=255)
                                                                                             gCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGctc        >  2:1148053/1‑100 (MQ=255)
                                                                                                  cTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCg   <  2:3071737/100‑1 (MQ=255)
                                                                                                   tGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCgg  <  1:2262334/100‑1 (MQ=255)
                                                                                                    |                                                                                                   
TGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGG  >  NC_000913/2596779‑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: