Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,707,318 C→G 6.0% intergenic (‑20/+90) PP_4163 ← / ← PP_4164 hypothetical protein/alpha/beta hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,707,3180CG 6.0% 212.5 / 3.3 67intergenic (‑20/+90)PP_4163/PP_4164hypothetical protein/alpha/beta hydrolase
Reads supporting (aligned to +/- strand):  ref base C (31/32);  new base G (2/2);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

GTCGCCGGGCGTTGATCGACCACCTTGGCGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTG  >  NC_002947/4707174‑4707461
                                                                                                                                                |                                                                                                                                               
gTCGCCGGGCGTTGATCGACCACCTTGGCGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAgcg                                                                                                                                             >  1:199633/1‑149 (MQ=255)
 tCGCCGGGCGTTGATCGACCACCTTGGCGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAgcgc                                                                                                                                            >  3:292578/1‑149 (MQ=255)
      gggCGTTGATCGACCACCTTGGCGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTGAGCGCCgcct                                                                                                                                       <  7:253256/149‑4 (MQ=255)
          gTTGATCGACCACCTTGGCGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGaa                                                                                                                                   <  2:144128/149‑1 (MQ=255)
           ttGATCGACCACCTTGGCGTCCTGCGGGACCGGTAGGCGGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAg                                                                                                                                  <  8:253010/149‑1 (MQ=255)
            tGATCGACCACCTTGGCGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGc                                                                                                                                 >  2:27732/1‑149 (MQ=255)
                 gACCACCTTGGCGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAg                                                                                                                            >  8:24320/1‑149 (MQ=255)
                    caccTTGGCGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGAt                                                                                                                         >  6:171575/1‑149 (MQ=255)
                           gCGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTgg                                                                                                                   <  4:292578/148‑1 (MQ=255)
                           gCGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTgg                                                                                                                   <  2:199633/148‑1 (MQ=255)
                            cGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCa                                                                                                                 <  3:198124/149‑1 (MQ=255)
                              tCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCt                                                                                                                     >  4:7383/1‑143 (MQ=255)
                              tCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCt                                                                                                                     <  3:7383/143‑1 (MQ=255)
                                    gggACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATg                                                                                                         >  1:267942/1‑149 (MQ=255)
                                     ggACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATgg                                                                                                        <  7:197462/149‑1 (MQ=255)
                                      gACCGGTAGGCTGCCGGCCCACGACACCGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCCAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGt                                                                                                       >  4:293600/1‑149 (MQ=255)
                                       aCCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGtt                                                                                                      <  1:150651/149‑1 (MQ=255)
                                       aCCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGtt                                                                                                      <  3:131509/149‑1 (MQ=255)
                                         cGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGt                                                                                                    >  5:81076/1‑149 (MQ=255)
                                         cGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCACTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCa                                                                                                                 >  6:223791/1‑136 (MQ=255)
                                         cGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCACTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCa                                                                                                                 <  5:223791/136‑1 (MQ=255)
                                           gTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGtct                                                                                                  >  7:99917/1‑149 (MQ=255)
                                                cTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTg                                                                                             <  7:266092/149‑1 (MQ=255)
                                                 tGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGt                                                                                            >  4:31179/1‑149 (MQ=255)
                                                    cGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTTGCAGTTGCATGGTTGTCTCCCTGTTc                                                                                          >  7:121542/1‑148 (MQ=255)
                                                          aCGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGc                                                                                   >  2:17644/1‑149 (MQ=255)
                                                           cGACACGGGGGTGGCGATCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGcc                                                                                  <  7:24320/149‑1 (MQ=255)
                                                           cGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTCGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGcc                                                                                  <  8:121542/149‑1 (MQ=255)
                                                               acGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTc                                                                                          <  7:160160/137‑1 (MQ=255)
                                                               acGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTc                                                                                          >  8:160160/1‑137 (MQ=255)
                                                               acGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGc                                                                              >  1:199559/1‑149 (MQ=255)
                                                                    ggTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGa                                                                         <  1:27732/149‑1 (MQ=255)
                                                                    ggTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGa                                                                         <  5:171575/149‑1 (MQ=255)
                                                                     gTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGac                                                                        <  2:199559/149‑1 (MQ=255)
                                                                         cGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACtt                                                                    >  5:183199/1‑149 (MQ=255)
                                                                         cGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACtt                                                                    <  5:291193/149‑1 (MQ=255)
                                                                                caggcaggCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGt                                                             >  5:160963/1‑149 (MQ=255)
                                                                                caggcaggCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGt                                                             >  4:223502/1‑149 (MQ=255)
                                                                                   gcaggcagCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGc                                                           >  2:285598/1‑148 (MQ=255)
                                                                                     aggcagCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAg                                                        >  7:238570/1‑149 (MQ=255)
                                                                                       gcagCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGt                                                       >  2:263983/1‑148 (MQ=255)
                                                                                         agCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCaa                                                    <  6:81076/149‑1 (MQ=255)
                                                                                         agCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCaa                                                    <  2:267942/149‑1 (MQ=255)
                                                                                          gCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAg                                                   <  6:183199/149‑1 (MQ=255)
                                                                                          gCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAg                                                   <  1:17644/149‑1 (MQ=255)
                                                                                           ccAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAg                                                   >  5:176646/1‑148 (MQ=255)
                                                                                              aGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCaaa                                               <  3:293600/149‑1 (MQ=255)
                                                                                                  aGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGcc                                            >  7:282873/1‑148 (MQ=255)
                                                                                                     gCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCaaa                                         <  1:263983/148‑1 (MQ=255)
                                                                                                            gAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGAt                                 <  8:99917/149‑1 (MQ=255)
                                                                                                              aaCAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTg                               <  8:238570/149‑1 (MQ=255)
                                                                                                               aCAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGa                              <  8:282873/149‑1 (MQ=255)
                                                                                                                   gCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATa                            >  5:271022/1‑147 (MQ=255)
                                                                                                                     aCGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGTGATTGATACAg                         >  4:217291/1‑148 (MQ=255)
                                                                                                                       gCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGGTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCt                      <  6:160963/149‑1 (MQ=255)
                                                                                                                         tCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATa                            <  2:13340/141‑1 (MQ=255)
                                                                                                                         tCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATa                            >  1:13340/1‑141 (MQ=255)
                                                                                                                           aTTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGca                   >  3:110420/1‑148 (MQ=255)
                                                                                                                               gCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAg              <  7:1734/149‑1 (MQ=255)
                                                                                                                               gCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAg              >  7:4716/1‑149 (MQ=255)
                                                                                                                                 gAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAggg            <  3:223502/149‑1 (MQ=255)
                                                                                                                                     cAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTgc        <  3:31179/149‑1 (MQ=255)
                                                                                                                                     aattcggcgCTGAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTgc        >  8:224610/5‑149 (MQ=255)
                                                                                                                                      attcggcgCTGAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGAt                                 <  1:232217/120‑1 (MQ=255)
                                                                                                                                      attcggcgCTGAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGAt                                 >  2:232217/4‑123 (MQ=255)
                                                                                                                                      aGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTgcg       >  3:118989/1‑149 (MQ=255)
                                                                                                                                           gcgcTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAATTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTg  >  2:253856/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                               
GTCGCCGGGCGTTGATCGACCACCTTGGCGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTG  >  NC_002947/4707174‑4707461

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑