Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,963,068 G→A 6.3% A211V (GCC→GTC)  fleS ← two‑component system sensor histidine kinase FleS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,963,0680GA 6.3% 181.8 / 2.1 63A211V (GCC→GTC) fleStwo‑component system sensor histidine kinase FleS
Reads supporting (aligned to +/- strand):  ref base G (29/30);  new base A (2/2);  total (31/32)
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.62e-01

GGGGTCACCCGGTCCCCCAGTGGCAACTCGCCACGGGCAAACACCAGCATGTCGCGCACCTGGTGTTCCAGTTCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGCCCCAACGACGACAGCCGCTCGTGGCGAGCCAGCTGGTCCTGGAGGCGACGGGTTTCAGTCAGGTCGTTGAGCAGC  >  NC_002947/4962924‑4963217
                                                                                                                                                |                                                                                                                                                     
ggggTCACCCGGTCCCCCAGTGGCAACTCGCCACGGGCAAACACCAGCATGTCGCGCACCTGGTGTTCCAGTTCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAgg                                                                                                                                                  >  1:1095843/1‑150 (MQ=255)
    tCACCCGGTCCCCCATTGGCAACTCGCCACGGGCAAACACCAGCATGTCGCGCACCTGGTGTTCCAGTTCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGc                                                                                                                                              >  1:854158/1‑150 (MQ=255)
      aCCCGGTCCCCCAGTGGCAACTCGCCACGGGCAAACACCAGCATGTCGCGCACCTGGTGTTCCAGTTCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGCCGGCCACCTGGCCg                                                                                                                                            >  2:999111/1‑150 (MQ=255)
          ggTCCCCCAGTGGCAACTCGCCACGGGCAAACACCAGCATGTCGCGCACCTGGTGTTCCAGTTCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGCCGGCCACCTGGCCGACaa                                                                                                                                        <  1:999111/150‑2 (MQ=255)
           gTCCCCCAGTGGCAACTCGCCACGGGCAAACACCAGCATGTCGCGCACCTGGTGTTCCAGTTCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATa                                                                                                                                       <  2:478864/150‑1 (MQ=255)
              ccccAGTGGCAACTCGCCACGGGCAAACACCAGCATGTCGCGCACCTGGTGTTCCAGTTCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAg                                                                                                                                    >  1:845248/1‑150 (MQ=255)
                ccAGTGGCAACTCGCCACGGGCAAACACCAGCATGTCGCGCACCTGGTGTTCCAGTTCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCAAGGTGGCCGGCATACAGCa                                                                                                                                  >  2:68459/1‑150 (MQ=255)
                  aGTGGCAACTCGCCACGGGCAAACACCAGCATGTCGCGCACCTGGTGTTCCAGTTCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATg                                                                                                                                <  1:122886/150‑1 (MQ=255)
                         aCTCGCCACGGGCAAACACCAGCATGTCGCGCACCTGGTGTTCCAGTTCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCg                                                                                                                         >  1:927255/1‑150 (MQ=255)
                                   ggCAAACACCAGCATGTCGCGCACCTGGTGTTCCAGTTCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGt                                                                                                               <  2:58847/150‑1 (MQ=255)
                                         cacCAGCATGTCGCGCACCTGGTGTTCCAGTTCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGAt                                                                                                         <  2:384045/150‑1 (MQ=255)
                                            cAGCATGTCGCGCACCTGGTGTTCCAGTTCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCtg                                                                                                      >  1:192825/1‑150 (MQ=255)
                                               cATGTCGCGCACCTGGTGTTCCAGTTCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCa                                                                                                                                         <  2:786270/112‑1 (MQ=255)
                                               cATGTCGCGCACCTGGTGTTCCAGTTCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCa                                                                                                                                         >  1:786270/1‑112 (MQ=255)
                                                aTGTCGCGCACCTGGCGTTCCAGTTCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGAGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCtggtgg                                                                                                  <  2:228690/150‑1 (MQ=255)
                                                      cgcACCTGGTGTTCCAGTTCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGc                                                                                            <  2:663030/150‑1 (MQ=255)
                                                      cgcACCTGGTGTTCCAGTTCGTGCAGGCGCTCCCTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGc                                                                                            <  1:680864/150‑1 (MQ=255)
                                                         aCCTGGTGTTCCAGTTCGTGCAGGCGCTCCTTCAAGGTCGCGGCGACACAAGTACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAg                                                                                         <  1:979650/150‑1 (MQ=255)
                                                         aCCTGGTGTTCCAGTTCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAg                                                                                         >  2:261623/1‑150 (MQ=255)
                                                           cTGGTGTTCCAGTTCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATa                                                                                                                                       >  2:1147601/1‑102 (MQ=255)
                                                           cTGGTGTTCCAGTTCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATa                                                                                                                                       <  1:1147601/102‑1 (MQ=255)
                                                           cTGGTGTTCCAGTTCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGCCGGCCACCTGGCCGACaacga                                                                                                                                     >  2:478857/1‑100 (MQ=255)
                                                            cggTGTTCCAGTTCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGATTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGCGGCa                                                                                      >  2:328839/2‑150 (MQ=255)
                                                             ggTGTTCCAGTTCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCaa                                                                                     <  2:192825/150‑1 (MQ=255)
                                                                 ttCCAGTTCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCAt                                                                                 >  2:857125/1‑150 (MQ=255)
                                                                   ccAGTTCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATgc                                                                               <  1:205089/150‑1 (MQ=255)
                                                                     aGTTCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATgcgc                                                                             >  2:1100916/1‑150 (MQ=255)
                                                                      gTTCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGGCGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGcc                                                                            >  1:345628/1‑150 (MQ=255)
                                                                       ttCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGccc                                                                           >  1:432246/1‑150 (MQ=255)
                                                                       ttCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGATGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGccc                                                                           <  2:845248/150‑1 (MQ=255)
                                                                          gTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGCCCCaa                                                                        >  2:1044519/1‑150 (MQ=255)
                                                                             cAGGCGCTCCTTCAAGGCCGCGGCGAAGAGCTGTCGGGTTTGCGAGCTCAGCGCATGCTCGTTGTAGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGCCCCAacga                                                                     <  2:854158/150‑1 (MQ=255)
                                                                                gcgcTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGCCCCAAcgacga                                                                  >  2:658548/1‑150 (MQ=255)
                                                                                gcgcTCCTTCAAGGTCCTGTTGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGATCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTCCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGCCCCAAcgacga                                                                  <  1:367770/150‑1 (MQ=255)
                                                                                  gcTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGCCCCAACGACGACa                                                                <  1:857125/150‑1 (MQ=255)
                                                                                  gcTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGCCCCAACGACGACa                                                                >  1:256519/1‑150 (MQ=255)
                                                                                     ccTTCAGGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGACCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGCCCCAACGACGACa                                                                >  2:901747/1‑147 (MQ=255)
                                                                                     ccTTCAGGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGACCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGCCCCAACGACGACa                                                                <  1:901747/147‑1 (MQ=255)
                                                                                     ccTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGACCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGCCCCAACGACGACa                                                                >  2:901763/1‑147 (MQ=255)
                                                                                     ccTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGACCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGCCCCAACGACGACa                                                                <  1:901763/147‑1 (MQ=255)
                                                                                      cTNCAAGGTCGCGGCGAAGCGCTGACGNGTTTGCGCGCTCNGCGCTTNCNCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCGGCCGACAGCGGCGGGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGCCGCAACGCCGACAGCCg                                                            >  1:848727/1‑150 (MQ=255)
                                                                                          aaGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTgg                                                                                                                                               <  2:242993/63‑1 (MQ=255)
                                                                                          aaGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTgg                                                                                                                                               >  1:242993/1‑63 (MQ=255)
                                                                                            ggTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGCCCCAACGACGACAGCCGCTCGTg                                                      <  2:256519/150‑1 (MQ=255)
                                                                                            ggTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGCCCCAACGACGACAGCCGCTCGTg                                                      <  2:432246/150‑1 (MQ=255)
                                                                                             gTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGCCCCAACGACGACAGCCGCTCGTgg                                                     <  2:345628/150‑1 (MQ=255)
                                                                                                  ggcgAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGCCCCAACGACGACAGCCGCTCGTGGCGAgt                                                <  2:927255/150‑2 (MQ=255)
                                                                                                           cTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCtggt                                                                                                    <  1:735927/89‑1 (MQ=255)
                                                                                                           cTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCtggt                                                                                                    >  2:735927/1‑89 (MQ=255)
                                                                                                           cTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGCCCCAACGACGACAGCCGCTCGTGGCGAGCCAGCTGGTc                                       <  2:814858/150‑1 (MQ=255)
                                                                                                            tGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGCCCCAACGACGACAGCCGCTCGTGGCGCGCCAGGTGGTcc                                      >  2:445233/1‑150 (MQ=255)
                                                                                                              aCGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGCCCCAACGACGACAGCCGCTCGTGGCGAGCCAGCTGGTCCTg                                    >  1:1074730/1‑150 (MQ=255)
                                                                                                                gggTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGCCCCAACGACGACAGCCGCTCGTGGCGAGCCAGCTGGTCCTGGa                                  <  1:511859/150‑1 (MQ=255)
                                                                                                                gggTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGACGACAGCAGCGTGCGGATCTGGTGGGCCAGCCAGGCAACCATGCGCCCCAACGACGACAGCCGCTCGTGGCGAGCCAGCTGGTCCTGGa                                  <  2:546530/150‑1 (MQ=255)
                                                                                                                   tttGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGCCCCAACGACGACAGCCGCTCGTGGCGAGCCAGCTGGTCCTGGAGGc                               <  2:1074730/150‑1 (MQ=255)
                                                                                                                            tCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGCCCCAACGACGACAGCCGCTCGTGGCGAGCCAGCTGGTCCTGGAGGCGACGGGttt                      <  1:261623/150‑1 (MQ=255)
                                                                                                                                cgcTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGCGGCAACCACGGGCCCATGCCACAGCACTGTCTCGTTTCGGTCCAGCTGTTCCTGGAGGCGACGGGTTtcagt                  <  1:328839/150‑1 (MQ=255)
                                                                                                                                cgcTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGCCCCAACGACGACAGCCACTCGTGGCGAGCCAGCTGGTCCTGGAGGCGACGGGTTtcagt                  >  2:874043/1‑150 (MQ=255)
                                                                                                                                      cTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGGCAGCGAGGCAACCATGCGCCCCAACGACAACAGCCGCACGTGTCGAGCCAGCTGGTCCTGGAGGCGACGGGGTTCAGTGAGGTc            >  1:560107/1‑150 (MQ=255)
                                                                                                                                         tatgccGGCCACCTGGCCGACATACAGCATGGCTACCGACAGCGGCGTGCGGAACTGGTGGGCAAGCGAGGCAACCATGCGCCCCAACGACGACAGCCGCTCGTGGCGAGCCAGCTGGACCTGGAGGCGACGGGTTTCAGTCAGGTCGtt         >  2:194397/6‑150 (MQ=255)
                                                                                                                                          ttGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGCCCCAACGACGACAGCCGCTCGTGGCGAGCCAGCTGGTCCTGGAGGCGACGGGTTTCAGTCAGGTCGTTg        >  1:974213/1‑150 (MQ=255)
                                                                                                                                           tgccGGCCACCTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGCCCCAACGACGACAGCCGCTCGTTGCGAGCCATCTGGTCTTGTAGGCGACGGGTTTCAGTCAGGTCGTTGa       >  2:854872/4‑150 (MQ=255)
                                                                                                                                              cGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGCCCCAACGACGACAGCCGCTCGTGGCGAGCCAGCTGGTCCTGGAGGCGACGGGTTTCAGTCAGGTCGTTGagca    <  1:538911/150‑1 (MQ=255)
                                                                                                                                                gCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGCCCCAACGACGACAGCCGCTCGTGGCGAGCCAGCTGGTCCTGGAGGCGACGGGTTTCAGTCAGGTCGTTGagcagc  >  1:94371/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GGGGTCACCCGGTCCCCCAGTGGCAACTCGCCACGGGCAAACACCAGCATGTCGCGCACCTGGTGTTCCAGTTCGTGCAGGCGCTCCTTCAAGGTCGCGGCGAAGCGCTGACGGGTTTGCGCGCTCAGCGCTTGCTCGTTGTCGGCCAGGTGGCCGGCATACAGCATGGCTGCCGACAGCGGCGTGCGGATCTGGTGGGCCAGCGAGGCAACCATGCGCCCCAACGACGACAGCCGCTCGTGGCGAGCCAGCTGGTCCTGGAGGCGACGGGTTTCAGTCAGGTCGTTGAGCAGC  >  NC_002947/4962924‑4963217

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: ‑