Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,676,202 C→G 7.2% intergenic (+22/‑38) hom → / → thrC homoserine dehydrogenase/threonine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,676,2020CG 7.2% 276.1 / 10.5 93intergenic (+22/‑38)hom/thrChomoserine dehydrogenase/threonine synthase
Reads supporting (aligned to +/- strand):  ref base C (43/42);  new base G (3/4);  total (46/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

CCAATGATCCTGGTGACCCATGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGC  >  NC_002947/1676058‑1676350
                                                                                                                                                |                                                                                                                                                     
ccAATGATCCTGGTGACCCATGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCcgcgc                                                                                                                                                  >  1:918532/1‑150 (MQ=255)
ccAATGATCCTGGTGACCCATGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCcgcgc                                                                                                                                                  <  2:717614/150‑1 (MQ=255)
  aaTGATCCTGGTGACGCATGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCGCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTGGAGCAGCTCAACTAACAAATATTGCCGACGAGCCGCCCGCGcgg                                                                                                                                                >  2:797038/1‑150 (MQ=255)
  aaTGATCCTGGTGACCCATGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGcgg                                                                                                                                                >  1:150531/1‑150 (MQ=255)
     gATCCTGGTGACCCATGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGcggccc                                                                                                                                             >  2:1861082/1‑150 (MQ=255)
        ccTGGTGACCCATGGCGTGGTCGAGCAGCGTCTCCATGACGCATAAGGCTTTCGCCGGGCCCGGCAGGACGTGGTCGGCAAGGTTGTGCGCCTTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGcggcccgtc                                                                                                                                          <  2:989377/150‑4 (MQ=255)
        ccTGGTGACCCATGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCGCGGGcggcccgtc                                                                                                                                          <  1:1307289/150‑4 (MQ=255)
        ccTGGTGACCCATGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGcggcccgtc                                                                                                                                          <  1:1255126/150‑4 (MQ=255)
         cTGGTGACCCATGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGcggccccggc                                                                                                                                         >  1:1620957/1‑150 (MQ=255)
         cTGGTGACCCATGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATGGCCGACGGGCCGCGCGGGcggcccgtcg                                                                                                                                         <  1:1209583/150‑5 (MQ=255)
             tGACCCATGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGCCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCa                                                                                                                                     <  2:1599215/150‑1 (MQ=255)
                 ccATGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCAt                                                                                                                                 >  1:1807470/1‑150 (MQ=255)
                  cATGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATc                                                                                                                                <  2:1587689/150‑1 (MQ=255)
                  cATGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAACAGCTCACATAACCAATTTTGTCGTGGATCAGCGCCTGCGGCGCCTGCCCCACCATc                                                                                                                                >  2:206306/1‑150 (MQ=255)
                      gCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGACGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCCTTCGCGCCGGGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAg                                                                                                                            <  2:1307851/150‑1 (MQ=255)
                           gTCGAGCAGCGTATCAATGACGCCATCGTCGCTATGGAAGCCCTGCAGGACGTGGGCGGCAAGGTTGTGCGCATTAGAGCTGGTATTCACAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTg                                                                                                                      >  1:432173/1‑150 (MQ=255)
                              gAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCAc                                                                                                                   >  1:434810/1‑150 (MQ=255)
                              gAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCAc                                                                                                                   <  2:1704023/150‑1 (MQ=255)
                                  agcGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCAt                                                                                                               <  2:1627036/150‑1 (MQ=255)
                                   gcGTATCAATTACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATg                                                                                                              <  2:794159/150‑1 (MQ=255)
                                    cGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATgc                                                                                                             <  2:1828918/150‑1 (MQ=255)
                                        tCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCta                                                                                                         >  1:448496/1‑150 (MQ=255)
                                         cAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCtat                                                                                                        >  1:273351/1‑150 (MQ=255)
                                          aaTGCCGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGATGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCtata                                                                                                       >  2:1120865/1‑150 (MQ=255)
                                          aaTGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGAGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGTGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCtata                                                                                                       >  1:566264/1‑150 (MQ=255)
                                           aTGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCTACGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTatat                                                                                                      <  1:206306/150‑1 (MQ=255)
                                           aTGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTatat                                                                                                      >  1:1233260/1‑150 (MQ=255)
                                             gACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCa                                                                                                    <  2:1620957/150‑1 (MQ=255)
                                             gACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCa                                                                                                    <  1:1160756/150‑1 (MQ=255)
                                                gCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTa                                                                                                 >  1:653800/1‑150 (MQ=255)
                                                gCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACATAAGCGCTATATCAgaa                                                                                                 >  2:1543409/1‑148 (MQ=255)
                                                gCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCCGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTa                                                                                                 >  2:189593/1‑150 (MQ=255)
                                                gCCATCCTCACTCTAAAAGCCCTGCACGACCTGCTCGGCAAGGCTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTa                                                                                                 >  2:21258/1‑150 (MQ=255)
                                                      gtcgCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGaa                                                                                                                             <  2:1233768/117‑1 (MQ=255)
                                                      gtcgCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGaa                                                                                                                             >  1:1233768/1‑117 (MQ=255)
                                                        cgCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGc                                                                                         <  1:541974/150‑1 (MQ=255)
                                                          ctctGGAAGCCCTGCAGGACGTGGTCGTCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCa                                                                                       <  1:1861082/150‑1 (MQ=255)
                                                          ctctGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCACGCGCGGCCCCGTCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCa                                                                                       <  2:1053697/150‑1 (MQ=255)
                                                           tctGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAg                                                                                      >  1:1884303/1‑150 (MQ=255)
                                                               gAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTTCCCGCGGCCAGGCAc                                                                                  <  1:1808429/150‑1 (MQ=255)
                                                               aaaGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGGTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGACCCGGCCCCAGCATCGAAGGTTTGCCCCCATGCGCTATATCAGTACCCGCGGCCAGGCAc                                                                                  <  1:76177/149‑1 (MQ=255)
                                                                aaGCCCTGCAGGACGTGGTCGGGAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGCGCCGCCCGCGCGGCCCCCGTTTAATCGAGAATTTGTACCACCCTTGCCCTATATCAGTACCCGCGGCCAGGCAcc                                                                                 >  2:880410/1‑150 (MQ=255)
                                                                 aGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCTTCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCg                                                                                <  1:346871/150‑1 (MQ=255)
                                                                 aGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCg                                                                                >  2:162896/1‑150 (MQ=255)
                                                                   cccTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGc                                                                              >  1:29602/1‑150 (MQ=255)
                                                                   cccTGCAGGACGCGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCcgcg                                                                                                                                                   <  1:1327831/82‑1 (MQ=255)
                                                                   cccTGCAGGACGCGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCcgcg                                                                                                                                                   >  2:1327831/1‑82 (MQ=255)
                                                                    ccTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGcc                                                                             <  1:1898867/150‑1 (MQ=255)
                                                                     cTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTTCCGACGGTCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGccc                                                                            >  1:233186/1‑150 (MQ=255)
                                                                          ggACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCa                                                                                                    >  1:1704765/1‑121 (MQ=255)
                                                                          ggACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCa                                                                                                    <  2:1704765/121‑1 (MQ=255)
                                                                           gACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAAtt                                                                      <  2:1665496/150‑1 (MQ=255)
                                                                            aCGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAAttt                                                                     <  2:273351/150‑1 (MQ=255)
                                                                             cGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGCGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTc                                                                    <  2:150531/150‑1 (MQ=255)
                                                                                ggTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGaa                                                                 <  2:653800/150‑1 (MQ=255)
                                                                                 gTCGGCAAGGTTGTGCGCATTCGCGTCGAACAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAg                                                                >  2:1875316/1‑150 (MQ=255)
                                                                                  tCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGa                                                               >  2:1869809/1‑150 (MQ=255)
                                                                                   cGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGAc                                                              <  1:880410/150‑1 (MQ=255)
                                                                                     gCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACCAATATTGCCGACGGGCCGCGCGggcggc                                                                                                                                               <  1:1484662/68‑1 (MQ=255)
                                                                                     gCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCGCGggcggc                                                                                                                                               >  2:1484662/1‑68 (MQ=255)
                                                                                       aaGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGtgc                                                          >  2:1827745/1‑150 (MQ=255)
                                                                                       aaGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGtgc                                                          >  2:261309/1‑150 (MQ=255)
                                                                                       aaGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCACTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGAAGtgg                                                          >  2:480690/1‑149 (MQ=255)
                                                                                        aGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAGTATTGCCGACGGGCCGCCCGCGCGGCCACGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGGCGtgct                                                         <  2:124901/150‑1 (MQ=255)
                                                                                         ggTTGTGCGCATTTGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGtgctg                                                        <  2:1917464/150‑1 (MQ=255)
                                                                                           ttGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTgctgct                                                      >  2:1340900/1‑150 (MQ=255)
                                                                                                   aTTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGGCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCTGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGctggcc                                              <  1:717207/150‑1 (MQ=255)
                                                                                                    ttCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGTGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGGACCCGCGGCCAGGCGCCGGCCCTGAATTTCGAAGACGTGCTGCTGGctggcct                                             <  2:566264/150‑1 (MQ=255)
                                                                                                       gcgTCGAGCAGCTCAACTAACAAATGTTGCCGCCGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGACAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGctggcctggc                                          <  1:434570/150‑1 (MQ=255)
                                                                                                         gTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTCTATCAGTACCCGCGGCCAGGCAACGGCCCTGAATTTCGAAGCAGTGCGGCTGGCTGGCCGGGCCa                                        >  2:315637/1‑150 (MQ=255)
                                                                                                            gAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCta                                                                                                         <  2:782950/82‑1 (MQ=255)
                                                                                                            gAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCta                                                                                                         >  1:782950/1‑82 (MQ=255)
                                                                                                                   tCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGAcggcgg                              <  1:1120865/150‑1 (MQ=255)
                                                                                                                   ccAACTAACCAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGAcggcgg                              <  2:1580658/149‑1 (MQ=255)
                                                                                                                    cAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACggcggc                             >  2:720825/1‑150 (MQ=255)
                                                                                                                     aaCTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCAcc                                                                                                                  <  1:836310/64‑1 (MQ=255)
                                                                                                                     aaCTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCAcc                                                                                                                  >  2:836310/1‑64 (MQ=255)
                                                                                                                        tAACAAACATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGt                         >  1:385958/1‑150 (MQ=255)
                                                                                                                         aaCAAATATTGCCGCCGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTa                        <  1:1133598/150‑1 (MQ=255)
                                                                                                                            aaaTATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGt                     >  2:364525/1‑150 (MQ=255)
                                                                                                                            aaaTATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGt                     >  1:301862/1‑150 (MQ=255)
                                                                                                                              atatTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGTGc                   <  1:1614153/150‑1 (MQ=255)
                                                                                                                              atatTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGTGc                   <  1:600925/150‑1 (MQ=255)
                                                                                                                                   gCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGTGCCGGaa              <  2:29602/150‑1 (MQ=255)
                                                                                                                                    ccGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATCTCAGTACCCGCGGGCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGGGCCGGaaa             >  1:483833/1‑150 (MQ=255)
                                                                                                                                    ccGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGTGCCGGaaa             >  2:590739/1‑150 (MQ=255)
                                                                                                                                      gACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGTGCCGGAAAAc           <  2:1233260/150‑1 (MQ=255)
                                                                                                                                      gACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGTGCCGGAAAAc           <  2:918532/150‑1 (MQ=255)
                                                                                                                                        ggggCCGCGCGGGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGTGCCGGaaaa            <  2:597816/146‑1 (MQ=255)
                                                                                                                                        ggggCCGCGCGGGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGTGCCGGaaaa            >  1:597816/2‑147 (MQ=255)
                                                                                                                                         gggCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGTGCCGGAAAACCTg        >  1:146184/1‑150 (MQ=255)
                                                                                                                                         gggCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGTGCCGGAAAACCTg        >  2:68100/1‑150 (MQ=255)
                                                                                                                                         gggCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGAACCCATGCGATATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGTGCCGGAAAACCTg        >  1:800054/1‑150 (MQ=255)
                                                                                                                                             cgcccgcGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGTCTATACGTGCCGGAAAACCTGCCAc    >  2:951013/1‑150 (MQ=255)
                                                                                                                                               cccgcGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGAc                                                              <  1:915960/90‑1 (MQ=255)
                                                                                                                                               cccgcGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGAc                                                              >  2:915960/1‑90 (MQ=255)
                                                                                                                                               cccgcGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGc  <  2:385958/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
CCAATGATCCTGGTGACCCATGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGC  >  NC_002947/1676058‑1676350

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