Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,421,570 G→C 7.7% G127G (GGG→GGC) ‡ moaB‑I → molybdenum cofactor biosynthesis protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,421,5700GC 7.7% 285.9 / 9.5 90G127G (GGG→GGC) ‡moaB‑Imolybdenum cofactor biosynthesis protein B
Reads supporting (aligned to +/- strand):  ref base G (40/43);  new base C (5/2);  total (45/45)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

GCGGTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCGGCCTGTGAC  >  NC_002947/2421422‑2421717
                                                                                                                                                     |                                                                                                                                                    
gcggTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCggg                                                                                                                                                      <  8:144971/149‑1 (MQ=255)
  ggTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCggg                                                                                                                                                      >  2:86436/1‑147 (MQ=255)
   gTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTg                                                                                                                                                   <  6:381734/149‑1 (MQ=255)
   gTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTg                                                                                                                                                   <  1:251412/149‑1 (MQ=255)
     aCCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGc                                                                                                                                                 >  3:219609/1‑149 (MQ=255)
        gggTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCaa                                                                                                                                              >  5:357425/1‑149 (MQ=255)
        gggTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCaa                                                                                                                                              >  6:340708/1‑149 (MQ=255)
         ggTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAAc                                                                                                                                             >  6:315643/1‑149 (MQ=255)
         ggTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAAc                                                                                                                                             <  4:300874/149‑1 (MQ=255)
              aCTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCGTTGGCCAGCCCGGCCAAGGCAc                                                                                                                                        <  5:188181/148‑1 (MQ=255)
                tGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGc                                                                                                                                      >  2:273811/1‑149 (MQ=255)
                  gTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTg                                                                                                                                    >  7:376628/1‑149 (MQ=255)
                        aCAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGtgc                                                                                                                              >  4:386663/1‑149 (MQ=255)
                         cAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGtgct                                                                                                                             >  5:93657/1‑149 (MQ=255)
                          aGCACCCCCGAGGCGGTGGCCTGCGTGCTGGAGAAGCAAGTAGACGGCTTTGGGGCGCTGTGCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCCTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGtgctg                                                                                                                            <  3:126116/149‑1 (MQ=255)
                          aGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGGCTTGGCCGGGCGGGGCCACCGCACGCCGgtgttgtg                                                                                                                            >  2:204767/1‑145 (MQ=255)
                           gCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGGCCGCTGGTGtgctgc                                                                                                                           >  1:195126/1‑149 (MQ=255)
                               ccccGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGctgcctgc                                                                                                                       >  8:210858/1‑149 (MQ=255)
                                 ccGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGctgcctgc                                                                                                                       >  1:49138/1‑147 (MQ=255)
                                 ccGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCa                                                                                                                     <  8:177766/149‑1 (MQ=255)
                                   gtggCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAgg                                                                                                                   <  6:93657/147‑1 (MQ=255)
                                        ggTGGCCTGCCTGCTGGACAAGCAAGTAGAGGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTcca                                                                                                              <  5:75915/149‑1 (MQ=255)
                                        ggTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCTGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTcca                                                                                                              <  8:376478/149‑1 (MQ=255)
                                         gTGGCCGGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGGGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTccac                                                                                                             <  2:53036/149‑1 (MQ=255)
                                          tGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTccacc                                                                                                            <  3:309033/149‑1 (MQ=255)
                                             cctgcctgCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAAc                                                                                                        >  7:6575/1‑149 (MQ=255)
                                              ctgcctgcTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACg                                                                                                       >  1:177653/1‑149 (MQ=255)
                                              ctgcctgcTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACg                                                                                                       >  7:337827/1‑149 (MQ=255)
                                                    gctgGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGGGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGc                                                                                                 <  1:220001/149‑1 (MQ=255)
                                                     ctgGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCtgtg                                                                                                  <  6:357425/147‑1 (MQ=255)
                                                        gACAAGCAAGTAGACTGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCtgt                                                                                                   >  8:315719/1‑143 (MQ=255)
                                                        gACAAGCAAGTAGACTGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCtgt                                                                                                   <  7:315719/143‑1 (MQ=255)
                                                          caagcaagTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCg                                                                                           >  2:288739/1‑149 (MQ=255)
                                                          caagcaagTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCg                                                                                           <  2:177653/149‑1 (MQ=255)
                                                           aagcaagTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAAAGGATGTGCGTACCg                                                                                           >  1:369041/1‑149 (MQ=255)
                                                             gcaagTAGACGTCTTTGGGGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTTGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGtt                                                                                        <  6:323696/149‑1 (MQ=255)
                                                              caagTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGc                                                                                                 >  7:113631/1‑139 (MQ=255)
                                                              caagTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGc                                                                                                 <  8:113631/139‑1 (MQ=255)
                                                               aattcgACGGCTTTGGGTCGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTgg                                                                                      <  6:137657/144‑1 (MQ=255)
                                                                  tAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGAt                                                                                   >  3:187028/1‑149 (MQ=255)
                                                                    gACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATgg                                                                                 <  8:6575/149‑1 (MQ=255)
                                                                    gACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATg                                                                                  <  6:213850/148‑1 (MQ=255)
                                                                    gACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATg                                                                                  <  7:259270/148‑1 (MQ=255)
                                                                      cGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCa                                                                               <  1:86436/149‑1 (MQ=255)
                                                                       ggCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCAt                                                                              >  2:127108/1‑149 (MQ=255)
                                                                          tttGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCt                                                                           >  6:416974/1‑149 (MQ=255)
                                                                            tGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTg                                                                          <  1:127108/148‑1 (MQ=255)
                                                                              gTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGc                                                                        <  3:386663/148‑1 (MQ=255)
                                                                               tGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCg                                                                      <  8:337827/149‑1 (MQ=255)
                                                                                 aGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTgg                                                                                      >  6:113424/1‑131 (MQ=255)
                                                                                 aGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTgg                                                                                      <  5:113424/131‑1 (MQ=255)
                                                                                 aGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGag                                                                    >  4:362772/1‑149 (MQ=255)
                                                                                 aGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGag                                                                    <  2:50436/149‑1 (MQ=255)
                                                                                    tGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGagcag                                                                 >  7:240742/1‑149 (MQ=255)
                                                                                    tGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCAGGCCGagcag                                                                 >  3:365896/1‑149 (MQ=255)
                                                                                         cGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGa                                                            >  1:316244/1‑149 (MQ=255)
                                                                                          ggCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTTCCTTGGCCGGGCTGGCCAACCGCACCCTGGTGTGCTGCCAGCCCGGCTACACCAACGCTGTGCGTACCCGGTGGGATGGCATACTGGCCGAGCAGCTCg                                                             >  4:179278/1‑147 (MQ=255)
                                                                                          ggCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGAt                                                           >  5:308700/1‑149 (MQ=255)
                                                                                             aGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTTGCCGGGCATGCCAACCGCCCGCGGGTGTTGTGCCAGCCAGGCGCCACCAACCCGTTGCGTACCCGTTGGGATGGCATGCTTGCCGAGCAGCACGaatcc                                                        >  2:367053/1‑145 (MQ=255)
                                                                                              gATTTCCGTGGCCGAATTTGGCACCTCGACCGTGCAGTCGCGGGCCTTGGCCGGGGTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGccc                                                       <  5:315643/149‑1 (MQ=255)
                                                                                                 ttCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATg                                                                                  <  2:55209/119‑1 (MQ=255)
                                                                                                 ttCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCGCCACCAACGCTGTGCGTACCGGTTGGGATg                                                                                  >  1:55209/1‑119 (MQ=255)
                                                                                                    cGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCAt                                                  <  2:49138/148‑1 (MQ=255)
                                                                                                    cGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATc                                                 <  1:288739/149‑1 (MQ=255)
                                                                                                     gTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGTCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCCAGCAGCTCGATGCCCGTCATCg                                                <  7:271479/149‑1 (MQ=255)
                                                                                                        gCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGccc                                             >  1:261503/1‑149 (MQ=255)
                                                                                                        gCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGcc                                              >  8:1780/1‑148 (MQ=255)
                                                                                                               ttGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAAc                                      >  5:363824/1‑149 (MQ=255)
                                                                                                                    aCCTCCACCGTGCAGTCGCGTGCCTTTGCCGGGCGGGCCAACCGCACCCGGGTGTGCTGCCTGCCAGGCTCCACCAACCCGGTGCGTACCGGTTGGGGTGGCATCCAGGCCCACCAGCAAGAAGCCCGTCATAGCCCATGCAAcata                                   >  6:30338/1‑144 (MQ=255)
                                                                                                                    aCCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTc                                   >  7:278719/1‑147 (MQ=255)
                                                                                                                    aCCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTc                                   <  5:30338/147‑1 (MQ=255)
                                                                                                                    aCCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTc                                   <  8:278719/147‑1 (MQ=255)
                                                                                                                     ccTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCATGCTCCACCAACGCTGTGCGTACCGGTTGGTATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTg                                <  3:238306/149‑1 (MQ=255)
                                                                                                                     ccTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTg                                <  4:219609/149‑1 (MQ=255)
                                                                                                                     ccTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTg                                <  3:362772/149‑1 (MQ=255)
                                                                                                                           ccGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGccc                                                       >  2:361261/1‑120 (MQ=255)
                                                                                                                           ccGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGccc                                                       <  1:361261/120‑1 (MQ=255)
                                                                                                                           ccGTGCAGTCGCGTGCCATGGCCGGGCTGGCCCACGGCACGCGGGTGTGCTGCCTGCCACGCTCCACCAACCCTGTGCGTACCGGGTGGGGTGGCCTCCTGGCCGCGCCGCTCGATTCCCGGCATCGCCCTTGCAACTTCGTGgcgccc                          >  2:383228/1‑147 (MQ=255)
                                                                                                                                 aGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAg                    >  7:17136/1‑149 (MQ=255)
                                                                                                                                  ccagcgTGCCGTTGGCCAGCCCGGCCAAGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAg                    >  1:236247/4‑148 (MQ=38)
                                                                                                                                   tCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAgca                  <  8:376628/149‑1 (MQ=255)
                                                                                                                                     gcgTGCCGTTGGCCAGCCCGGCCAAGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAgcagg                >  1:284604/1‑149 (MQ=255)
                                                                                                                                        tgccgttGGCCAGCCCGGCCAAGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAgc                   >  4:422232/6‑143 (MQ=255)
                                                                                                                                        tgccgttGGCCAGCCCGGCCAAGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAgc                   <  3:422232/138‑1 (MQ=255)
                                                                                                                                         gccgttGGCCAGCCCGGCCAAGGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAgcaggcag             >  1:149421/5‑149 (MQ=255)
                                                                                                                                          ccTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGc            <  8:240742/148‑1 (MQ=255)
                                                                                                                                          ccTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCg           >  1:282811/1‑149 (MQ=255)
                                                                                                                                           cgttGGCCAGCCCGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCg           >  8:237891/3‑149 (MQ=255)
                                                                                                                                           cTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCgg          <  1:273811/149‑1 (MQ=255)
                                                                                                                                                gCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCGGCCtg      <  2:261503/149‑1 (MQ=255)
                                                                                                                                                 ccGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCGGCCtgt     <  1:204767/149‑1 (MQ=255)
                                                                                                                                                 ccGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCGGCCtgt     <  2:316244/149‑1 (MQ=255)
                                                                                                                                                    ggCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCGGccgggcct  <  6:363824/149‑7 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCGGTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCGGCCTGTGAC  >  NC_002947/2421422‑2421717

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