Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,313,576 A→T 6.3% D238V (GAC→GTC)  hemB → 5‑aminolevulinate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,313,5760AT 6.3% 324.0 / 8.7 106D238V (GAC→GTC) hemB5‑aminolevulinate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (49/50);  new base T (2/5);  total (51/55)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 114 total aligned reads displayed.

ACCGCGATCATGGCCTACTCGACCAAGTTCGCCTCGTCGTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGCGCCGCTTCGCGCCTGCCACTGGCGGCGTACCAGGTGAGTGGCGAGTATGCGATGATCAAGTTCGGTGGCCTGGCCG  >  NC_002947/3313428‑3313725
                                                                                                                                                     |                                                                                                                                                     
acCGCGATCATGGCCTACTCGACCAAGTTCGCCTCGTCGTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTcgac                                                                                                                                                      <  1:438800/150‑1 (MQ=255)
  cgcgATCATGGCCTACTCGACCAAGTTCGCCTCGTCGTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGGTGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTcgacga                                                                                                                                                    >  2:940718/1‑150 (MQ=255)
  cgcgATCATGGCCTACTCGACCAAGTTCGCCTCGTCGTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTcgacga                                                                                                                                                    >  2:1586741/1‑150 (MQ=255)
     gATCATGGCCTACTCGACCAAGTTCGCCTCGTCGTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGCGTCGCTGCTCGACGAGCa                                                                                                                                                 >  1:292057/1‑150 (MQ=255)
     gATCATGGCCTACTCGACCAAGTTCGCCTCGTCGTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCa                                                                                                                                                 <  2:1822651/150‑1 (MQ=255)
     gATCATGGCCTACTCGACCAAGTTCGCCTCGTCGTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCa                                                                                                                                                 >  2:1919882/1‑150 (MQ=255)
         aTGGCCTACTCGACCAAGTTCGCCTCGTCGTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGTCGAGCAgcga                                                                                                                                             <  2:595349/150‑4 (MQ=255)
         aTGGCCTACTCGACCAAGTTCGCCTCGTCGTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCCCTGCTCGTCGAGCAgcga                                                                                                                                             <  1:2187998/150‑4 (MQ=255)
          tGGCCTACTCGACCAAGTTCGCCTCGTCGTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTcgacga                                                                                                                                                    >  1:1050539/1‑142 (MQ=255)
          tGGCCTACTCGACCAAGTTCGCCTCGTCGTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTcgacga                                                                                                                                                    <  2:1050539/142‑1 (MQ=255)
             ccTACTGGACCAAGTTCGCCTCGTCGTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAgggg                                                                                                                                         <  2:2074336/150‑1 (MQ=255)
             ccTACTCGACAAACTGCGCCTCGTCGTTGTATGGCCCGTTCCGTGAGACAGGCCAGCAGGCGGTGGTCGGGGCTGGTAAAAGCTACCAGATGAACCCGATGAACCGCCGCGATGTGCTGCGGGAGTCGCTGCTCGACGAGCAGGAAgggg                                                                                                                                         <  1:1151526/150‑1 (MQ=255)
                    gACCAAGTTCGCCTCGTCGTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGGGCGACAGACACGGCGCGACTCGCATGGGGGGGCGGACGt                                                                                                                                  >  2:737205/1‑150 (MQ=255)
                     aCCAAGTTCGCCTCGTCGTTGTATGGCCCGTTCCGTAAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTg                                                                                                                                 <  2:2037941/150‑1 (MQ=255)
                      ccAAGTTCGCCTCGTCGTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGc                                                                                                                                >  2:1165951/1‑150 (MQ=255)
                      ccAAGTTCGCCTCGTCGTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGc                                                                                                                                >  2:1197270/1‑150 (MQ=255)
                      ccAAGTTCGCCTCGTCGTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGc                                                                                                                                <  2:1659728/150‑1 (MQ=255)
                      ccAAGTTCGCCTCGTCGTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGc                                                                                                                                <  1:1197270/150‑1 (MQ=255)
                      ccAAGTTCGCCTCGTCGTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGc                                                                                                                                <  1:1165951/150‑1 (MQ=255)
                       cAAGTTCGCCTCGTCGTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCt                                                                                                                               >  2:1160852/1‑150 (MQ=255)
                       cAAGTTCGCCTCGTCGTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCt                                                                                                                               <  2:356315/150‑1 (MQ=255)
                        aaGTTCGCCTCGTCGAGGGAATTCCAGATCCGGTAGGCCGGCTGGACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTg                                                                                                                              >  2:1954859/1‑150 (MQ=255)
                         aGTTCGCCTCGTCGTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGCAATGGCGGGAAGCTGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGa                                                                                                                             <  2:2393277/150‑1 (MQ=255)
                          gTTCGCCTCGTCGTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGAt                                                                                                                            >  1:2367904/1‑150 (MQ=255)
                           ttCGCCTCGTCGTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATg                                                                                                                           >  1:1083410/1‑150 (MQ=255)
                           ttCGCCTCGTCGTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATg                                                                                                                           <  2:1083410/150‑1 (MQ=255)
                             cGCCTCGTCGTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGt                                                                                                                         <  2:2167550/150‑1 (MQ=255)
                                cTCGTCGTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAgggg                                                                                                                                         <  2:57508/131‑1 (MQ=255)
                                cTCGTCGTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAgggg                                                                                                                                         >  1:57508/1‑131 (MQ=255)
                                cTCGTCGTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCaa                                                                                                                      <  2:2183163/150‑1 (MQ=255)
                                 tcgtcgTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAg                                                                                                                     >  1:1288678/1‑150 (MQ=255)
                                            tGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGc                                                                                                          <  1:1919882/150‑1 (MQ=255)
                                            tGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGc                                                                                                          >  2:1752557/1‑150 (MQ=255)
                                               cccGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATg                                                                                                                           <  1:123622/130‑1 (MQ=255)
                                               cccGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATg                                                                                                                           >  2:123622/1‑130 (MQ=255)
                                                ccGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAgccggcctgtt                                                                                                           >  1:1962368/1‑141 (MQ=255)
                                                ccGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCTGtt                                                                                                           <  2:1962368/145‑2 (MQ=255)
                                                 cGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGGCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCGGATGGTAAAGACGGCCGGCGCGTAcc                                                                                                     >  1:1005422/1‑150 (MQ=255)
                                                  gTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCAGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCt                                                                                                    <  2:1288678/150‑1 (MQ=255)
                                                     ccGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGa                                                                                                 <  2:888178/150‑1 (MQ=255)
                                                      cGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGAt                                                                                                >  2:2254214/1‑150 (MQ=255)
                                                      cGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGAt                                                                                                <  1:2178296/150‑1 (MQ=255)
                                                       gTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGTGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGt                                                                                                            >  1:1063779/1‑137 (MQ=255)
                                                       gTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGTGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGt                                                                                                            <  2:1063779/137‑1 (MQ=255)
                                                             ccGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAgg                                                                                                                                           <  2:2387525/100‑1 (MQ=255)
                                                             ccGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAgg                                                                                                                                           >  1:2387525/1‑100 (MQ=255)
                                                             ccGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAgccg                                                                                                                  >  2:2290830/1‑125 (MQ=255)
                                                             ccGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAgccg                                                                                                                  <  1:2290830/125‑1 (MQ=255)
                                                                gcggcACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCg                                                                                      <  1:251285/150‑1 (MQ=255)
                                                                  cggcACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGa                                                                                     <  2:1622386/150‑1 (MQ=255)
                                                                  cagcACCGCACTGAAGAATGACCGGAAAAGCTACCTGATGATCCCTACGAACCGCCGCGAAGCGATTCGTGAGTCGCTGCTCTACGACCAGGAAGGGGCGGACGTGCTGCCGGTCAATCCGGCCGGCGCGTACCTTGATGTGATCGCCGa                                                                                     <  2:1005422/148‑1 (MQ=255)
                                                                            cTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGc                                                                                        >  2:1178153/1‑137 (MQ=255)
                                                                            cTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGc                                                                                        <  1:1178153/137‑1 (MQ=255)
                                                                             tGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGcgccg                                                                          >  2:344167/1‑150 (MQ=255)
                                                                                     gACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGAt                                                                                                >  1:726936/1‑120 (MQ=255)
                                                                                     gACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGAt                                                                                                <  2:726936/120‑1 (MQ=255)
                                                                                      aCCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCcgc                                                                              >  1:1816943/1‑137 (MQ=255)
                                                                                         gCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGTCGAGCAgcga                                                                                                                                             <  2:1347193/71‑4 (MQ=255)
                                                                                         gCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGTCGAGCAgcga                                                                                                                                             >  1:1347193/1‑68 (MQ=255)
                                                                                         gCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGCGCCGCTTCGCGCCTGc                                                              >  2:2165332/1‑150 (MQ=255)
                                                                                         gCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGCGCCGCTTCGCGCCTGc                                                              >  2:2004976/1‑150 (MQ=255)
                                                                                              aGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGCGCCGCTTCGCGCCTGCCACTg                                                         >  1:1111557/1‑150 (MQ=255)
                                                                                               gCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTAc                                                                                                      >  1:513460/1‑104 (MQ=255)
                                                                                               gCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTAc                                                                                                      <  2:513460/104‑1 (MQ=255)
                                                                                                ttACCAGATGAACCCGATGAACCGCCGCGAAGCGTTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCGAGCCGGCCGGTGCGTACCTTGATGTGATCGCCTACATCCGCGCCGCTTCGCGCCTGCCACTggc                                                       <  2:292057/149‑1 (MQ=255)
                                                                                                 tACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGCGCCGCTTCGCGCCTGCCACTggcg                                                      >  2:458702/1‑150 (MQ=255)
                                                                                                 tACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGAGCGCCGACATCCGCGCCGCTTCGCGCCTGCCACTggcg                                                      <  1:1875788/150‑1 (MQ=255)
                                                                                                   ccAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGCGCCGCTTCGCGCCTGCCACTggcggc                                                    <  1:1591160/150‑1 (MQ=255)
                                                                                                     aGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGt                                                                                                            >  1:2136498/1‑92 (MQ=255)
                                                                                                     aGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGt                                                                                                            <  2:2136498/92‑1 (MQ=255)
                                                                                                      gATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGAt                                                                                                <  1:61112/103‑1 (MQ=255)
                                                                                                      gATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGCGCCGCTTCGCGCCTGCCACTggt                                                       >  1:1342754/1‑143 (MQ=255)
                                                                                                            cccGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGCGCCGCTTCGCGCCTGCCACTGGCGGCGTACCAGGt                                           >  1:1925840/1‑150 (MQ=255)
                                                                                                             ccGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGAt                                                                                                <  1:568870/96‑1 (MQ=255)
                                                                                                             ccGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGAt                                                                                                >  2:568870/1‑96 (MQ=255)
                                                                                                             ccGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGCGCCGCTTCGCGCCTGCCACTGGCGGCGTACCAGGTg                                          >  2:146284/1‑150 (MQ=255)
                                                                                                              cGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGCGCCGCTTCGCGCCTGCCACTGGCGGCGTACCAGGTGa                                         <  1:2352562/150‑1 (MQ=255)
                                                                                                               gATGAACCGCCGCGAAGCGGTGCGTGAGTCTCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGCGCCGCTTCGCGCCTGCCACTGGCGGCGTACCAGGTGAg                                        >  1:2362702/1‑150 (MQ=255)
                                                                                                                    aCCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGCGCCGCTTCGCGCCTGCCACTGGCGGCGTACCAGGTGAGTGGCg                                   >  2:1638140/1‑150 (MQ=255)
                                                                                                                     ccgccgCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGAt                                                                                                >  1:1083619/1‑88 (MQ=255)
                                                                                                                     ccgccgCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGCGCCGCTTCGCGCCTGCCACTGGCGGCGTACCAGGTGAGTGGCGa                                  <  2:950880/150‑1 (MQ=255)
                                                                                                                      cgccgcGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGCGCCGCTTCGCGCCTGCCACTGGCGGCGTACCAGGAGAGTGGCGAg                                 >  2:2315815/1‑150 (MQ=255)
                                                                                                                        ccgcGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGCGCCGCTTCGCGCCTGCCACTGGCGGCGTACCAGGTGAGTGGCGAGTa                               >  2:1869142/1‑150 (MQ=255)
                                                                                                                        ccgcGAAGCGGTGCGTGAGTCGCTGCTCGAAGAGCAGGAAGGGGCGGCCGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTTATGTGATCGCCGACATCCGNGCCGCTTCGCGCCTGCCACTGGCGGCGTACCAGGTGAGTGGCGAGTa                               <  1:986707/150‑1 (MQ=255)
                                                                                                                          gcgAAGCGGTGCGTGAGTCGCTGCTCGACGACCAGAAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGCGCCGCTTCGCGCCTGCCACTGGCGGCGTACCAGGTGAGTGGCGAGTATg                             >  1:1666236/1‑150 (MQ=255)
                                                                                                                             aaGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGCGCCGCTTCGCGCCTGCCACTGGCGGCGTACCAGGTGAGTGGCGAGTATGCga                          <  1:940718/150‑1 (MQ=255)
                                                                                                                                  gTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGCGCCGCTTCGCGCCt                                                                >  2:2190087/1‑107 (MQ=255)
                                                                                                                                  gTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGCGCCGCTTCGCGCCt                                                                <  1:2190087/107‑1 (MQ=255)
                                                                                                                                  gTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGCGCCGCTTCGCGCCTGCCACTGGCGGCGTACCAGGTGAGTGGCGAGTATGCGATGATc                     >  2:535971/1‑150 (MQ=255)
                                                                                                                                   tGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGAAATCCAAGGTGTTTAGCTCCGTCCACGGCCGGCGTACCAGGTGAGTGGCGAGTATGCGATGATCa                    >  2:529755/1‑150 (MQ=255)
                                                                                                                                          gTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGCGCCGCTTCGCGCCTGCCACTGGCGGCGTACCAGGTGAGTGGCGAGTATGCga                          >  1:939084/1‑137 (MQ=255)
                                                                                                                                          gTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGCGCCGCTTCGCGCCTGCCACTGGCGGCGTACCAGGTGAGTGGCGAGTATGCga                          <  2:939084/137‑1 (MQ=255)
                                                                                                                                          cttcctgctCGTCGAGCAGGAAGGGGCGGACGTGCTGATGGTCAGTCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGCGCCGCTTCGCGCCTGCCACTGGCGGCGTAGCAGGTGAGTGGCGAGTATGCGATGATCAAGTTCgg             >  2:1290641/5‑150 (MQ=255)
                                                                                                                                           tCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGCGCCGCTTCGCGCCTGCCACTGGCGGCGTACCAGGTGAGTGGCGAGTATGCGATGATCAAGTTCGGt            <  1:1638140/150‑1 (MQ=255)
                                                                                                                                           tCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGCGCCGCTTCGCGCCTGCCACTGGCGGCGTACCAGGTGAGTGGCGAGTATGCGATGATCAAGTTCGGt            <  1:1112230/150‑1 (MQ=255)
                                                                                                                                             gttgctGGTCGACCAGGAAGGGGCGGCCGTGCTGATGGTAAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCCACATCCGCGCGGTTTCGCGCCTGCCACTGGCGGCGTACCAGGTGAGTGGCGAGTATGCGATAATCAAGTTCGGtgg          <  2:1111557/148‑1 (MQ=255)
                                                                                                                                                  tCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGCGCCGCTTCGCGACTGCCACTGGCGGCGTACCAGGTGAGTGGCGAGTATGCGATGATCAAGTTCGGtggcctgg     >  1:162618/1‑150 (MQ=255)
                                                                                                                                                    gacgaGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGCACCGCTTCGCGCCTGCCACTGGCGGCGTACCAGGTGAGTGGCGAGTATGCGATGATCAAGTTCGGtggcctggcc   >  1:772323/1‑150 (MQ=255)
                                                                                                                                                    gacgaGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCAGACATCCGCGCCGCTTCGCGCCTGCCACTGGCGGCGTACCAGGTGAGTGGCGAGTATGCGATGATCAAGTTCGGtggcctggcc   >  2:798601/1‑150 (MQ=255)
                                                                                                                                                     acgaGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGCGCCGCTTCGCGCCTGCCACTGGCGGCGTACCAGGTGAGTGGCGAGTATGCGATGATCAAGTTCGGTGGCCTGGCCg  <  1:1752557/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCGCGATCATGGCCTACTCGACCAAGTTCGCCTCGTCGTTGTATGGCCCGTTCCGTGAGGCCGGCGGCACCGCACTGAAGGGCGACCGCAAAAGCTACCAGATGAACCCGATGAACCGCCGCGAAGCGGTGCGTGAGTCGCTGCTCGACGAGCAGGAAGGGGCGGACGTGCTGATGGTCAAGCCGGCCGGTGCGTACCTTGATGTGATCGCCGACATCCGCGCCGCTTCGCGCCTGCCACTGGCGGCGTACCAGGTGAGTGGCGAGTATGCGATGATCAAGTTCGGTGGCCTGGCCG  >  NC_002947/3313428‑3313725

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: