Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,524,302 A→G 7.4% T336A (ACC→GCC)  murG → UDP‑N‑acetylglucosamine‑‑N‑acetylmuramyl‑ (pentapeptide) pyrophosphoryl‑undecaprenol N‑acetylglucosamine transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,524,3020AG 7.4% 192.4 / 5.3 69T336A (ACC→GCC) murGUDP‑N‑acetylglucosamine‑‑N‑acetylmuramyl‑ (pentapeptide) pyrophosphoryl‑undecaprenol N‑acetylglucosamine transferase
Reads supporting (aligned to +/- strand):  ref base A (33/30);  new base G (2/3);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-01

CGATCACCAGACTCACAACGCCCAATATCTGGCCCGCGAAGGCGCTGCCTTCCTGATGCCACAAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCAGAAAGCCATGCCCCAGCCGAAGATGGGCCGTATTCGCCGCATTCACTTCGTCGGCATCGGCGGCGTGG  >  NC_002947/1524157‑1524450
                                                                                                                                                 |                                                                                                                                                    
cgaTCACCAGACTCACAACGCCCAATATCTGGCCCGCGAAGGCGCTGCCTTCCTGATGCCACAAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCg                                                                                                                                                   >  7:105424/1‑149 (MQ=255)
   tCACCAGACTCACAACGCCCAATATCTGGCCCGCGAAGGCGCTGCCTTCCTGATGCCACAAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTc                                                                                                                                                <  2:141645/149‑1 (MQ=255)
     aCCAGACTCACAACGCCCAATATCTGGCCCGCGAAGGCGCTGCCTTCCTGATGCCACAAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCgg                                                                                                                                              <  6:304867/149‑1 (MQ=255)
            tCACAACGCCCAATATCTGGCCCGCGAAGGCGCTGCCTTCCTGATGCCACAAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGCCAGGCGCCGAGCGGTGCCTgc                                                                                                                                       >  8:195582/1‑148 (MQ=255)
              acaACGCCCAATATCTGGCCCGCGAAGGCGCTGCCTTCCTGATGCCACAAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGc                                                                                                                                      >  6:343472/1‑148 (MQ=255)
              acaACGCCCAATATCTGGCCCGCGAAGGCGCTGCCTTCCTGATGCCACAAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGc                                                                                                                                      >  3:222133/1‑148 (MQ=255)
                 aCGCCCAATATCTGGCCCGCGAAGGCGCTGCCTTCCTGATGCCACAAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGCCAGGCGCCGAGCGGTGCCTgccatga                                                                                                                                  <  7:195582/149‑7 (MQ=255)
                 aCGCCCAATATCTGGCCAGCGAAGGCGCTGCCTTCCTGATGCCACAAGCGACTACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCTTGCCAGGCGCCGAGCGGTGCCTgccatga                                                                                                                                  <  1:313792/149‑7 (MQ=255)
                  cGCCCAATATCTGGCCCGCGAAGGCGCTGCCTTCCTGATGCCACAAGCGACAACTGGCCCAGCGCAACTCGCGGAACGCCTGATCGAGGTGCTGATGCAAGCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAAc                                                                                                                                 <  5:391044/149‑1 (MQ=255)
                      cAATATCTGGCCCGCGAAGGCGCTGCCTTCCTGATGCCACAAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACctgc                                                                                                                             >  2:11249/1‑149 (MQ=255)
                      cAATATCTGGCCCGCGAAGGCGCTGCCTTCCTGATGCCACAAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACctgc                                                                                                                             >  7:255939/1‑149 (MQ=255)
                      cAATATCTGGCCCGCGAAGGCGCTGCCTTCCTGATGCCACAAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACctgc                                                                                                                             <  7:288605/149‑1 (MQ=255)
                      cAATATCTGGCCCGCGAAGGCGCTGCCTTCCTGATGCCACAAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACctgc                                                                                                                             >  3:323859/1‑149 (MQ=255)
                      cAATATCTGGCCCGCGAAGGCGCTGCCTTCCTGATGCCACAAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACccgc                                                                                                                             >  7:340416/1‑149 (MQ=255)
                      cAATATCTGGCCCGCGAAGGCGCTGCCTTCCTGATGCCACAAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAAAGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGAACCGCTCGGCGCCTGGCCAACCctgc                                                                                                                             >  7:200339/1‑149 (MQ=255)
                          atCTGGCCCGCGAAGGGGCTGCCTTCCTGATGCCACAAGCGAGAACTGTCGCAGCGCTACTCGCTGAACGCCTTAACGAGGTGCTGATGCAACCCGAGAAAGTCATCGTCTTGGCAGGCACCGCTCGGCGCCTGGCCAAACCtgctgc                                                                                                                          <  8:280852/148‑1 (MQ=255)
                              ggCCCGCGAAGGCGCTGCCTTCCTGATGCCACAAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCGCCGAGCGGTGCCTgcca                                                                                                                                     <  5:10999/133‑4 (MQ=255)
                              ggCCCGCGAAGGCGCTGCCTTCCTGATGCCACAAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCGCCGAGCGGTGCCTgcca                                                                                                                                     >  6:10999/1‑130 (MQ=255)
                              ggCCCGCGAAGGCGCTGCCTTCCTGATGCCACAAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAAc                                                                                                                       >  7:291888/1‑147 (MQ=255)
                                    cgAAGGCGCTGCCTTCCTGATGCCACTATAGACAACGGGCTCAGCGCAACTCGCTGAACGACTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGTCTCCTGTCCAAACCTGCTGCAACCAGCACCg                                                                                                               <  8:65620/149‑1 (MQ=255)
                                        ggCGCTGCCTTCCTGATGCCACAAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGt                                                                                                           <  8:291888/149‑1 (MQ=255)
                                           gcTGCCTTCCTGATGCCACAAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGa                                                                                                        >  6:366968/1‑149 (MQ=255)
                                           gcTGCCTTCCTGATGCCACAAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGa                                                                                                        >  7:195633/1‑149 (MQ=255)
                                           gcTGCCTTCCTGATGCCACAAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCCGCTGCAACCAGCACCGTGGTCGa                                                                                                        <  8:340416/149‑1 (MQ=255)
                                                 ttCCTGATGCCACAAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTg                                                                                                  >  7:308023/1‑149 (MQ=255)
                                                 ttCCTGATGCCACAAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTg                                                                                                  >  5:282796/1‑149 (MQ=255)
                                                   ccTGATGCCACAAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGc                                                                                                 <  8:105424/148‑1 (MQ=255)
                                                    cTGATGCCACAAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCt                                                                                               >  1:370256/1‑149 (MQ=255)
                                                       aTGCCACAAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTg                                                                                              >  4:217844/1‑147 (MQ=255)
                                                          ccACAAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGcc                                                                                                                                     <  8:225854/105‑1 (MQ=255)
                                                          ccACAAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGcc                                                                                                                                     >  7:225854/1‑105 (MQ=255)
                                                           cacaAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTg                                                                                        >  2:121952/1‑149 (MQ=255)
                                                            acaAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTTGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTgg                                                                                       <  7:348217/149‑1 (MQ=255)
                                                            acaAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTgg                                                                                       <  4:222133/149‑1 (MQ=255)
                                                               aGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAg                                                                                           <  8:257232/142‑1 (MQ=255)
                                                               aGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAg                                                                                           >  7:257232/1‑142 (MQ=255)
                                                                gCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGcccg                                                                                   >  4:335670/1‑148 (MQ=255)
                                                                   aCAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATgg                                                                                <  8:195633/149‑1 (MQ=255)
                                                                    cAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGt                                                                               <  2:17408/149‑1 (MQ=255)
                                                                     aaCTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGtt                                                                              >  1:25148/1‑149 (MQ=255)
                                                                          gcgcagcgcaACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTgag                                                                           <  4:323859/147‑1 (MQ=255)
                                                                           cgcagcgcaACTCGCTGAACGCCTGAACGAGGTGCCGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGc                                                                        >  2:224039/1‑149 (MQ=255)
                                                                            gcagcgcaACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGc                                                                        >  3:114495/1‑148 (MQ=255)
                                                                             cagcgcaACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCa                                                                      <  2:25148/149‑1 (MQ=255)
                                                                                cgcaACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCAGaa                                                                   <  1:140289/149‑1 (MQ=255)
                                                                                cgcaACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCAGaa                                                                   <  3:299804/149‑1 (MQ=255)
                                                                                 gcaACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCAGaaa                                                                  <  8:255939/149‑1 (MQ=255)
                                                                                         cTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCAGAAAGCCATGcc                                                          >  2:315785/1‑149 (MQ=255)
                                                                                         cTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCAGAAAGCCATGcc                                                          >  6:357457/1‑149 (MQ=255)
                                                                                              cGCCTGAACGCGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCAGAAAGCCATGCCCCAGc                                                     <  8:308023/149‑1 (MQ=255)
                                                                                                    aaCGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCAGAAAGCCATGCCCCAGc                                                     <  4:63751/143‑1 (MQ=255)
                                                                                                    aaCGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCAGAAAGCCATGCCCCAGc                                                     >  3:63751/1‑143 (MQ=255)
                                                                                                     aCGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCAGAAAGCCATGCCCCAGCCGAAGAt                                              <  5:366968/149‑1 (MQ=255)
                                                                                                         ggTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCGATGGTTGAGAGCCAGAAAGCCATGCCCCAGCCGAAGATGGGc                                          >  4:195100/1‑149 (MQ=255)
                                                                                                            gCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCAGAAAGCCATGCCCCAGCCGAAGATGGGCCg                                        >  4:277301/1‑148 (MQ=255)
                                                                                                                aTGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCAGaaa                                                                  >  5:72383/1‑118 (MQ=255)
                                                                                                                aTGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCAGaaa                                                                  <  6:72383/118‑1 (MQ=255)
                                                                                                                   cAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCAGAAAGCCATGCCCCAGCCGAAGATGGGCCGTATTcgc                                 <  4:220432/148‑1 (MQ=255)
                                                                                                                    aaCCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCAGAAAGCCATGCCCCAGCCGAAGATGGGCCGTATTcgccg                               >  8:290098/1‑149 (MQ=255)
                                                                                                                       ccGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTg                                                                                        <  7:228049/89‑1 (MQ=255)
                                                                                                                       ccGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTg                                                                                        >  8:228049/1‑89 (MQ=255)
                                                                                                                         gagaAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCAGAAAGCCATGCCCCAGCCGAAGATGGGCCGTATTCGCCGCAtt                           >  2:108067/1‑148 (MQ=255)
                                                                                                                           gaAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCAGAAAGCCATGCCCCAGCCGAAGATGGGCCGTATTCGCCGCATTCAc                        <  5:378054/149‑1 (MQ=255)
                                                                                                                              aCTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCAGAAAGCCATGCCCCAGCCGAAGATGGGCCGTATTCGCCGCATTCACTtc                     >  2:343332/1‑149 (MQ=255)
                                                                                                                               cTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCAGAAAGCCATGCCCCAGCCGAAGATGGGCCGTATTCGCCGCATTCACTtcg                    <  1:11249/149‑1 (MQ=255)
                                                                                                                                tCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCAGAAAGCCATGCCCCAGCCGAAGATGGGCCGTa                                      >  2:266033/1‑130 (MQ=255)
                                                                                                                                tCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCAGAAAGCCATGCCCCAGCCGAAGATGGGCCGTa                                      <  1:266033/130‑1 (MQ=255)
                                                                                                                                  aaCGTCATGGCAGGCTCCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCAGAAAGCCATGCCCCAGCCGAAGATGGGCCGTATTCGCCGCATTCACTtcgtc                  <  3:217844/148‑1 (MQ=255)
                                                                                                                                    cGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCGTGGTTGAGAGCCAGAAAGCCATGCCCCAGCCGAAGATGGGCCGTATTCGCCGCATTCACTTCGTCGGc               <  3:335670/149‑1 (MQ=255)
                                                                                                                                     gtttggccaggcgCCGAGCGGTGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCAGAAAGCCATGCCCCAGCCGAAGATGGGCCGTATTCGCCGCATTCACTTCGTCGGCa              >  1:251389/8‑149 (MQ=255)
                                                                                                                                           gcaggcaCCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCAGAAAGCCATGCCCCAGCCGAAGATGGGCCGTATTCGCCGCATTCACTTCGTCGGCATcggcg        >  5:377621/1‑149 (MQ=255)
                                                                                                                                             aggcaCCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCAGAAAGCCATGCCCCAGCCGAAGATGGGCCGTATTc                                   >  7:46110/1‑120 (MQ=255)
                                                                                                                                             aggcaCCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCAGAAAGCCATGCCCCAGCCGAAGATGGGCCGTATTc                                   <  8:46110/120‑1 (MQ=255)
                                                                                                                                                caCCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCAGAAAGCCATGCCCCAGCCGAAGATGGGCCGTATTCGCCGCATTCACTTCGTCGGCATCGGCGGCGTg   <  1:224039/149‑1 (MQ=255)
                                                                                                                                                caCCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCAACGTCGTCGATATCTGCATGGAGGTGGCCCAAGGATGAGAGCCAGAAAGCCATGCCCCAGCCGAAGATGGGCCGTATTAGCCGCATTCCCTTCGTCGGCATCGGCGGCGTg   >  7:273486/1‑149 (MQ=255)
                                                                                                                                                 aCCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCAGAAAGCCATGCCCCAGCCGAAGATGGGCCGTATTCGCCGCATTCACTTCGTCGGCATCGGCGGCGTgg  <  2:370256/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CGATCACCAGACTCACAACGCCCAATATCTGGCCCGCGAAGGCGCTGCCTTCCTGATGCCACAAGCGACAACTGGCGCAGCGCAACTCGCTGAACGCCTGAACGAGGTGCTGATGCAACCCGAGAAACTCAACGTCATGGCAGGCACCGCTCGGCGCCTGGCCAAACCTGCTGCAACCAGCACCGTGGTCGATATCTGCCTGGAGGTGGCCCATGGTTGAGAGCCAGAAAGCCATGCCCCAGCCGAAGATGGGCCGTATTCGCCGCATTCACTTCGTCGGCATCGGCGGCGTGG  >  NC_002947/1524157‑1524450

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: