Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 93,195 G→A 6.8% R133R (CGC→CGT trpB ← tryptophan synthase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294793,1950GA 6.8% 1004.5 / 40.3 325R133R (CGC→CGTtrpBtryptophan synthase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (162/140);  new base A (11/11);  total (174/151)
Fisher's exact test for biased strand distribution p-value = 8.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 339 total aligned reads displayed.

GCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTT  >  NC_002947/93047‑93344
                                                                                                                                                    |                                                                                                                                                     
gCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAgc                                                                                                                                                      >  1:6330804/1‑150 (MQ=255)
 cAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAgcg                                                                                                                                                     <  2:4951753/150‑1 (MQ=255)
  agcgCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAAGGCAGGCCGAAacgg                                                                                                                                                    <  2:3067658/147‑2 (MQ=255)
  agcgCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAAGGCAGGCCGAAacgg                                                                                                                                                    <  1:6462275/147‑2 (MQ=255)
    ggCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAgcgcgc                                                                                                                                                  >  1:666353/1‑150 (MQ=255)
      cTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGTGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCgcgg                                                                                                                                                <  2:7221731/150‑1 (MQ=255)
       ttCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCgcggc                                                                                                                                               <  2:5618624/150‑1 (MQ=255)
       ttCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCgcggc                                                                                                                                               <  2:7106204/150‑1 (MQ=255)
           ttCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACa                                                                                                                                           >  2:3123040/1‑150 (MQ=255)
            tCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAg                                                                                                                                          <  2:3390805/150‑1 (MQ=255)
                 gCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtggcg                                                                                                                                     >  2:2687058/1‑150 (MQ=255)
                 gCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtggcg                                                                                                                                     >  1:5003196/1‑150 (MQ=255)
                   gTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtggcggt                                                                                                                                   <  1:4052337/150‑1 (MQ=255)
                      tttCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtggcggtggc                                                                                                                                <  2:7901989/150‑1 (MQ=255)
                      tttCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtggcggtggc                                                                                                                                >  2:3322675/1‑150 (MQ=255)
                      tttCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtggcggtggc                                                                                                                                <  1:2204617/150‑1 (MQ=255)
                        tCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGa                                                                                                                              <  1:6564655/150‑1 (MQ=255)
                          aGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGaca                                                                                                                            >  1:7014016/1‑150 (MQ=255)
                            ggTACCGGTGCTGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACAcc                                                                                                                          >  1:4130722/1‑150 (MQ=255)
                                      ccGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCg                                                                                                                >  2:3920997/1‑150 (MQ=255)
                                       cggcggTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCgg                                                                                                               >  2:4834710/1‑150 (MQ=255)
                                       cggcggTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCgg                                                                                                               >  2:249806/1‑150 (MQ=255)
                                            gTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCg                                                                                                          <  2:8069077/150‑1 (MQ=255)
                                                  ggCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCg                                                                                                    <  1:3122614/150‑1 (MQ=255)
                                                   gCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCgg                                                                                                   >  1:5886897/1‑150 (MQ=255)
                                                    cACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGc                                                                                                  <  2:6745159/150‑1 (MQ=255)
                                                     tttacttcGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCcgcg                                                                                                 >  1:3475872/7‑150 (MQ=255)
                                                     aCGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCg                                                                                                 >  1:7996352/1‑150 (MQ=255)
                                                     aCGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCg                                                                                                 >  1:8808671/1‑150 (MQ=255)
                                                     aCGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCg                                                                                                 >  1:3180612/1‑150 (MQ=255)
                                                     aCGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCg                                                                                                 >  2:5646191/1‑150 (MQ=255)
                                                        aTCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATc                                                                                              >  2:667998/1‑150 (MQ=255)
                                                         tctcGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGc                                                                                                                  <  2:5193663/129‑1 (MQ=255)
                                                          ttcGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGAtgat                                                                                            <  2:7582739/149‑4 (MQ=255)
                                                          ttcGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGAtgat                                                                                            >  1:8210192/2‑147 (MQ=255)
                                                          ctcGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAg                                                                                            <  2:2413890/150‑1 (MQ=255)
                                                           tcGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAgg                                                                                           >  2:3725444/1‑150 (MQ=255)
                                                           tcGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAgg                                                                                           >  2:4479406/1‑150 (MQ=255)
                                                           tcGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAgg                                                                                           >  2:5784786/1‑150 (MQ=255)
                                                           tcGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAgg                                                                                           >  1:1831253/1‑150 (MQ=255)
                                                            cGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGgc                                                                                          >  2:1599125/1‑150 (MQ=255)
                                                              gcacccAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGATGATgcgg                                                                                        <  2:824640/147‑2 (MQ=255)
                                                                gcCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCtt                                                                                      >  2:4670041/1‑150 (MQ=255)
                                                                  ccAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCtttt                                                                                    <  1:6086820/150‑1 (MQ=255)
                                                                  ccAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCtttt                                                                                    >  1:7181447/1‑150 (MQ=255)
                                                                  ccAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCtttt                                                                                    >  2:4428031/1‑150 (MQ=255)
                                                                     gcagcTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGc                                                                                 <  2:1471541/150‑1 (MQ=255)
                                                                          ttCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATg                                                                            >  1:963530/1‑150 (MQ=255)
                                                                           tCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATgc                                                                           <  1:5175214/150‑1 (MQ=255)
                                                                            cATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATgcg                                                                          <  2:680292/150‑1 (MQ=255)
                                                                              tGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCt                                                                        >  1:7878875/1‑150 (MQ=255)
                                                                               gCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCtt                                                                       >  2:1378032/1‑150 (MQ=255)
                                                                               gCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCtt                                                                       <  2:4870269/150‑1 (MQ=255)
                                                                               gCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCtt                                                                       <  1:65078/150‑1 (MQ=255)
                                                                               gCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCtt                                                                       >  1:840296/1‑150 (MQ=255)
                                                                                cGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTg                                                                      <  1:8327493/150‑1 (MQ=255)
                                                                                 ggAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTgg                                                                     >  2:2070435/1‑150 (MQ=255)
                                                                                  gAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGc                                                                    <  2:4332529/150‑1 (MQ=255)
                                                                                  gAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGc                                                                    <  2:6352882/150‑1 (MQ=255)
                                                                                  gAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGc                                                                    <  1:1030622/150‑1 (MQ=255)
                                                                                     cacGGTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCcag                                                                 <  2:5940887/150‑1 (MQ=255)
                                                                                       cGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCcagca                                                               <  1:1528476/150‑1 (MQ=255)
                                                                                         ttCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCagcagc                                                             >  2:5524638/1‑150 (MQ=255)
                                                                                             ccTGCTGGCGCACGATGTCGGTGGCACCAATGTAGATCAAGCAAGGCAGGCCGAAGCGCGCGGCGACCGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGATTCGGCGATCAGGCGCTTTTTGCCCATGCGCTAGGCCAGCAGCACtt                                                         >  2:1220727/1‑150 (MQ=255)
                                                                                                gctgGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGc                                                      >  2:2662860/1‑150 (MQ=255)
                                                                                                gctgGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGc                                                      <  1:2151621/150‑1 (MQ=255)
                                                                                                 ctgGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGcc                                                     >  2:8485544/1‑150 (MQ=255)
                                                                                                      gcTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACCGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGc                                                >  2:5359715/1‑150 (MQ=255)
                                                                                                      gcTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTGTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGc                                                >  1:2534379/1‑150 (MQ=255)
                                                                                                        tCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAg                                              <  2:1026822/150‑1 (MQ=255)
                                                                                                         caacgTCAGTGGCACCCATATAGACCACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGATGATGCGGGTTTTGCCCATGCGCTTGGCAAGCAGGATCTGGCCGAGGCaat                                             >  2:2308211/5‑148 (MQ=255)
                                                                                                          gATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAgtt                                            >  1:7659719/1‑150 (MQ=255)
                                                                                                              tcagTGGCACCCATATAGACCACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGATGATGCGGGTTTTGCCCATGCGCTTGGCAAGCAGGATCTGGCCGAGGCAattgtta                                        <  1:1237086/147‑2 (MQ=255)
                                                                                                                ggTGGCACCCCTGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACCGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCAc                                      >  2:6566462/1‑150 (MQ=255)
                                                                                                                      aCCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTtgtg                                >  1:5267566/1‑150 (MQ=255)
                                                                                                                      aCCCATATAGACCACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGATGATGCGGGTTTTGCCCATGCGCTTGGCAAGCAGGATCTGGCCGAGGCAATTGTTAATCTtgtg                                <  1:5969531/150‑1 (MQ=255)
                                                                                                                       cccATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTgc                               <  1:4215104/150‑1 (MQ=255)
                                                                                                                           tGTAGTTGACGCAAGGCAGGCCGATGCGCGCGGGGACTGTGGCGGTGGCGAGAGCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGGTTTTTGCGCATGCGCTTGGCCAGCAGCACTTGGCGGTTGCAGTTGTTCACCTTGTgcgcgc                           >  2:8823132/1‑150 (MQ=255)
                                                                                                                           tGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGacac                                                                                                                           <  1:4681908/54‑1 (MQ=255)
                                                                                                                           tGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTgcgcgc                           >  2:7662015/1‑150 (MQ=255)
                                                                                                                             tAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCg                         <  1:8544910/150‑1 (MQ=255)
                                                                                                                              aGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCgg                        >  2:7204913/1‑150 (MQ=255)
                                                                                                                               gATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGt                       >  1:3657402/1‑150 (MQ=255)
                                                                                                                                aTCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGtc                      >  2:6983783/1‑149 (MQ=255)
                                                                                                                                aTCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTa                      >  2:6144893/1‑150 (MQ=255)
                                                                                                                                 tCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTAt                     >  1:6247851/1‑150 (MQ=255)
                                                                                                                                 tCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTAt                     >  2:3497898/1‑150 (MQ=255)
                                                                                                                                 ccACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGATGATGCGGGTTTTGCCCATGCGCTTGGCAAGCAGGATCTGGCCGAGGCAATTGTTAATCTTGTGGCTGCCGGTAt                     >  1:7711197/2‑150 (MQ=255)
                                                                                                                                   aCGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATgg                   >  1:6047003/1‑150 (MQ=255)
                                                                                                                                     gCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGtt                 >  2:7773930/1‑150 (MQ=255)
                                                                                                                                      cAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTg                >  2:6910363/1‑150 (MQ=255)
                                                                                                                                       aaGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGa               <  1:3874330/150‑1 (MQ=255)
                                                                                                                                       aaGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGa               >  2:7316670/1‑150 (MQ=255)
                                                                                                                                          gcaggcCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGct            >  1:5643268/1‑150 (MQ=255)
                                                                                                                                             ggcCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCtt         <  1:3289138/150‑1 (MQ=255)
                                                                                                                                              gacGAATCGCGAGGCGAAAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCTCTTTCCGCGATAAGGCGCTATATGCCCATGCGCTTGGACAGAAGCACTTGGCCGATGCAGTTGTTCACCGTGTGCGCGCCGGAATGGTATAGCTCTTc        >  2:866559/3‑150 (MQ=255)
                                                                                                                                                 gAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACg     >  1:2714285/1‑150 (MQ=255)
                                                                                                                                                  aaGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGc    <  1:5547698/150‑1 (MQ=255)
                                                                                                                                                    gcgcgcGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCtt  <  2:1747828/150‑1 (MQ=255)
                                                                                                                                                    gcgcgcGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCtt  <  2:3263994/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTT  >  NC_002947/93047‑93344

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: