Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,718,691 G→A 8.1% intergenic (+35/+96) fabB → / ← PP_4176 3‑oxoacyl‑ACP synthase/5‑oxo‑L‑prolinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,718,6910GA 8.1% 249.8 / 8.9 87intergenic (+35/+96)fabB/PP_41763‑oxoacyl‑ACP synthase/5‑oxo‑L‑prolinase
Reads supporting (aligned to +/- strand):  ref base G (39/40);  new base A (4/3);  total (44/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.62e-01

CTGCGTAAAACCGAAGAAAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCCAGAAAATTCTCCCGGGCCACCTTGT  >  NC_002947/4718549‑4718839
                                                                                                                                              |                                                                                                                                                    
cTGCGTAAAACCGAAGAAAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCggg                                                                                                                                                <  8:222570/149‑1 (MQ=255)
     tAAAACCGAAGAAAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGtt                                                                                                                                           >  5:184257/1‑149 (MQ=255)
       aaaCCGAAGAAAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGAACCAGTCGGGGCGttt                                                                                                                                          <  3:296965/148‑1 (MQ=255)
        aaCCGAAGAAAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGttttt                                                                                                                                        <  6:184257/149‑1 (MQ=255)
          ccGAAGAAAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTtttgt                                                                                                                                      <  7:301111/149‑1 (MQ=255)
            gaagaaAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTtttgttt                                                                                                                                    >  8:148813/1‑149 (MQ=255)
              agaaAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTgg                                                                                                                                  <  6:145989/149‑1 (MQ=255)
                aaaaCGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGt                                                                                                                                <  3:40773/149‑1 (MQ=255)
                aaaaCGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGt                                                                                                                                >  3:232972/1‑149 (MQ=255)
                   aCGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACTCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAATTGAAAACGCCCCGACTGGTTCggg                                                                                                                                                <  7:146678/130‑1 (MQ=255)
                   aCGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACCCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTCCTGTAAATGAAAACGCCCCGACGCGTTCggg                                                                                                                                                >  8:146678/1‑130 (MQ=255)
                       cAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGATGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCtgt                                                                                                                         >  3:227855/1‑149 (MQ=255)
                        aaGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCtgtg                                                                                                                        >  5:47007/1‑149 (MQ=255)
                               acacGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGcc                                                                                                                  >  6:309167/1‑148 (MQ=255)
                               acacGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCct                                                                                                                 >  5:82085/1‑149 (MQ=255)
                                cacGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCctc                                                                                                                >  4:135741/1‑149 (MQ=255)
                                     tCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTcgc                                                                                                           <  6:260345/149‑1 (MQ=255)
                                     tCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTcgc                                                                                                           >  7:87465/1‑149 (MQ=255)
                                       aTGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCCAGCGGTGGGAAGGCATGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCgg                                                                                                         <  4:2662/149‑1 (MQ=255)
                                         gAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGTCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCggg                                                                                                        >  8:5716/1‑148 (MQ=255)
                                            cAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTaaa                                                                                                    >  1:156515/1‑149 (MQ=255)
                                             aagaGCTTCGGCTTCGGCGGCACCTACGCCACGGTGGTGCTCATGCGCTGGGAAGGCATGTAAGACGCTGCTGCTGTAAATGAATACGCCCCGACTGGGTGGGGGCGTTTTTGTTTGGGTTGCGTGTGCTGGCCGCTTCGCGGGTAAAc                                                                                                   <  8:66995/146‑1 (MQ=255)
                                             aaCAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAc                                                                                                   >  6:266835/1‑149 (MQ=255)
                                             aaCAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAc                                                                                                   <  1:157299/149‑1 (MQ=255)
                                             aaCAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAc                                                                                                   >  4:196947/1‑149 (MQ=255)
                                               cAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAccc                                                                                                 >  8:243669/1‑149 (MQ=255)
                                               cAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAccc                                                                                                 <  6:47007/149‑1 (MQ=255)
                                                 gcttcggcttcgGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGc                                                                                               <  7:5716/149‑1 (MQ=255)
                                                 gcttcggattcgGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGc                                                                                               <  2:342041/149‑1 (MQ=255)
                                                  cttcggcttcggCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGt                                                                                                                                            >  7:312409/1‑103 (MQ=255)
                                                  cttcggcttcggCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGt                                                                                                                                            <  8:312409/103‑1 (MQ=255)
                                                  cttcggcttcggCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGt                                                                                                                                            <  4:142835/103‑1 (MQ=255)
                                                  cttcggcttcggCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGt                                                                                                                                            >  3:142835/1‑103 (MQ=255)
                                                     cggcttcggcGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTcc                                                                                            >  5:265816/1‑148 (MQ=255)
                                                      ggcttcggcGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCt                                                                                           <  4:232972/148‑1 (MQ=255)
                                                        cttcggcGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCCGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTaca                                                                                        <  7:105471/149‑1 (MQ=255)
                                                           cggcggCACCAACGCCACCCTGGTCCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTAcaca                                                                                      <  7:243669/148‑1 (MQ=255)
                                                            ggcggcACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAg                                                                                     <  3:196947/148‑1 (MQ=255)
                                                                gcACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCAt                                                                                >  3:229490/1‑149 (MQ=255)
                                                                 cACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTc                                                                                             >  1:176520/1‑135 (MQ=255)
                                                                 cACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTc                                                                                             <  7:304247/135‑1 (MQ=255)
                                                                 cACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTc                                                                                             >  8:304247/1‑135 (MQ=255)
                                                                 cACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTc                                                                                             <  2:176520/135‑1 (MQ=255)
                                                                 cACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCAtt                                                                               <  3:135741/149‑1 (MQ=255)
                                                                    cAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTCTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTg                                                                             >  3:47578/1‑148 (MQ=255)
                                                                    cAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACCCACCGACTGTTTCGGGGCGTTATTGTTTGGGTTGCCTGTGCTGGCCTATTAGCGGCTAAACCAGCTGCTACACAGGTCATTTg                                                                             >  7:229674/1‑148 (MQ=255)
                                                                      aCGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAg                                                                          >  1:47818/1‑149 (MQ=255)
                                                                       cGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGt                                                                         <  5:309167/149‑1 (MQ=255)
                                                                        gCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTg                                                                        <  3:215468/149‑1 (MQ=255)
                                                                           aCCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCtt                                                                     >  2:160676/1‑149 (MQ=255)
                                                                           aCCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCtt                                                                     <  8:87465/149‑1 (MQ=255)
                                                                               tGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTtgt                                                                   <  7:7534/147‑1 (MQ=255)
                                                                               tGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTtgt                                                                   >  8:7534/1‑147 (MQ=255)
                                                                                 gTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTgg                                                                 <  4:229490/147‑1 (MQ=255)
                                                                                  tGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCa                                                              <  6:82085/149‑1 (MQ=255)
                                                                                     tCAAGCGCTGGGAAGGCAAGTAAGTCGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGGTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTg                                                           <  3:19305/149‑1 (MQ=255)
                                                                                     tCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAc                                                                                                   >  7:156155/1‑109 (MQ=255)
                                                                                     tCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAc                                                                                                   <  8:156155/109‑1 (MQ=255)
                                                                                       aaGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTg                                                           >  8:313038/1‑147 (MQ=255)
                                                                                           gcTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGc                                                      >  1:164946/1‑148 (MQ=255)
                                                                                               ggAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATg                                                 >  1:126621/1‑149 (MQ=255)
                                                                                                     cAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGgcgcgc                                           >  6:268944/1‑149 (MQ=255)
                                                                                                       aGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCt                                         >  8:88932/1‑149 (MQ=255)
                                                                                                          aaGAGGCTTCTGCTGTAATTGAATACGCCACGAGTGGTTCGTTGCGTTTTTGTTTGGTTTTCCTGTGCTGGCCTCTTCGAGGTTAAACCCGCTCCTACAAAGTTCTTTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTg                                       <  8:229674/148‑1 (MQ=255)
                                                                                                           aGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTggg                                     <  2:47818/149‑1 (MQ=255)
                                                                                                            gACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGc                                    <  5:266835/149‑1 (MQ=255)
                                                                                                              cGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCaa                                  <  2:164946/149‑1 (MQ=255)
                                                                                                               gctgctgctGTAACTGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAac                                 <  7:313038/149‑1 (MQ=255)
                                                                                                                    tgctgTAAATGAAAACGCCCCGACTCGTTCGGGGCGATTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTCGCCATTCCAGGCAAATGGCGCGCATTGGGCAACAgcgc                            >  8:274541/1‑149 (MQ=38)
                                                                                                                      ctgTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGc                                                      >  2:72260/1‑121 (MQ=255)
                                                                                                                      ctgTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGc                            <  4:20869/147‑1 (MQ=255)
                                                                                                                          aaaTGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCCAGaa                       >  6:7572/1‑148 (MQ=255)
                                                                                                                           aacaaaaaCGCCCCGAACCAGTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGc                                    >  4:21815/5‑134 (MQ=37)
                                                                                                                           aacaaaaaCGCCCCGAACCAGTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGc                                    <  3:21815/130‑1 (MQ=37)
                                                                                                                            aTGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCCAGAAAAt                    >  5:332262/1‑149 (MQ=255)
                                                                                                                               aaaaCGCCCCGAACCAGTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCtt                                        <  3:12315/126‑1 (MQ=38)
                                                                                                                               aaaaCGCCCCGAACCAGTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCtt                                        >  4:12315/1‑126 (MQ=38)
                                                                                                                               aaaaCGCCCCGAACCAGTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCAAGAAAATtct                 >  3:216182/1‑149 (MQ=38)
                                                                                                                                aaaCGCCCCGAACCAGTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCCAGAAAATtctc                >  2:87277/1‑149 (MQ=38)
                                                                                                                                    gCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAaca                                >  7:102319/1‑129 (MQ=255)
                                                                                                                                    gCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAaca                                <  8:102319/129‑1 (MQ=255)
                                                                                                                                    gCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCCAGAAAATTCTCCCgg            >  8:253975/1‑149 (MQ=255)
                                                                                                                                      cccGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCGTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCCAGAAAATTCTCCCGGGc          >  2:344792/1‑149 (MQ=255)
                                                                                                                                      cccGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCCAGAAAATTCTCCCGGGc          >  6:183085/1‑149 (MQ=255)
                                                                                                                                         gACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCCAGAAAATTCTCCCGGGCCAc       <  1:160676/149‑1 (MQ=255)
                                                                                                                                          aCTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCCAGAAAATTCTCCCGGGCCAcc      <  1:344792/149‑1 (MQ=255)
                                                                                                                                          aCTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCCAGAAAATTCTCCCGGGCCAcc      <  4:47578/149‑1 (MQ=255)
                                                                                                                                              gTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACGCCGCCAGAAAATTCTCCCGGGCCACCTtg   <  2:156515/148‑1 (MQ=255)
                                                                                                                                              gTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCCAGAAAATTCTCCCGGGCCACCTtgt  <  6:332262/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
CTGCGTAAAACCGAAGAAAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCCAGAAAATTCTCCCGGGCCACCTTGT  >  NC_002947/4718549‑4718839

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: