Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,421,572 T→C 7.8% L128P (CTG→CCG)  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,5720TC 7.8% 246.2 / 9.0 92L128P (CTG→CCG) moaB‑Imolybdenum cofactor biosynthesis protein B
Reads supporting (aligned to +/- strand):  ref base T (37/45);  new base C (5/2);  total (45/47)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.59e-01

GTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCGGCCTGTGACAG  >  NC_002947/2421425‑2421719
                                                                                                                                                    |                                                                                                                                                    
gTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTg                                                                                                                                                     <  6:381734/149‑1 (MQ=255)
gTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTg                                                                                                                                                     <  1:251412/149‑1 (MQ=255)
  aCCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGc                                                                                                                                                   >  3:219609/1‑149 (MQ=255)
     gggTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCaa                                                                                                                                                >  6:340708/1‑149 (MQ=255)
     gggTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCaa                                                                                                                                                >  5:357425/1‑149 (MQ=255)
      ggTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAAc                                                                                                                                               <  4:300874/149‑1 (MQ=255)
      ggTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAAc                                                                                                                                               >  6:315643/1‑149 (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                                                                                                     <  7:315719/143‑1 (MQ=255)
                                                     gACAAGCAAGTAGACTGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCtgt                                                                                                     >  8:315719/1‑143 (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                                                                                    <  7:259270/148‑1 (MQ=255)
                                                                 gACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATg                                                                                    <  6:213850/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                                     <  8:278719/147‑1 (MQ=255)
                                                                                                                 aCCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTc                                     >  7:278719/1‑147 (MQ=255)
                                                                                                                 aCCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTc                                     <  5:30338/147‑1 (MQ=255)
                                                                                                                  ccTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCATGCTCCACCAACGCTGTGCGTACCGGTTGGTATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTg                                  <  3:238306/149‑1 (MQ=255)
                                                                                                                  ccTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTg                                  <  3:362772/149‑1 (MQ=255)
                                                                                                                  ccTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTg                                  <  4:219609/149‑1 (MQ=255)
                                                                                                                        ccGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGccc                                                         <  1:361261/120‑1 (MQ=255)
                                                                                                                        ccGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGccc                                                         >  2:361261/1‑120 (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       <  2:316244/149‑1 (MQ=255)
                                                                                                                                              ccGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCGGCCtgt       <  1:204767/149‑1 (MQ=255)
                                                                                                                                                 ggCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCGGccgggcct    <  6:363824/149‑7 (MQ=255)
                                                                                                                                                   cTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCGGCCTGTGACAg  <  5:416974/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCGGCCTGTGACAG  >  NC_002947/2421425‑2421719

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: