Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,165,540 C→T 100% L133F (CTC→TTC)  gtsC → mannose/glucose ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,165,5400CT99.0% 1023.6 / ‑2.1 286L133F (CTC→TTC) gtsCmannose/glucose ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base T (145/137);  total (145/141)
Fisher's exact test for biased strand distribution p-value = 1.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 286 total aligned reads displayed.

GCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCCTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATCCCCGACGCGCTGGTCAAGGCCGCGCGGTTGGACGGTGCC  >  NC_002947/1165393‑1165686
                                                                                                                                                   |                                                                                                                                                  
gcgGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTc                                                                                                                                                  <  1:388941/150‑1 (MQ=255)
    tGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCa                                                                                                                                              >  2:2301954/1‑150 (MQ=255)
      cTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAg                                                                                                                                            >  2:4420846/1‑150 (MQ=255)
        gATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCt                                                                                                                                          >  2:1841340/1‑150 (MQ=255)
         aTCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCtt                                                                                                                                         >  1:1748905/1‑150 (MQ=255)
          tctcCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTtg                                                                                                                                        <  2:3617773/150‑1 (MQ=255)
           ctcCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTtgg                                                                                                                                       >  2:3877543/1‑150 (MQ=255)
            tcCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTtggc                                                                                                                                      <  1:4422097/150‑1 (MQ=255)
               aCCGCGGTTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTtggcctg                                                                                                                                   <  1:4114588/150‑1 (MQ=255)
               aCCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTtggcctg                                                                                                                                   <  2:778859/150‑1 (MQ=255)
               aCCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTtggcctg                                                                                                                                   >  2:5145813/1‑150 (MQ=255)
               aCCGCGATGGGTGCACTGAACGGCTAGGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTtggcctg                                                                                                                                   >  2:266584/1‑150 (MQ=255)
                 cgcgATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTtggcctggc                                                                                                                                 >  1:5188406/1‑150 (MQ=255)
                 cgcgATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTtggcctggc                                                                                                                                 >  1:4208041/1‑150 (MQ=255)
                 cgcgATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTtggcctggc                                                                                                                                 <  1:88814/150‑1 (MQ=255)
                 cgcgATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTtggcctggc                                                                                                                                 <  1:3505936/150‑1 (MQ=255)
                   cgATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCa                                                                                                                               <  1:1706511/150‑1 (MQ=255)
                         gTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCa                                                                                                                         >  1:722667/1‑150 (MQ=255)
                         gTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCa                                                                                                                         >  1:849157/1‑150 (MQ=255)
                                gAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACtggtgg                                                                                                                  >  1:104156/1‑150 (MQ=255)
                                 aaCGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGc                                                                                                                 <  1:2955652/150‑1 (MQ=255)
                                 aaCGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGc                                                                                                                 >  2:5512436/1‑150 (MQ=255)
                                 aaCGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGc                                                                                                                 <  2:2808276/150‑1 (MQ=255)
                                      cTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCctggt                                                                                                            >  1:4595609/1‑150 (MQ=255)
                                       tATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCCTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCctggtg                                                                                                           <  1:3628619/150‑1 (MQ=255)
                                        aTGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCctggtgc                                                                                                          >  1:2012147/1‑150 (MQ=255)
                                           tgCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCctggtgctgg                                                                                                       <  2:314207/150‑1 (MQ=255)
                                            gCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCctggtgctggt                                                                                                      <  1:5437231/150‑1 (MQ=255)
                                              gcacGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCtggtgctggtgc                                                                                                    <  2:5208940/147‑1 (MQ=255)
                                                tCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCAt                                                                                                  <  1:3849632/150‑1 (MQ=255)
                                                tCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCAt                                                                                                  <  2:4290373/150‑1 (MQ=255)
                                                tCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCGGGTGCAt                                                                                                  >  2:678162/1‑150 (MQ=255)
                                                tCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTACCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCAt                                                                                                  <  2:2604433/150‑1 (MQ=255)
                                                     gtgGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGt                                                                                             <  1:3964070/150‑1 (MQ=255)
                                                      tgGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTc                                                                                            >  2:5103775/1‑150 (MQ=255)
                                                       ggCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCt                                                                                           <  1:1888991/150‑1 (MQ=255)
                                                        gcgcTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTa                                                                                          >  1:5335150/1‑150 (MQ=255)
                                                          gcTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACg                                                                                        >  1:1062856/1‑150 (MQ=255)
                                                          gcTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACg                                                                                        >  1:1262219/1‑150 (MQ=255)
                                                          gcTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACg                                                                                        >  1:3067747/1‑150 (MQ=255)
                                                            ttCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACggc                                                                                      <  2:1639199/150‑1 (MQ=255)
                                                            ttCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACggc                                                                                      >  2:1500603/1‑150 (MQ=255)
                                                              ccGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACggcct                                                                                    <  2:1913929/150‑1 (MQ=255)
                                                               cGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACggcctg                                                                                   <  1:5470457/150‑1 (MQ=255)
                                                                 aGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACggcctggc                                                                                 >  1:4936355/1‑150 (MQ=255)
                                                                   gTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACggcctggcct                                                                               >  2:2543484/1‑150 (MQ=255)
                                                                     tCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTc                                                                             <  1:4732386/150‑1 (MQ=255)
                                                                      cGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTca                                                                            >  2:859332/1‑150 (MQ=255)
                                                                       gCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTcac                                                                           >  1:4723582/1‑150 (MQ=255)
                                                                        cAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTcacc                                                                          <  2:4595609/150‑1 (MQ=255)
                                                                         aaCTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTcacca                                                                         >  1:4859828/1‑150 (MQ=255)
                                                                          aCTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGGGGTCTACGGCCTGGCCTTcaccac                                                                        <  2:3385049/150‑1 (MQ=255)
                                                                            tGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGc                                                                      >  2:5411335/1‑150 (MQ=255)
                                                                                 ttcGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCTGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGttc                                                                 >  1:2715450/1‑150 (MQ=255)
                                                                                 ttcGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGttc                                                                 >  2:4153924/1‑150 (MQ=255)
                                                                                 ttcGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGttc                                                                 >  2:104130/1‑150 (MQ=255)
                                                                                 ttcGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGttc                                                                 <  2:3134397/150‑1 (MQ=255)
                                                                                 ttcGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGttc                                                                 >  1:159128/1‑150 (MQ=255)
                                                                                   cGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTtgg                                                                                                                                       <  2:2913189/78‑1 (MQ=255)
                                                                                   cGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTtgg                                                                                                                                       >  1:2913189/1‑78 (MQ=255)
                                                                                   cGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGttctt                                                               <  1:2301954/150‑1 (MQ=255)
                                                                                     gCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTcc                                                             <  1:4027677/150‑1 (MQ=255)
                                                                                     gCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTcc                                                             >  1:4154486/1‑150 (MQ=255)
                                                                                     gCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTcc                                                             <  1:598550/150‑1 (MQ=255)
                                                                                             ctgTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTc                                                     >  2:3036947/1‑150 (MQ=255)
                                                                                                ttCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTAc                                                  <  2:1140378/150‑1 (MQ=255)
                                                                                                  cGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGt                                                >  1:3476573/1‑150 (MQ=255)
                                                                                                   ggCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTg                                               >  1:4708914/1‑150 (MQ=255)
                                                                                                   ggCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTg                                               >  1:2540997/1‑150 (MQ=255)
                                                                                                    gCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGa                                              <  2:3844860/150‑1 (MQ=255)
                                                                                                          tCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATcc                                        >  1:1301416/1‑150 (MQ=255)
                                                                                                          tCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATcc                                        <  2:3774339/150‑1 (MQ=255)
                                                                                                           ccTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATccc                                       >  1:1154931/1‑150 (MQ=255)
                                                                                                            cTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATcccc                                      >  1:5540890/1‑150 (MQ=255)
                                                                                                            cTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATcccc                                      <  2:2121938/150‑1 (MQ=255)
                                                                                                            cTGCCGGTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATcccc                                      <  1:3388362/150‑1 (MQ=255)
                                                                                                             tGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTTACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATCCCCg                                     >  2:4932912/1‑150 (MQ=255)
                                                                                                             tGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATCCCCg                                     >  2:192726/1‑150 (MQ=255)
                                                                                                             tGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATCCCCg                                     >  2:4035813/1‑150 (MQ=255)
                                                                                                                 gTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATCCCCGAcgc                                 <  1:4100917/150‑1 (MQ=255)
                                                                                                                 gTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATCCCCGAcgc                                 >  2:4228120/1‑150 (MQ=255)
                                                                                                                   tCCAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCctggtgctggtg                                                                                                     >  2:5276437/1‑80 (MQ=255)
                                                                                                                     cAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATCCCCGACGCGCTg                             >  2:4670856/1‑150 (MQ=255)
                                                                                                                     cAGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATCCCCGACGCGCTg                             <  1:5114991/150‑1 (MQ=255)
                                                                                                                      aGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATCCCCGACGCGCTgg                            >  2:575853/1‑150 (MQ=255)
                                                                                                                      aGACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATCCCCGACGCGCTgg                            >  1:2126225/1‑150 (MQ=255)
                                                                                                                       gACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATCCCCGACGCGCTGGt                           <  1:2176674/150‑1 (MQ=255)
                                                                                                                       gACCGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATCCCCGACGCGCTGGt                           <  1:5007090/150‑1 (MQ=255)
                                                                                                                          cGTGCTGCTGCCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATCCCCGACGCGCTGGTCaa                        >  1:5267270/1‑150 (MQ=255)
                                                                                                                                gctgcCGGCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATCCCCGACGCGCTGGTCAAGGCcgc                  >  2:3279870/1‑150 (MQ=255)
                                                                                                                                     cGGCCTCGTTCACCCTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATCCCCGACGCGCTGGTCAAGGCCGCGCGGt             <  2:4698230/150‑1 (MQ=255)
                                                                                                                                      ggCCTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATCCCCGACGCGCTGGTCAAGGCCGCGCGGtt            >  2:4977313/1‑150 (MQ=255)
                                                                                                                                        ccTCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATCCCCGACGCGCTGGTCAAGGCCGCGCGGTTgg          >  2:4361864/1‑150 (MQ=255)
                                                                                                                                          tCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCGACTTCTACGTGAGCATCCCCGACGCGCTGGTCAAGGCCGCGCGGTTGGAc        <  1:5238614/150‑1 (MQ=255)
                                                                                                                                          tCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATCCCCGACGCGCTGGTCAAGGCCGCGCGGTTGGAc        >  2:13880/1‑150 (MQ=255)
                                                                                                                                          tCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATCCCCGACGCGCTGGTCAAGGCCGCGCGGTTGGAc        >  2:1598983/1‑150 (MQ=255)
                                                                                                                                          tCGTTCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATCCCCGACGCGCTGGTCAAGGCCGCGCGGTTGGAc        >  1:3162704/1‑150 (MQ=255)
                                                                                                                                             ttCACCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATCCCCGACGCGCTGGTCAAGGCCGCGCGGTTGGACGGt     >  1:1469252/1‑150 (MQ=255)
                                                                                                                                                aCCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATCCCCGACGCGCTGGTCAAGGCCGCGCGGTTGGACGGTGcc  <  2:5188406/150‑1 (MQ=255)
                                                                                                                                                aCCTTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATCCCCGACGCGCTGGTCAAGGCCGCGCGGTTGGACGGTGcc  <  2:1735899/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCCTCGGCAAGCTTGGCCTGGCCAGTACCACTGGTGGCCTGGTGCTGGTGCATGTGGTCTACGGCCTGGCCTTCACCACGCTGTTCTTCCGCAACTTCTACGTGAGCATCCCCGACGCGCTGGTCAAGGCCGCGCGGTTGGACGGTGCC  >  NC_002947/1165393‑1165686

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: