Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,950,411 G→A 19.7% intergenic (+777/+21) PP_5673 → / ← flhB D‑alanine‑‑D‑alanine ligase/flagellin export apparatus substrate specificity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,950,4110GA19.7% 135.8 / 28.4 67intergenic (+777/+21)PP_5673/flhBD‑alanine‑‑D‑alanine ligase/flagellin export apparatus substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base G (45/8);  new base A (7/6);  total (53/14)
Fisher's exact test for biased strand distribution p-value = 2.35e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.12e-02

TACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGC  >  NC_002947/4950298‑4950546
                                                                                                                 |                                                                                                                                       
tACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGATACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCCCCGGCGGCCCTCGATCTCTAAAGGAACGTCAACCTAACCCCCAACAGAta                                                                                                                   >  1:454231/1‑135 (MQ=18)
   aaacaaGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCGCGATCTCACAGGCAACGCAAAACTATAGCCGAACACATCAAc                                                                                                                >  1:632930/1‑136 (MQ=255)
           cGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCgcgg                                                                                                        >  2:593525/1‑136 (MQ=255)
             gTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCgcggag                                                                                                      >  2:345752/1‑134 (MQ=255)
             gTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCGCGATCTCACAGGCAACGCAAACCTATAGGCCAACAGATAAACAGTCgcggcc                                                                                                      >  1:569599/1‑135 (MQ=21)
                 cgcgcgcCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAgg                                                                                                  >  1:639522/1‑136 (MQ=255)
                                 gaataaCTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCa                                                                                  >  1:343028/6‑136 (MQ=255)
                                   ataaCTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAgg                                                                                >  1:621134/4‑136 (MQ=255)
                                         ggCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGAACAACTGTCGCGGCGCAGGTCCGGCCGGATCGGCAGGtcctcc                                                                          >  1:398758/1‑136 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGCTCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                                     >  1:616730/1‑136 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGACAAAAACCTCACGGCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                                     >  1:109863/1‑136 (MQ=18)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGACAAAAACCTATCGCCGCACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGAACGGCAGGGAGTCCTTCaa                                                                     >  2:553320/1‑135 (MQ=18)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGACAAAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCggg                                                                                         >  2:502585/1‑116 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGACAAAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCggg                                                                                         <  1:502585/116‑1 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGACAAAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                                     >  2:635221/1‑136 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGACAAAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                                     >  2:446298/1‑136 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGACAAAAACCTATCGCCGAACAGATCAACTGACGCGGCGCAGGACCGGCGGGATCGGCAGGTCGTCCTTCaa                                                                     >  2:77382/1‑135 (MQ=21)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCCACGCAAACCAAACGCCGAACACATCAACAGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                                     >  2:394586/1‑136 (MQ=25)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAAc                                                                                                                        <  2:675255/85‑1 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAAc                                                                                                                        >  1:675255/1‑85 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCaa                                                                     >  1:409342/1‑135 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                                     >  1:219757/1‑136 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                                     >  1:142390/1‑136 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                                     >  2:143322/1‑136 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                                     >  2:7251/1‑136 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                                     >  2:255697/1‑136 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                                     >  1:216104/1‑136 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                                     >  2:577590/1‑136 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                                     >  2:673363/1‑136 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                                     >  2:6359/1‑136 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                                     >  1:298004/1‑136 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                                     >  1:58784/1‑136 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                                     >  2:630410/1‑136 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                                     >  2:32430/1‑136 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                                     >  2:622493/1‑136 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATAAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                                     >  1:493650/1‑136 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCAATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                                     >  2:435514/1‑136 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACCCAAACCAATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                                     >  1:660604/1‑136 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATATCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                                     >  2:586884/1‑136 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGAACTCACAGGCAACGCAAAACTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                                     >  1:657952/1‑136 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCGCGATCTCCCAGGCAACGCAAACCAATCGCCGCAAACAACAACTGTCCCCGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTAAg                                                                     >  2:515786/1‑136 (MQ=21)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCGCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGCTCGGCAGGGCGTCCGTCAg                                                                     >  2:25101/1‑136 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCGCGATCTCACAGGCAACGCAAACCAATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                                     >  2:578399/1‑136 (MQ=255)
                                              gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCCCGCGGGCGGCGCTCGAAATCACACGCAACGCCAACCGATCGCCGAACAGATAAAATGTCGCCGCGCCGGGCCGGGGGGATCGGCAGGGCGTCCttaac                                                                     >  1:95955/1‑133 (MQ=21)
                                              gAAACAAATGTAGGAGCGAGCTAGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCAAACAGGCCACGCAAACCTATAGCCGAACAGAAGAACAGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTAGTCCTTCAg                                                                     >  2:476515/1‑136 (MQ=25)
                                               aaaCAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGt                                                                    >  1:119023/1‑136 (MQ=255)
                                                  cAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGc                                                                 >  2:260461/1‑136 (MQ=255)
                                                   aaaTGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGct                                                                <  2:409342/136‑1 (MQ=255)
                                                        tAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGACAAAAACCTCACGGCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTggg                                                           <  2:109863/136‑1 (MQ=21)
                                                             gcgagcTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGACAAAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCActt                                                      >  1:26959/1‑136 (MQ=255)
                                                                agcTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGcttgcc                                                   >  2:247347/1‑136 (MQ=255)
                                                                agcTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGcttgcc                                                   >  2:610207/1‑136 (MQ=255)
                                                                   ttGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGCTCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGcttgccctt                                                >  2:336276/1‑136 (MQ=255)
                                                                   ttGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGcttgccctt                                                >  2:332890/1‑136 (MQ=255)
                                                                   ttGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGcttgccctt                                                >  1:510859/1‑136 (MQ=255)
                                                                   ttGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGcttgccctt                                                >  1:648878/1‑136 (MQ=255)
                                                                     gCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGACAAAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCActtgcccttgc                                              <  2:26959/136‑1 (MQ=255)
                                                                     gCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGcttgcccttgc                                              <  1:622493/136‑1 (MQ=255)
                                                                     gCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGcttgcccttgc                                              <  1:673363/136‑1 (MQ=255)
                                                                     gCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGcttgcccttgc                                              <  1:213998/136‑1 (MQ=255)
                                                                             tGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCgg                                      >  1:196645/1‑136 (MQ=255)
                                                                                      cgggcggCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGa                             <  1:260461/136‑1 (MQ=255)
                                                                                                tCGATCTCACAGGCGACAAAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAaca                   <  1:553320/136‑1 (MQ=255)
                                                                                                tCGATCTCACAGGCGACAAAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAaca                   <  1:77382/136‑1 (MQ=255)
                                                                                                tCGATCTCACAGGCGACAAAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAaca                   <  1:446298/136‑1 (MQ=255)
                                                                                                tCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAaca                   <  1:32430/136‑1 (MQ=255)
                                                                                                tCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAaca                   <  1:247347/136‑1 (MQ=255)
                                                                                                                 gCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTgc  >  2:352956/1‑136 (MQ=255)
                                                                                                                 |                                                                                                                                       
TACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGC  >  NC_002947/4950298‑4950546

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