Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,596,874 Δ1 bp 15.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.15.3% 398.1 / 74.3 137intergenic (‑136/+32)ypfJ/purCputative neutral zinc metallopeptidase/phosphoribosylaminoimidazole‑succinocarboxamide synthetase
Reads supporting (aligned to +/- strand):  ref base C (62/52);  new base . (11/10);  total (74/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 141 total aligned reads displayed.

GTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCT  >  NC_000913/2596774‑2596970
                                                                                                   |                                                                                                  
gTAGGTGGTTCAGAGGAAGATTACCTCAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGc                                                                                                    >  2:2301654/1‑100 (MQ=255)
gTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGc                                                                                                    >  1:2237034/1‑100 (MQ=255)
 tAGGTGGATCAAAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGcc                                                                                                   >  2:2184885/1‑100 (MQ=255)
  aGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCa                                                                                                  <  2:639979/100‑1 (MQ=255)
    gtggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGc                                                                                                >  1:2883487/1‑100 (MQ=255)
      ggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGCGAGCGGTTGTCTTGCCGTGCCAGCTg                                                                                             >  2:483797/1‑100 (MQ=255)
       gtggagagGAAGATTACCAGAGGAAAGCGGGATGAAATGGGGGGGGAATGCGCGGTCAATGAATATATTCAGAGAGCGGTTGGCTTGCCGTGcga                                                                                                  <  2:230498/91‑3 (MQ=255)
       gtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGc                                                                                            >  1:927986/1‑100 (MQ=255)
          cAGAGGAAGATTACCACCGGAAAGCGGGATGACATGGCGGAGGAATGCGCGTTCACTGCATATATTCCGAGCGCGGTTGTCTTGCCGTGACAGCTGAACg                                                                                         >  2:850198/1‑100 (MQ=255)
          cAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCAAGAGCGGTTGTCTTGCCGTGCCAGCTGCACg                                                                                         <  2:3223852/100‑1 (MQ=255)
           agagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACgg                                                                                        <  1:4016413/100‑1 (MQ=255)
           agagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACg                                                                                         <  1:1572029/99‑1 (MQ=255)
            gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGc                                                                                       >  2:2359378/1‑100 (MQ=255)
             agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCa                                                                                      <  1:3432255/100‑1 (MQ=255)
              ttccgatctccACAGGAAAGCGGGATGAAATGGCCGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGATGGCACGGCaa                                                                                     <  2:2191190/91‑1 (MQ=255)
              ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCaa                                                                                     <  1:2135970/100‑1 (MQ=255)
              ggAAGATTAACACAGGAAAGCGGGATGAAATGGCGCAGCAATGAGCGTTCAATGCATATATTCCGAGAGGGGTTGTCTTGCCGTGCAGCTGGCACGGCaa                                                                                     <  2:1856338/100‑1 (MQ=255)
               gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCAAg                                                                                    <  2:2764994/100‑1 (MQ=39)
                aaGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCAAga                                                                                   <  1:1642133/100‑1 (MQ=39)
                 aGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGACCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgat                                                                                  >  2:1042850/1‑100 (MQ=255)
                  gATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgatg                                                                                 <  1:799984/100‑1 (MQ=255)
                  gATTAACACAGAAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTAAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCAAgaca                                                                                 <  2:2483844/100‑3 (MQ=37)
                     tACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGatgatg                                                                              <  1:397193/100‑1 (MQ=255)
                     tACCACACGAAAGCGGGATGAAATGGCTGATGAATGCGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCTTGCCAGCTGCACGGCAAGatgatg                                                                              <  1:3970458/100‑1 (MQ=255)
                      aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGAtgatga                                                                             <  1:3522248/100‑1 (MQ=255)
                      aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGAtgatga                                                                             <  2:2453928/100‑1 (MQ=255)
                       ccACAGGGAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCaaagatcggaa                                                                            >  1:2191204/1‑91 (MQ=255)
                       ccACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGaa                                                                            >  1:23764/1‑100 (MQ=255)
                          caGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAg                                                                         <  2:2074587/100‑1 (MQ=255)
                           aGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGa                                                                        >  1:173502/1‑100 (MQ=255)
                            ggAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGaa                                                                       >  2:3549044/1‑100 (MQ=255)
                             gAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGaaa                                                                      <  1:1634593/100‑1 (MQ=255)
                             gAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGaaa                                                                      <  1:603128/100‑1 (MQ=255)
                             gAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGaaa                                                                      <  1:3716500/100‑1 (MQ=255)
                              aaaGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGaaaa                                                                     <  1:3647855/100‑1 (MQ=255)
                              aaaGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGaaaa                                                                     <  1:1175299/100‑1 (MQ=255)
                               aaGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGaaaaa                                                                    >  2:1091602/1‑100 (MQ=255)
                               aaGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGACGTGCCAGCTGCACTGCAAGATGATGAACAGaaaaa                                                                    <  2:1207219/100‑1 (MQ=255)
                                 gCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGaaaaa                                                                    >  1:722141/1‑98 (MQ=255)
                                 gCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATc                                                                  <  2:1640352/100‑1 (MQ=255)
                                 gCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATc                                                                  <  1:3345616/100‑1 (MQ=255)
                                    ggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGt                                                               >  1:1915523/1‑100 (MQ=255)
                                     gATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTc                                                              <  1:2543311/100‑1 (MQ=255)
                                     gATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTATTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTc                                                              <  2:3354076/100‑1 (MQ=255)
                                        gAAATGGCGGAGGAATGAGCGTTCGATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAg                                                           >  1:3151532/1‑100 (MQ=255)
                                         aaaTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGc                                                          >  2:322435/1‑100 (MQ=255)
                                          aaTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCt                                                         >  2:3039517/1‑100 (MQ=255)
                                           aTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTg                                                        <  2:2759616/100‑1 (MQ=255)
                                              gCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTac                                                     >  2:783244/1‑100 (MQ=255)
                                               cGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTaca                                                    >  2:1451847/1‑100 (MQ=255)
                                               cGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTaca                                                    <  1:2336294/100‑1 (MQ=255)
                                                  aggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACAccc                                                 <  1:2439425/100‑1 (MQ=255)
                                                  aggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACAccc                                                 <  1:3728310/100‑1 (MQ=255)
                                                  aggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACAccc                                                 <  2:2515431/100‑1 (MQ=255)
                                                   ggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCa                                                <  1:893606/100‑1 (MQ=255)
                                                     aaTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCa                                                >  1:2402636/1‑98 (MQ=255)
                                                     aaTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAgg                                              <  1:1864123/100‑1 (MQ=255)
                                                     aaTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAgg                                              >  2:2115005/1‑100 (MQ=255)
                                                        gAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcgg                                           >  1:2295096/1‑100 (MQ=255)
                                                           cGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGCTGAACAGAAAAATCACCCCAGCTGTACACCCa                                                >  2:3973062/1‑92 (MQ=255)
                                                           cGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGcggcgg                                        <  1:2116624/100‑1 (MQ=255)
                                                            gTTCAATGCATATATTCCTAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGCCCAGCTGTACACCCAGGcggcgg                                        >  2:1289732/1‑99 (MQ=255)
                                                            gTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcggg                                       <  2:1951/100‑1 (MQ=255)
                                                             ttCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcgggc                                      <  1:3241792/100‑1 (MQ=255)
                                                              tCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcgggcg                                     <  1:323521/100‑1 (MQ=255)
                                                                 aTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACg                                  <  1:322438/100‑1 (MQ=255)
                                                                 aTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACg                                  >  2:1930094/1‑100 (MQ=255)
                                                                  tGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACgg                                 >  1:2644148/1‑100 (MQ=255)
                                                                   gCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGc                                <  2:871933/100‑1 (MQ=255)
                                                                    cATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACgg                                 >  2:3263590/1‑98 (MQ=255)
                                                                    cATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCt                               <  1:1639095/100‑1 (MQ=255)
                                                                    cATATATTCCGAGAGCGGTTGTCTTACCGTGCCAGCTGCACAGCAAGATGATGAACAGAAAAATCACTCCAGCTATACACCCAAGCCGCCAGCGa                                    >  2:491740/1‑95 (MQ=255)
                                                                       atatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCa                            <  1:1894068/100‑1 (MQ=255)
                                                                        tatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCAt                           >  2:866690/1‑100 (MQ=255)
                                                                         atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                          >  2:3539835/1‑100 (MQ=255)
                                                                          ttCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAg                         <  2:2237613/100‑1 (MQ=255)
                                                                           tCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAgg                        >  1:317610/1‑100 (MQ=255)
                                                                            ccGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGc                       >  1:1210958/1‑100 (MQ=255)
                                                                               agagCGGTTGTCATTCCTTGACAGCTGCACTTCAAGTTGATGAACAGAAAAATCAGTCCAGCTGTACTCCCAGGCTGCGGGCGACGGCTTCATAg                         <  2:3931556/95‑1 (MQ=255)
                                                                                gagCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCg                   >  2:1555599/1‑100 (MQ=255)
                                                                                      ttGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAgg             <  1:1037227/100‑1 (MQ=255)
                                                                                      tcatCTTGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAgg             >  1:738614/4‑100 (MQ=38)
                                                                                         tCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCa          >  1:332206/1‑100 (MQ=255)
                                                                                         tCTTGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCa          >  1:564599/1‑100 (MQ=39)
                                                                                          cTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAc         >  1:2694277/1‑100 (MQ=255)
                                                                                          cTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAc         >  2:1703938/1‑100 (MQ=255)
                                                                                          cTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCCACGGCTTCATAGGCTTCGATCAGGCCAc         >  2:2303729/1‑100 (MQ=255)
                                                                                          cTTGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAc         >  2:2536396/1‑100 (MQ=39)
                                                                                          cTTGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAc         >  1:3604180/1‑100 (MQ=39)
                                                                                          cTTGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAc         >  2:2470503/1‑100 (MQ=39)
                                                                                           ttGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAcc        >  1:2325123/1‑100 (MQ=255)
                                                                                           ttGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGAGGGCGACGGCTTCATAGGCTTCGATCAGGCCAcc        >  1:292647/1‑100 (MQ=255)
                                                                                            tGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTGCGATCAGGCCACCg       >  1:3395662/1‑100 (MQ=255)
                                                                                              ccGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTTATAGGCGTCGATCAGGCCACCGAg     <  2:3466510/100‑1 (MQ=255)
                                                                                              ccGTGCCAGCTGCACGCCAAGATGATGAACAGAAAAGTCAGTCCAGCTGTACACCCAGGCGGCGGGAGATGGCTTTATAGGCTTCGATGCGGCCACCGAg     <  2:2400798/100‑1 (MQ=255)
                                                                                               cGTGCCAGCTGCACGGCAAGATGATGAACAGAACAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAgg    <  1:2202416/100‑1 (MQ=255)
                                                                                               cGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAgg    <  1:1930104/100‑1 (MQ=255)
                                                                                               cGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAgg    >  1:2435892/1‑100 (MQ=255)
                                                                                                gTGCCAGATGCACGGCAAGACGATGAGCAGAAAAATCAGTCCAGCTGTACCCCCAGGCGGCGGGCGACGGCTTCATGGGCTTCGATCAGGCCACCGAgg    <  2:1858582/99‑1 (MQ=255)
                                                                                                 tGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGct  >  1:1568577/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
GTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCT  >  NC_000913/2596774‑2596970

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: