Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,865,699 T→A 6.1% M77L (ATG→TTG) ‡ purN ← phosphoribosylglycinamide formyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,865,6990TA 6.1% 1346.5 / 51.7 410M77L (ATG→TTG) ‡purNphosphoribosylglycinamide formyltransferase
Reads supporting (aligned to +/- strand):  ref base T (195/190);  new base A (8/17);  total (203/207)
Fisher's exact test for biased strand distribution p-value = 9.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.18e-01

Only 100 of 413 total aligned reads displayed.

GCGTATGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGCCCTT  >  NC_002947/1865552‑1865848
                                                                                                                                                    |                                                                                                                                                     
  gTATGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAgc                                                                                                                                                   >  2:5163864/1‑150 (MQ=255)
    aTGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAATCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAgcgc                                                                                                                                                 <  2:7481815/150‑1 (MQ=255)
    aTGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAgcgc                                                                                                                                                 >  2:2866802/1‑150 (MQ=255)
      gCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCgg                                                                                                                                               >  2:4299950/1‑150 (MQ=255)
       cAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGc                                                                                                                                              <  2:2745278/150‑1 (MQ=255)
        aGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCTAATGGcgcgcat                                                                                                                                             <  1:1907863/150‑4 (MQ=255)
        aGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCa                                                                                                                                             >  1:8261483/1‑150 (MQ=255)
        aGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCa                                                                                                                                             >  1:8079640/1‑150 (MQ=255)
          gCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATc                                                                                                                                           >  2:3992306/1‑150 (MQ=255)
          gCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATc                                                                                                                                           >  2:8118507/1‑150 (MQ=255)
          gCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATc                                                                                                                                           <  2:6972502/150‑1 (MQ=255)
             cTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCGAAGCCTTCGATGCCGCGCTAATGGCGC‑GCATCGac                                                                                                                                        <  1:96831/150‑2 (MQ=38)
             cTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGaa                                                                                                                                        <  2:3037892/150‑1 (MQ=255)
             cTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGaa                                                                                                                                        >  1:5113467/1‑150 (MQ=255)
                gTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGc                                                                                                                                     <  1:3146639/150‑1 (MQ=255)
                 tACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCt                                                                                                                                    >  1:5276713/1‑150 (MQ=255)
                   cTTGGGCAGCAGCGCCGGGTGGATGTTGAGCAGGCGGCCCTGATAGGGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCACAGCCGTCGATG‑CGCGCCCTTAGCGCGGCATCGAAGGCTTc                                                                                                                                  <  2:4921745/150‑1 (MQ=255)
                    ttGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGTCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCa                                                                                                                                 >  1:1386697/1‑150 (MQ=255)
                    ttGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCTTCGATGCCGCGCTAATGGCGC‑GCATCGACGGCtttg                                                                                                                                 <  2:8522558/150‑3 (MQ=38)
                    ttGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCa                                                                                                                                 >  1:3790283/1‑150 (MQ=255)
                    ttGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCa                                                                                                                                 <  1:2944911/150‑1 (MQ=255)
                    tagggCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGTCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCa                                                                                                                                 >  1:3605551/3‑150 (MQ=255)
                        gcagcaGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGc                                                                                                                             >  2:631891/1‑150 (MQ=255)
                           gcagcGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGt                                                                                                                          <  1:7790806/150‑1 (MQ=255)
                            cagcGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTc                                                                                                                         >  1:2794325/1‑150 (MQ=255)
                            cagcGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTc                                                                                                                         <  2:3232851/150‑1 (MQ=255)
                               cGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGaa                                                                                                                      >  2:5406671/1‑150 (MQ=255)
                                 aCGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACt                                                                                                                    <  1:4299950/150‑1 (MQ=255)
                                   gggTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTgg                                                                                                                  <  2:7864737/150‑1 (MQ=255)
                                   gggTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTgg                                                                                                                  <  2:2291261/150‑1 (MQ=255)
                                   gggTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTgg                                                                                                                  <  1:3915597/150‑1 (MQ=255)
                                    ggTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTggg                                                                                                                 >  1:6719771/1‑150 (MQ=255)
                                     gTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGgt                                                                                                                >  2:3146425/1‑150 (MQ=255)
                                        gATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGtgtg                                                                                                             >  2:7304573/1‑150 (MQ=255)
                                           gTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTc                                                                                                          >  2:6485735/1‑150 (MQ=255)
                                             tGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGa                                                                                                        <  1:6230417/150‑1 (MQ=255)
                                              gAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAg                                                                                                       <  1:6153786/150‑1 (MQ=255)
                                              gAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAg                                                                                                       <  1:1797209/150‑1 (MQ=255)
                                              gAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAg                                                                                                       <  1:8807583/150‑1 (MQ=255)
                                              gAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAg                                                                                                       >  1:3249434/1‑150 (MQ=255)
                                                    gcggcCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCacgac                                                                                                 >  2:126918/1‑150 (MQ=255)
                                                       gcCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCCCGGCCGTCGAACTGGGGGTGGTCGAGCCCGACGCt                                                                                              >  2:5904423/1‑150 (MQ=255)
                                                          cTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTc                                                                                           >  2:4037356/1‑150 (MQ=255)
                                                          cTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTc                                                                                           >  2:2579980/1‑150 (MQ=255)
                                                            gATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGa                                                                                         >  1:7242213/1‑150 (MQ=255)
                                                               aGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATgc                                                                                      <  2:188207/150‑1 (MQ=255)
                                                                 tGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATgccg                                                                                    >  2:5003116/1‑150 (MQ=255)
                                                                 tGGCGCACGAAGCCGCCACTGAGGAGGCGCATGAAACCGGCCAGCACCCCCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATgccg                                                                                    <  2:5760677/150‑1 (MQ=255)
                                                                   gcgcACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGa                                                                                                                       <  2:1950561/113‑1 (MQ=255)
                                                                     gcACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATgccggccg                                                                                >  1:5789546/1‑150 (MQ=255)
                                                                      cACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGc                                                                               >  2:8646010/1‑150 (MQ=255)
                                                                       aCGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCt                                                                              >  1:7030169/1‑150 (MQ=255)
                                                                         gAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTgc                                                                            <  1:4375011/150‑1 (MQ=255)
                                                                         gAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTgc                                                                            <  2:2433434/150‑1 (MQ=255)
                                                                              cgccACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGgcgc                                                                       <  2:7734614/150‑1 (MQ=255)
                                                                                 cACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGtt                                                                    >  1:3475802/1‑150 (MQ=255)
                                                                                  aCTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTg                                                                   >  2:4396055/1‑150 (MQ=255)
                                                                                  aCTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTg                                                                   >  1:3963818/1‑150 (MQ=255)
                                                                                      aGGATGCGCATGAACCCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAgg                                                               <  1:8162449/150‑1 (MQ=255)
                                                                                            cgcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAg                                                         <  2:3568422/150‑1 (MQ=255)
                                                                                                tGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCAt                                                     >  1:1321600/1‑150 (MQ=255)
                                                                                                tGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCAt                                                     >  1:2580009/1‑150 (MQ=255)
                                                                                                 gAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATc                                                    >  2:3467358/1‑150 (MQ=255)
                                                                                                  aaaCCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCg                                                   <  1:4871827/150‑1 (MQ=255)
                                                                                                   aaCCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCgg                                                  >  2:3349646/1‑150 (MQ=255)
                                                                                                       ggCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGcgcg                                              <  1:6422992/150‑1 (MQ=255)
                                                                                                        gCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCgg                                             >  2:6653623/1‑150 (MQ=255)
                                                                                                         ccAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGt                                            <  2:8272243/150‑1 (MQ=255)
                                                                                                          cAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGtt                                           <  2:5606535/150‑1 (MQ=255)
                                                                                                            gCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTgg                                         >  2:112899/1‑150 (MQ=255)
                                                                                                            gCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACGGGGTGTGGTCGAGCACCAAGCAGTAGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCAACGTCGCGGTTgg                                         >  1:2161940/1‑150 (MQ=255)
                                                                                                              accaccAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGaa                                       <  1:445721/150‑1 (MQ=255)
                                                                                                               ccaccaGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGaaa                                      <  2:4508439/150‑1 (MQ=255)
                                                                                                                caccaGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGTATCGGCGCGGTTGGAAAg                                     >  2:7180877/1‑150 (MQ=255)
                                                                                                                caccaGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAg                                     >  1:6603573/1‑150 (MQ=255)
                                                                                                                caccaGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAg                                     >  2:7910407/1‑150 (MQ=255)
                                                                                                                caccaGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAg                                     >  2:3840385/1‑150 (MQ=255)
                                                                                                                 accaGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGc                                    >  1:5359517/1‑150 (MQ=255)
                                                                                                                  ccaGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCa                                   >  1:8427497/1‑150 (MQ=255)
                                                                                                                  ccaGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCa                                   <  2:6572951/150‑1 (MQ=255)
                                                                                                                   caGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCTTGAGCCCGGCATCGAAGCCTTCACGGCCGACGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGACAGCAc                                  <  1:1499371/150‑1 (MQ=255)
                                                                                                                       tCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACcgcg                              >  2:8585642/1‑150 (MQ=255)
                                                                                                                         cGGGGCAACGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTTGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACcgcgcg                            <  2:2078226/150‑1 (MQ=255)
                                                                                                                          ggggCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCgg                           >  2:5240018/1‑150 (MQ=255)
                                                                                                                               aaaGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACg                      >  1:6989060/1‑150 (MQ=255)
                                                                                                                                  gCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCAc                   >  2:3617720/1‑150 (MQ=255)
                                                                                                                                   ccGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCAcc                  >  1:5213330/1‑150 (MQ=255)
                                                                                                                                    cGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCg                 <  1:7605809/150‑1 (MQ=255)
                                                                                                                                     ttCGATGCCGCGCTAATGGCGC‑GCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCgg                >  2:618750/2‑150 (MQ=255)
                                                                                                                                     gTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCgg                <  1:232195/150‑1 (MQ=255)
                                                                                                                                      tCGATGCCGCGCTAATGGCGC‑GCATCGACGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCggg               >  1:1013458/1‑150 (MQ=255)
                                                                                                                                           g‑cgcgcCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTc          <  1:8585642/150‑1 (MQ=255)
                                                                                                                                              gcgcCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCt        <  1:561443/150‑1 (MQ=255)
                                                                                                                                               cgcCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTg       >  1:5662972/1‑150 (MQ=255)
                                                                                                                                                 ccATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGcc     <  1:5464509/150‑1 (MQ=255)
                                                                                                                                                 ccATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGcc     >  1:6651636/1‑150 (MQ=255)
                                                                                                                                                 ccATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGcc     <  2:1322153/150‑1 (MQ=255)
                                                                                                                                                 ccATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGcc     <  2:1567371/150‑1 (MQ=255)
                                                                                                                                                    ttAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGCCCtt  <  1:3395233/150‑1 (MQ=255)
                                                                                                                                                    ttAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGCCCtt  <  1:685964/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCGTATGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGCCCTT  >  NC_002947/1865552‑1865848

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