Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,935,087 G→C 5.0% L103L (CTG→CTC htrB → lipid A biosynthesis lauroyl acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,935,0870GC 5.0% 1082.8 / 32.6 318L103L (CTG→CTChtrBlipid A biosynthesis lauroyl acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (153/149);  new base C (4/12);  total (157/161)
Fisher's exact test for biased strand distribution p-value = 6.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 327 total aligned reads displayed.

CGTCGGCGCATCGCCGCGCGCAACCTGGAGCTGTGCTTCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGCGCCGCGCTGCTCGGCCAGCAGCACACCATCGACGGCATGTACC  >  NC_002947/1934938‑1935235
                                                                                                                                                     |                                                                                                                                                    
cGTCGGCGCATCGCCGCGCGCAACCTGGAGCTGTGCTTCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTg                                                                                                                                                      <  2:4911227/150‑1 (MQ=255)
 gTCGGCGCATCGCCGCGCGCAACCTGGAGCTGTGCTTCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTgg                                                                                                                                                     >  2:7944265/1‑150 (MQ=255)
       gcATCGCCGCGCGCAACCTGGAGCTGTGCTTCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGcc                                                                                                                                               <  2:118983/150‑1 (MQ=255)
        cATCGCCGCGCGCAACCTGGAGCTGTGCTTCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCt                                                                                                                                              >  1:2069486/1‑150 (MQ=255)
         aTCGCCGCGCGCAACCTGGAGCTGTGCTTCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTg                                                                                                                                             >  2:2741310/1‑150 (MQ=255)
          tCGCCGCGCGCAACCTGGAGCTGTGCTTCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTgg                                                                                                                                            <  1:46087/150‑1 (MQ=255)
             ccgcGCGCAACCTGGAGCTGTGCTTCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGccc                                                                                                                                         >  1:8197375/1‑150 (MQ=255)
              cgcgcgCAACCTGGAGCTGTGCTTCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCa                                                                                                                                        >  1:9140549/1‑150 (MQ=255)
              cgcgcgCAACCTGGAGCTGTGCTTCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCa                                                                                                                                        <  1:3381118/150‑1 (MQ=255)
                cgcgcAACCTGGAGCTGTGCTTCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCata                                                                                                                                      <  2:3777012/150‑1 (MQ=255)
                cgcgcAACCTGGAGCTGTGCTTCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCGCGCCAGCCTGGccttgg                                                                                                                                      <  2:5344776/150‑5 (MQ=255)
                 gcgcAACCTGGAGCTGTGCTTCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCatat                                                                                                                                     >  1:1533794/1‑150 (MQ=255)
                    cAACCTGGAGCTGTGCTTCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGa                                                                                                                                  <  2:6598643/150‑1 (MQ=255)
                    cAACCTGGAGCTGTGCTTCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGa                                                                                                                                  <  2:5957940/150‑1 (MQ=255)
                     aaCCTGGAGCTGTGCTTCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAg                                                                                                                                 <  1:4238019/150‑1 (MQ=255)
                     aaCCTGGAGCTGTGCTTCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAg                                                                                                                                 <  1:3776531/150‑1 (MQ=255)
                      aCCTGGAGCTGTGCTTCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAgg                                                                                                                                <  1:2650085/150‑1 (MQ=255)
                           gAGCTGTGCTTCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCtt                                                                                                                           <  2:6037973/150‑1 (MQ=255)
                           gAGCTGTGCTTCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCtt                                                                                                                           <  2:2472604/150‑1 (MQ=255)
                           gAGCTGTGCTTCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCtt                                                                                                                           <  2:41218/150‑1 (MQ=255)
                           gAGCTGTGCGTCCCGGAGCAGTCGGTCGACGAACCGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAAGGGCCATGAGCTGGTGGGGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCAGATCGAGGGCCtt                                                                                                                           <  2:1486041/150‑1 (MQ=255)
                            aGCTGTGCTTCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTg                                                                                                                          <  1:4183883/150‑1 (MQ=255)
                             gCTGTGCTTCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGa                                                                                                                         <  2:6649743/150‑1 (MQ=255)
                             gCTGTGCTTCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGa                                                                                                                         <  2:6706672/150‑1 (MQ=255)
                             gCTGTGCTTCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGa                                                                                                                         <  2:4363993/150‑1 (MQ=255)
                               tgtgCTTCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGc                                                                                                                       >  2:9005245/1‑150 (MQ=255)
                                 tgCTTCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCAc                                                                                                                     <  1:6454140/150‑1 (MQ=255)
                                   cTTCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCt                                                                                                                   >  1:1166528/1‑150 (MQ=255)
                                    ttCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTg                                                                                                                  <  1:216860/150‑1 (MQ=255)
                                               gtcggtcgACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCca                                                                                                       >  1:4132771/1‑150 (MQ=255)
                                               gtcggtcgACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCca                                                                                                       >  2:5031100/1‑150 (MQ=255)
                                               gtcggtcgACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCca                                                                                                       <  1:8922023/150‑1 (MQ=255)
                                                tcggtcgACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCcag                                                                                                      >  1:2986761/1‑150 (MQ=255)
                                                      gacgaACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCg                                                                                                >  2:1703196/1‑150 (MQ=255)
                                                      gacgaACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCg                                                                                                >  1:1003891/1‑150 (MQ=255)
                                                       acgaACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCctgccccgcc                                                                                               >  2:860235/1‑144 (MQ=255)
                                                       acgaACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCgg                                                                                               <  2:3967529/150‑1 (MQ=255)
                                                       acgaACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCgg                                                                                               <  2:5215311/150‑1 (MQ=255)
                                                       acgaACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCgg                                                                                               <  1:642747/150‑1 (MQ=255)
                                                       acgaACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCgg                                                                                               <  1:6513161/150‑1 (MQ=255)
                                                       acgaACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCgg                                                                                               <  1:8065427/150‑1 (MQ=255)
                                                         gaACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCgggg                                                                                             >  2:9122217/1‑150 (MQ=255)
                                                           aCGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGcc                                                                                           >  1:8469700/1‑149 (MQ=255)
                                                           aCGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCa                                                                                           <  1:3852112/150‑1 (MQ=255)
                                                             ggCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAgg                                                                                         >  1:3293512/1‑150 (MQ=255)
                                                                aGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGACATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGgcg                                                                                      <  2:36554/150‑1 (MQ=255)
                                                                 gCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGgcgc                                                                                     <  2:8251508/150‑1 (MQ=255)
                                                                   gTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGa                                                                                   >  1:248418/1‑150 (MQ=255)
                                                                    tttGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGAt                                                                                  <  2:331743/150‑1 (MQ=255)
                                                                       gCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCt                                                                               >  1:2283312/1‑150 (MQ=255)
                                                                         tCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCa                                                                             >  1:613108/1‑150 (MQ=255)
                                                                         tCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCa                                                                             >  1:4062233/1‑150 (MQ=255)
                                                                         tCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCa                                                                             >  1:366820/1‑150 (MQ=255)
                                                                         tCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCa                                                                             >  2:4279921/1‑150 (MQ=255)
                                                                         tCAAGGCCAACTTCGCCTCTACCGGCATCGCCTCCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAGGGCCAGGCGCGCCAGCCTGGccttgg                                                                                                                                      >  2:5834704/1‑89 (MQ=37)
                                                                             ggCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGc                                                                         >  2:4775600/1‑150 (MQ=255)
                                                                             ggCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGc                                                                         >  2:2983784/1‑150 (MQ=255)
                                                                                 aaCTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTg                                                                     >  2:1920499/1‑150 (MQ=255)
                                                                                 aaCTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTg                                                                     <  2:4401762/150‑1 (MQ=255)
                                                                                    ttCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCAc                                                                  >  2:2070284/1‑150 (MQ=255)
                                                                                       gCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTc                                                               >  1:5265161/1‑150 (MQ=255)
                                                                                        ccTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCa                                                              >  1:7100740/1‑150 (MQ=255)
                                                                                         ctctACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCAc                                                             >  1:4412882/1‑150 (MQ=255)
                                                                                              ccGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTAccc                                                        >  2:8188513/1‑150 (MQ=255)
                                                                                              ccGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTAccc                                                        <  1:1644620/150‑1 (MQ=255)
                                                                                                ggCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTg                                                      <  2:7100740/150‑1 (MQ=255)
                                                                                                         ttcttcGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGgc                                             <  1:4922608/150‑1 (MQ=255)
                                                                                                          tcttcGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGgcg                                            >  2:5679499/1‑150 (MQ=255)
                                                                                                           cttcGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGgcgc                                           <  1:310273/150‑1 (MQ=255)
                                                                                                             tcGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGcgccg                                         >  1:904396/1‑150 (MQ=255)
                                                                                                             tcGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGcgccg                                         <  1:7976557/150‑1 (MQ=255)
                                                                                                                 aaTGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGACATCGGCGCCGCgct                                     <  2:6872991/150‑1 (MQ=255)
                                                                                                                 aaTGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGCGCCGCgct                                     <  2:3119984/150‑1 (MQ=255)
                                                                                                                 aaTGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGCGCCGCgct                                     >  2:1235637/1‑150 (MQ=255)
                                                                                                                  aTGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCCCTTCACTACCCTGGAAATCGGCGCCGCgctg                                    >  1:7684956/1‑150 (MQ=255)
                                                                                                                      ccATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGCGCCGCGCTGCTCg                                <  2:5889075/150‑1 (MQ=255)
                                                                                                                      ccATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGCGCCGCGCTGCTCg                                >  1:4737143/1‑150 (MQ=255)
                                                                                                                       cATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGCGCCGCGCTGCTCgg                               >  2:4557908/1‑150 (MQ=255)
                                                                                                                       cATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTACAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGCGCCGCGCTGCTCgg                               >  2:3430118/1‑150 (MQ=255)
                                                                                                                        aTGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGCGCCGCGCTGCTCGGc                              <  1:3339805/150‑1 (MQ=255)
                                                                                                                        aTGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGCGCCGCGCTGCTCGGc                              <  1:2250432/150‑1 (MQ=255)
                                                                                                                           aGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGCGCCGCGCTGCTCGGCcag                           <  2:739245/150‑1 (MQ=255)
                                                                                                                              tggtggtggCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCCGGGCGCGATCCTCATGGCGGTGCCCTTCACTACCCGGGGCCGCGGCGCCGCGCTGCTCGGCcagcgg                        >  2:9027977/1‑148 (MQ=255)
                                                                                                                              tggtggtggCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGCGCCGCGCTGCTCGGCcagcag                        >  2:1669563/1‑150 (MQ=255)
                                                                                                                               ggtggtggCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGCGCCGCGCTGCTCGGCCagcagc                       <  2:7983077/150‑1 (MQ=255)
                                                                                                                                gtggtggCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGCGCCGCGCTGCTCGGCCAgcagca                      >  2:7468017/1‑150 (MQ=255)
                                                                                                                                 tggtggCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGCGCCGCGCTGCTCGGCCAGCAGcac                     <  2:7569520/150‑1 (MQ=255)
                                                                                                                                     ggCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGCGCCGCGCTGCTCGGCCAGCAGCACACCa                 <  1:1267447/150‑1 (MQ=255)
                                                                                                                                     gcagcatggCCAGGCGCGCCAGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGCGCCGCGCTGCTCGGCCAGCAGCACACCa                 <  2:2935369/143‑1 (MQ=255)
                                                                                                                                       cccAAGGCCAGGCTGGCGCGCCTGGCCCATCTCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGCGCCGCGCTGCTCGGCCAGCAGCACACCATc               <  2:8092399/150‑1 (MQ=255)
                                                                                                                                       cccAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGCGCCGCGCTGCTCGGCCAGCAGCACACCATc               <  1:1768946/150‑1 (MQ=255)
                                                                                                                                        ccAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGCGCCGCGCTGCTCGGCCAACAGCACACCATCg              >  1:2392377/1‑150 (MQ=255)
                                                                                                                                         cAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGCGCCGCGCTGCTCGGCCAGCAGCACACCATCGa             >  1:6858124/1‑150 (MQ=255)
                                                                                                                                          aaGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGCGCCGCGCTGCTCGGCCAGCAGCACACCATCGAc            <  1:5940434/150‑1 (MQ=255)
                                                                                                                                          aaGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGCGCCGCGCTGCTCGGCCAGCAGCACACCATCGAc            <  1:4509222/150‑1 (MQ=255)
                                                                                                                                            ggCCAGGCGCGCCAGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGCGCCGCGCTGCTCGGCCAGCAGCACACCATCGACgg          >  1:7072501/1‑150 (MQ=255)
                                                                                                                                              ccAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGCGCCGCGCTGCTCGGCCAGCAGCACACCATCGACGGCa        >  1:8703019/1‑150 (MQ=255)
                                                                                                                                              ccAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGCGCCGCGCTGCTCGGCCAGCAGCACACCATCGACGGCa        >  2:4142977/1‑150 (MQ=255)
                                                                                                                                               cAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGCGCCGCGCTGCTCGGCCAGCAGCACACCATCGACGGCAt       >  2:8421245/1‑150 (MQ=255)
                                                                                                                                                  gCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGCGCCGCGCTGCTCGGCCAGCAGCACACCATCGACGGCATGTa    >  1:3999298/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGTCGGCGCATCGCCGCGCGCAACCTGGAGCTGTGCTTCCCGGAGCTGTCGGTCGACGAACGGCAGCGTTTGCTCAAGGCCAACTTCGCCTCTACCGGCATCGCCTTCTTCGAAATGGCCATGAGCTGGTGGTGGCCCAAGGCCAGGCTGGCGCGCCTGGCCCATATCGAGGGCCTTGAGCACCTGCAGGCCGCCCAGCAGGCGGGGCAGGGCGCGATCCTCATGGCGGTGCACTTCACTACCCTGGAAATCGGCGCCGCGCTGCTCGGCCAGCAGCACACCATCGACGGCATGTACC  >  NC_002947/1934938‑1935235

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