Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,596,874 Δ1 bp 18.3% intergenic (‑137/+31) 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,8730C.18.3% 436.1 / 98.4 153intergenic (‑136/+32)ypfJ/purCputative neutral zinc metallopeptidase/phosphoribosylaminoimidazole‑succinocarboxamide synthetase
Reads supporting (aligned to +/- strand):  ref base C (58/66);  new base . (11/17);  total (69/84)
Fisher's exact test for biased strand distribution p-value = 5.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.88e-01

Only 100 of 158 total aligned reads displayed.

TAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTC  >  NC_000913/2596775‑2596971
                                                                                                  |                                                                                                   
tatgtggtTAAGAGGAAGATTACCACGGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCGATCCATATATTCCGAGAGCGGTTGTCTTGCCGTGcc                                                                                                    <  1:801223/97‑1 (MQ=255)
tAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGcc                                                                                                    <  1:801208/100‑1 (MQ=255)
 aGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGcc                                                                                                    <  1:966338/99‑1 (MQ=255)
  ggtggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAg                                                                                                  >  1:4084461/1‑100 (MQ=255)
    tggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCt                                                                                                >  1:2747059/1‑100 (MQ=255)
     ggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTg                                                                                              >  1:1504568/1‑100 (MQ=255)
      gtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGc                                                                                             >  2:494723/1‑100 (MQ=255)
      gtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGc                                                                                             >  2:1333348/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:4060082/1‑100 (MQ=255)
          agagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACgg                                                                                         >  2:253686/1‑100 (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)
             ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCaa                                                                                      <  2:4237984/100‑1 (MQ=255)
              gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCAAg                                                                                     <  1:1158992/100‑1 (MQ=39)
               aaGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAga                                                                                    >  2:3251921/1‑100 (MQ=255)
               aaGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCAAga                                                                                    <  1:2132609/100‑1 (MQ=39)
               aaGATTACCACAGGAAAGCGGGATGAAATGGCGGACGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCa                                                                                       <  1:2256807/97‑1 (MQ=255)
                aGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgat                                                                                   >  2:1592512/1‑100 (MQ=255)
                 gATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCATGCAGCTGGCACGGCAAgaca                                                                                  <  1:3396773/100‑3 (MQ=38)
                  aTTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgatga                                                                                 <  2:1404729/100‑1 (MQ=255)
                    tACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGatgatg                                                                               <  1:2278360/100‑1 (MQ=255)
                    tACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGatgatg                                                                               <  1:570466/100‑1 (MQ=255)
                     aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGAtgatga                                                                              >  2:1063862/1‑100 (MQ=255)
                     aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGAtgatga                                                                              <  2:3922605/100‑1 (MQ=255)
                     aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCAAgacaaccg                                                                              <  2:2027662/100‑7 (MQ=37)
                     aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCAAgacaac                                                                                <  1:3893976/98‑5 (MQ=38)
                         caGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAg                                                                          >  2:2950401/1‑100 (MQ=255)
                           ggAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGaa                                                                        >  1:2645795/1‑100 (MQ=255)
                           ggAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGaa                                                                        >  2:515771/1‑100 (MQ=255)
                              aaGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGaaaaa                                                                     <  2:448947/100‑1 (MQ=255)
                               aGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCTAGCTGCACGGCAAGATGCTGAACAGAAAAAt                                                                    >  2:1792833/1‑100 (MQ=255)
                                gCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATc                                                                   >  2:642203/1‑100 (MQ=255)
                                    gATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTc                                                               >  2:1860713/1‑100 (MQ=255)
                                        aaaTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGc                                                           >  2:1764450/1‑100 (MQ=255)
                                          aTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTg                                                         <  2:1309758/100‑1 (MQ=255)
                                          aTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTg                                                         <  1:1082956/100‑1 (MQ=255)
                                            ggCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGtt                                                       <  2:1980990/100‑2 (MQ=255)
                                              cGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTaca                                                     >  2:1530203/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)
                                                  ggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTTCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCa                                                 <  1:3088058/100‑1 (MQ=255)
                                                  ggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCa                                                 <  2:2796951/100‑1 (MQ=255)
                                                      tGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcg                                             <  2:3371583/100‑1 (MQ=255)
                                                         gCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCGGgcggcg                                          <  2:2315511/100‑1 (MQ=255)
                                                            ttCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACGCAGGCggcggg                                        >  2:4135327/1‑99 (MQ=255)
                                                              cAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcgggc                                       >  2:4175895/1‑98 (MQ=255)
                                                               aaTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGAc                                    >  2:3802741/1‑100 (MQ=255)
                                                                 ggCATATGTTCCGAGAGGGGATGTCGTGCCGTGACAGCGGCACGGCAATATGATGACCAGAAAAATGAGTCCAGCTGTACACCCAGGCGGCGGGCGACgg                                  <  1:3211163/99‑1 (MQ=255)
                                                                  gatcttatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGc                                 <  1:840085/95‑1 (MQ=255)
                                                                   cATATATTCCGAGATCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCt                                <  1:679548/100‑1 (MQ=255)
                                                                    atatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCtt                               >  2:3358550/1‑100 (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:4023556/1‑100 (MQ=255)
                                                                         ttCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAg                          <  1:4135296/100‑1 (MQ=255)
                                                                          tCCGAGGGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAg                          <  2:469214/99‑1 (MQ=255)
                                                                            cGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTATACCCAGGCGGCGGGCGATGGCTTCATAggtt                       <  2:2971689/100‑3 (MQ=255)
                                                                            cGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCt                       >  1:3241114/1‑100 (MQ=255)
                                                                            cGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCt                       <  2:3582357/100‑1 (MQ=255)
                                                                              agagCGGTTGTCTTGCCGTTCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTc                     >  1:1106317/1‑100 (MQ=255)
                                                                              agagCGGTTGTCTTGCCGTGCCAGCTGCCCGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTc                     <  1:3861644/100‑1 (MQ=255)
                                                                              agagCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTc                     <  1:178235/100‑1 (MQ=255)
                                                                               gagCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCg                    <  2:62488/100‑1 (MQ=255)
                                                                                  tcatcatCTTGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATc                 >  2:531289/7‑100 (MQ=37)
                                                                                  cGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGAtn                 <  2:481844/100‑2 (MQ=255)
                                                                                  cGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATc                 <  2:2406355/100‑1 (MQ=255)
                                                                                  cGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATc                 >  1:3613167/1‑100 (MQ=255)
                                                                                   ggTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCa                <  2:691063/100‑1 (MQ=255)
                                                                                   ggTTGCCTTGCCGTGCCAGCTGCACGTCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCgcac                 <  1:1302912/99‑4 (MQ=255)
                                                                                     ttGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAgg              <  1:4253488/100‑1 (MQ=255)
                                                                                     ttGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAgg              <  2:3040771/100‑1 (MQ=255)
                                                                                       gTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGcc            <  1:924456/100‑1 (MQ=255)
                                                                                       atCTTGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGcc            >  2:3119400/2‑100 (MQ=39)
                                                                                        tCTTGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCa           >  2:2143205/1‑100 (MQ=39)
                                                                                         cTTGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTGCAGCTTTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAc          >  1:2598471/1‑100 (MQ=38)
                                                                                         cTTGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAc          >  1:734904/1‑100 (MQ=39)
                                                                                          ttGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTCCACCCAGGCGGCGCCACCCGGCTTCATACGCTCCGATCACGcc            >  1:2782650/1‑97 (MQ=255)
                                                                                          ttGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAcc         >  1:1883874/1‑100 (MQ=255)
                                                                                           tGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCg        >  1:2039459/1‑100 (MQ=255)
                                                                                           tGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCg        >  2:2099802/1‑100 (MQ=255)
                                                                                           tGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCg        >  2:2921678/1‑100 (MQ=255)
                                                                                           tGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCg        >  2:4111783/1‑100 (MQ=255)
                                                                                             ccGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACTGCTTCATAGGCTTCGATCAGGCCACCGAg      <  1:2466999/100‑1 (MQ=255)
                                                                                             ccGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAg      >  2:2263998/1‑100 (MQ=255)
                                                                                              cGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCCGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCCTACAAGGCTTCGATCATGCCACCGAgg     >  1:3590750/1‑100 (MQ=255)
                                                                                               gTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGc    <  1:1528005/100‑1 (MQ=255)
                                                                                                 gCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGctc  >  2:1148053/1‑100 (MQ=255)
                                                                                                 gCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGctc  <  1:2208668/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
TAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTC  >  NC_000913/2596775‑2596971

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: