Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 479,502 C→T 7.3% intergenic (‑34/+22) cca ← / ← ycgB multifunctional tRNA nucleotidyltransferase/2',3'‑cyclic phosphodiesterase/2' nucleotidase/phosphatase/type IV piliation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947479,5020CT 7.3% 244.4 / 5.7 82intergenic (‑34/+22)cca/ycgBmultifunctional tRNA nucleotidyltransferase/2',3'‑cyclic phosphodiesterase/2' nucleotidase/phosphatase/type IV piliation protein
Reads supporting (aligned to +/- strand):  ref base C (39/37);  new base T (4/2);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.63e-01

CCATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGTCTCCAGGTGGATATC  >  NC_002947/479356‑479649
                                                                                                                                                   |                                                                                                                                                    
ccATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGt                                                                                                                                                    >  3:218053/1‑149 (MQ=255)
ccATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGt                                                                                                                                                    >  2:290564/1‑149 (MQ=255)
       cTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACGGGTTATCCTGTcggcag                                                                                                                                              <  5:213215/148‑3 (MQ=255)
       cTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGt                                                                                                                                             >  8:96697/1‑149 (MQ=255)
          gCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCAt                                                                                                                                          >  3:130944/1‑149 (MQ=255)
             tttCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATggg                                                                                                                                       <  4:139935/149‑1 (MQ=255)
                   cGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTATCCATAACTGCCCACAGGATAACCCGTCGGCGGTCATGGGGTTGCa                                                                                                                                 >  1:256971/1‑149 (MQ=255)
                       ggTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTGCATAACTGCCGACAGGATAACCCGTCGGGGGTCATGGGGTTGCATCag                                                                                                                             >  3:369035/1‑149 (MQ=255)
                         tGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCagag                                                                                                                           <  4:354209/149‑1 (MQ=255)
                              cccACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAggtggt                                                                                                                      <  6:98023/149‑1 (MQ=255)
                              cccACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAggtggt                                                                                                                      <  6:62481/149‑1 (MQ=255)
                                  ccaccaGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGAt                                                                                                                  <  4:351678/149‑1 (MQ=255)
                                   caccaGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGAt                                                                                                                  >  4:211927/1‑148 (MQ=255)
                                      cagccagTCGATATCCCGGACAGGGCGGCCGACCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACa                                                                                                              <  4:369035/149‑1 (MQ=255)
                                        gccagTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGc                                                                                                            >  6:341658/1‑149 (MQ=255)
                                           agTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAg                                                                                                         <  2:241755/149‑1 (MQ=255)
                                            gTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAgg                                                                                                        >  3:184682/1‑149 (MQ=255)
                                               gATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGTCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTcc                                                                                                     <  5:234710/149‑1 (MQ=255)
                                                atatCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCa                                                                                                    >  7:335530/1‑149 (MQ=255)
                                                 tatCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCa                                                                                                    <  8:218270/148‑1 (MQ=255)
                                                   tCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGc                                                                                                 <  1:172934/149‑1 (MQ=255)
                                                     gCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCgg                                                                                               <  4:184682/149‑1 (MQ=255)
                                                       tGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGGATGTGCATGTGCAGCTCTCCATAAATGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGc                                                                                              >  1:295102/1‑149 (MQ=255)
                                                       tGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGAATAACAGCCAGGTCCATGCGGc                                                                                              >  3:331294/1‑148 (MQ=255)
                                                         aCAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGcc                                                                                             <  7:96697/147‑1 (MQ=255)
                                                         aCAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCa                                                                                            >  8:283797/1‑148 (MQ=255)
                                                           aGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAg                                                                                         >  2:213913/1‑149 (MQ=255)
                                                              gcggcCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCg                                                                                      <  4:218053/149‑1 (MQ=255)
                                                              gcggcCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCg                                                                                      >  8:232613/1‑149 (MQ=255)
                                                                  ccGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGgcgc                                                                                  >  8:154090/1‑149 (MQ=255)
                                                                  ccGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGgcgc                                                                                  >  3:164857/1‑149 (MQ=255)
                                                                  ccGAGCAAACGGTCGCGCACGGCGCCGCCAACTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCt                                                                                 >  1:358384/1‑149 (MQ=255)
                                                                     aGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCttt                                                                               <  4:130944/149‑1 (MQ=255)
                                                                      gCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTc                                                                              >  5:172694/1‑149 (MQ=255)
                                                                        aaaCGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTGGGCGGTCATGGGGTTGCATCAGAAGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGAGCTTTCGc                                                                            >  1:128162/1‑149 (MQ=255)
                                                                        aaaCGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGc                                                                            <  5:49276/149‑1 (MQ=255)
                                                                          aCGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTa                                                                          <  7:154090/149‑1 (MQ=255)
                                                                           cGGTCGCGCACGGCGCCGCCAACTTTGTAGCTGTGCATGTGCAGCTCTCCATAACTGCCGACCGGATATCCCGGCGGCGGTCATGGGGTTGCCTCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTAt                                                                         <  8:250371/149‑1 (MQ=255)
                                                                           cGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTAt                                                                         <  3:360449/149‑1 (MQ=255)
                                                                           cGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTAt                                                                         >  2:377627/1‑149 (MQ=255)
                                                                                gcgcACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGCTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTcg                                                                    >  4:71852/1‑149 (MQ=255)
                                                                                        cgccgcCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCgcgg                                                            <  2:128162/149‑1 (MQ=255)
                                                                                        cgccgcCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCgcgg                                                            <  1:213913/149‑1 (MQ=255)
                                                                                            gccAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGcggcgg                                                        >  6:64421/1‑148 (MQ=255)
                                                                                              cAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCgcggcg                                                          <  7:56769/145‑1 (MQ=255)
                                                                                              cAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCgcggcg                                                          >  8:56769/1‑145 (MQ=255)
                                                                                              cAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCAt                                                     >  5:291017/1‑149 (MQ=255)
                                                                                               aaCTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCTGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCAtg                                                    >  2:367114/1‑149 (MQ=255)
                                                                                               aaCTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCAtg                                                    >  6:282549/1‑149 (MQ=255)
                                                                                                    tGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAg                                                                                                         >  8:131211/1‑92 (MQ=255)
                                                                                                    tGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAg                                                                                                         <  7:131211/92‑1 (MQ=255)
                                                                                                    tGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATgtggt                                                <  5:282549/148‑1 (MQ=255)
                                                                                                    tGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATgtggt                                                <  5:64421/148‑1 (MQ=255)
                                                                                                    tGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATgtgg                                                 <  1:351955/147‑1 (MQ=255)
                                                                                                    tGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATgtgg                                                 <  1:367114/147‑1 (MQ=255)
                                                                                                    tGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGAATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATgtggt                                                <  4:331294/148‑1 (MQ=255)
                                                                                                    tGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTAATGGGGTTGCATCAGAGGTGGTGTAAAACAGCCAGGTCCATGCGGCCATAGTCTGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATgtgg                                                 <  2:369067/147‑1 (MQ=255)
                                                                                                     gTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCaga                                                                                                                            <  5:209848/72‑1 (MQ=255)
                                                                                                     gTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCaga                                                                                                                            >  6:209848/1‑72 (MQ=255)
                                                                                                     gTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATgtggtg                                               <  7:283797/148‑1 (MQ=255)
                                                                                                      tAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAg                                             <  2:358384/149‑1 (MQ=255)
                                                                                                             tgcatgtgcaGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACtca                                       >  2:41253/1‑149 (MQ=255)
                                                                                                             tgcatgtgcaGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACtca                                       <  1:290564/149‑1 (MQ=255)
                                                                                                               catgtgcaGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGAGCAAGTGATGAGTACtcac                                      <  5:341658/149‑2 (MQ=255)
                                                                                                               catgtacaGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGACATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACtcat                                      >  1:80023/1‑148 (MQ=255)
                                                                                                                          ccccATGACCGCCGACGGGTTATCCTGTCGGCAGTTATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTa                                           <  1:139343/130‑1 (MQ=37)
                                                                                                                          ccccATGACCGCCGACGGGTTATCCTGTCGGCAGTTATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTa                                           >  2:139343/3‑132 (MQ=37)
                                                                                                                            ccATGACCGCCGACGGGTTATCCTGTCGGCAGTTATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCtt                        >  1:386712/1‑149 (MQ=38)
                                                                                                                            ccATGACCGCCGACGGGTTATCCTGTCGGCAGTTATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCtt                        >  6:239415/1‑149 (MQ=38)
                                                                                                                            ccATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCtt                        >  1:65839/1‑149 (MQ=255)
                                                                                                                              aTAACTGCCGACAGGATAACCCGTCGCCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCAAGGACCATGCGGCCATAGTCGGCGCTTTCGCTATCCTCGCCGAGCGCCGGCATGTGGTGAGTACTCAACACTTGCTCGACTTg                       >  5:235143/1‑148 (MQ=255)
                                                                                                                                aaCTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCAc                    >  2:229171/1‑149 (MQ=255)
                                                                                                                                 aCTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACg                   <  6:172694/149‑1 (MQ=255)
                                                                                                                                   tGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGAGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGt                 >  1:348605/1‑149 (MQ=255)
                                                                                                                                    gCCGACGGGTTATCCTGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGtc                >  6:259863/1‑149 (MQ=255)
                                                                                                                                     ccGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGtct               >  6:11042/1‑149 (MQ=255)
                                                                                                                                     ccGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGtc                <  6:291017/148‑1 (MQ=255)
                                                                                                                                     ccGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGtc                <  2:295102/148‑1 (MQ=255)
                                                                                                                                        aCAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGtctc              >  7:322521/1‑147 (MQ=255)
                                                                                                                                           ggATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGTCTCCAGGt         <  7:232613/149‑1 (MQ=255)
                                                                                                                                            gATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGTCTCCAGGTg        >  2:362776/1‑149 (MQ=255)
                                                                                                                                            gATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGTCTCCAGGTg        >  4:1986/1‑149 (MQ=255)
                                                                                                                                             ttatcccGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGTCTCCAGGTg        >  8:349957/5‑148 (MQ=255)
                                                                                                                                                 cccGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGTCTCCAGGTGGatat   <  4:164857/149‑1 (MQ=255)
                                                                                                                                                 cccGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGTCTCCAGGTGGatat   <  1:377627/149‑1 (MQ=255)
                                                                                                                                                  ccGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGTCTCCAGGTGGATATc  >  6:381907/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGTCTCCAGGTGGATATC  >  NC_002947/479356‑479649

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: