Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,811,285 G→A 5.2% E205K (GAA→AAA)  priA → primosome assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,811,2850GA 5.2% 232.1 / 2.0 77E205K (GAA→AAA) priAprimosome assembly protein
Reads supporting (aligned to +/- strand):  ref base G (35/38);  new base A (2/2);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

CCCACAGCCTGCTGGCCAAACTCAACCTGAACAAGGACAGCCTCGACCTGCTGCTGGCCAAGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGTGGCAAGACCGAGGTCTACCTGCAACTGATCCGTGAAACCCTGGAAGCCGGCAAGCAGGCGCTGGT  >  NC_002947/5811139‑5811433
                                                                                                                                                   |                                                                                                                                                    
cccACAGCCTGCTGGCCAAACTCAACCTGAACAAGGACAGCCTCGACCTGCTGCTGGCCAAGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGaa                                                                                                                                                    <  5:116092/149‑1 (MQ=255)
cccACAGCCTGCTGGCCAAACTCAACCTGAACAAGGACAGCCTCGACCTGCTGCTGGCCAAGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGaa                                                                                                                                                    <  5:98137/149‑1 (MQ=255)
 ccACAGCCTGCTGGCCAAACTCAACCTGAACAAGGACAGCCTCGATCTGCTGCTGGCCAAGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAag                                                                                                                                                   <  5:160755/149‑1 (MQ=255)
 ccACAGCCTGCTGGCCAAACTCAACCTGAACAAGGACAGCCTCGACCTGCTGCTGGCCAAGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACg                                                                                                                                                      >  8:30760/1‑146 (MQ=255)
  cacaGCCTGCTGGCCAAACTCAACCTGAACAAGGACAGCCTCGACCTGCTGCTGGCCAAGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAaga                                                                                                                                                  <  8:194547/149‑1 (MQ=255)
   acaGCCTGCTGGCCAAACTCAACCTGAACAAGGACAGCCTCGACCTGCTGCTGGCCAAGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAagag                                                                                                                                                 >  5:206294/1‑149 (MQ=255)
          gctgGCCAAACTCAACCTGAACAAGGACAGCCTCGACCTGCTGCTGGCCAAGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAAcgc                                                                                                                                           >  8:134927/1‑148 (MQ=255)
              gCCAAACTCAACCTGAACAAGGACAGCCTCGACCTGCTGCTGGCCAAGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAg                                                                                                                                       <  2:142160/148‑1 (MQ=255)
              gCCAAACTCAACCTGAACAAGGACAGCCTCGACCTGCTGCTGGCCAAGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGc                                                                                                                                      >  1:86127/1‑149 (MQ=255)
               ccAAACTCAACCTGAACAAGGACAGCCTCGACCTGCTGCTGGCCAAGGACCTGGTACAGATCGAGAAGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGcc                                                                                                                                     >  3:133197/1‑149 (MQ=255)
                cAAACTCAACCTGAACAAGGACAGCCTCGACCTGCTGCTGGCCAAGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGaa                                                                                                                                        <  2:322364/145‑1 (MQ=255)
                cAAACTCAACCTGAACAAGGACAGCCTCGACCTGCTGCTGGCCAAGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAgag                                                                                                                                     >  1:322364/1‑146 (MQ=255)
                  aaCTCAACCTGAACAAGGACAGCCTCGACCTGCTGCTGGCCAAGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAagag                                                                                                                                                 <  5:33465/134‑1 (MQ=255)
                  aaCTCAACCTGAACAAGGACAGCCTCGACCTGCTGCTGGCCAAGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAagag                                                                                                                                                 >  6:33465/1‑134 (MQ=255)
                       aaCCTGAACAAGGACAGCCTCGACCTGCTGCTGGCCAAGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGc                                                                                                                             >  4:157622/1‑149 (MQ=255)
                       aaCCTGAACAAGGACAGCCTCGACCTGCTGCTGGCCAAGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGc                                                                                                                             >  5:183556/1‑149 (MQ=255)
                        aCCTGATCAAGGACAGCCTCGACCTGCTGCTGGCCAAGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGcc                                                                                                                            <  6:335258/149‑1 (MQ=255)
                         ccTGAACAAGGACAGCCTCGACCTGCTGCTGGCCAAGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCg                                                                                                                           >  6:343874/1‑149 (MQ=255)
                         ccTGAACAAGGACAGCCTCGACCTGCTGCTGGCCAAGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCg                                                                                                                           <  5:344627/149‑1 (MQ=255)
                          cTGAACAAGGACAGCCTCGACCTGCTGCTGGCCAAGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCgt                                                                                                                          >  7:323800/1‑149 (MQ=255)
                                  ggACAGCCTCGACCTGCTGCTGGCCAAGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAgg                                                                                                                  <  3:160524/149‑1 (MQ=255)
                                  ggACAGCCTCGACCTGCTGCTGGCCAAGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACAAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAgg                                                                                                                  >  6:333475/1‑149 (MQ=255)
                                  ggACAGCCTCGACCTGCTGCTGGCCAAGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACAAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAgg                                                                                                                  >  6:327553/1‑149 (MQ=255)
                                        ccTCGACCTGCTGCTGGCCAAGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTg                                                                                                            >  1:140689/1‑149 (MQ=255)
                                        ccTCGACCTGCTGCTGGCCAAGGACCTGGTACAGATCGAGACGCGACGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTg                                                                                                            >  6:360999/1‑149 (MQ=255)
                                            gACCTGCTGCTGGCCAAGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGc                                                                                                          <  3:239460/147‑1 (MQ=255)
                                            gACCTGCTGCTGGCCAAGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGc                                                                                                          <  3:269331/147‑1 (MQ=255)
                                                  ctgctgTCCAAGGACCTGGTACAGATCGAGACGCGCCTCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAAGTACCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCg                                                                                                   <  1:345342/148‑1 (MQ=255)
                                                          cAAGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCt                                                                                          >  8:162260/1‑149 (MQ=255)
                                                           aaGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCtt                                                                                         <  7:207576/149‑1 (MQ=255)
                                                           aaGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCtt                                                                                         >  5:175748/1‑149 (MQ=255)
                                                            aGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTg                                                                                        <  8:323800/149‑1 (MQ=255)
                                                             ggACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGa                                                                                                                                         <  4:243334/99‑1 (MQ=255)
                                                             ggACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGa                                                                                                                                         >  3:243334/1‑99 (MQ=255)
                                                             ggACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCgg                                                                                                  <  8:163538/138‑1 (MQ=255)
                                                             ggACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCgg                                                                                                  >  7:163538/1‑138 (MQ=255)
                                                                ccTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTgg                                                                                    <  7:162260/149‑1 (MQ=255)
                                                                     tACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGc                                                                               <  2:140689/149‑1 (MQ=255)
                                                                     tACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGc                                                                               <  3:331286/149‑1 (MQ=255)
                                                                      aCAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCg                                                                              >  7:193101/1‑149 (MQ=255)
                                                                       cAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCg                                                                              <  8:193101/148‑1 (MQ=255)
                                                                         gATCGAGACGCGCCGCCACTTGCCGGGACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCa                                                                           >  1:203270/1‑149 (MQ=255)
                                                                          aTCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCAc                                                                          >  5:6063/1‑149 (MQ=255)
                                                                           tCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCAcc                                                                         <  8:103214/149‑1 (MQ=255)
                                                                            cGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCg                                                                        >  6:330626/1‑149 (MQ=255)
                                                                             gagaCGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCgg                                                                       >  1:124295/1‑149 (MQ=255)
                                                                              agaCGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCgg                                                                       >  2:317060/1‑148 (MQ=255)
                                                                                aCGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAg                                                                    >  2:123957/1‑149 (MQ=255)
                                                                                aCGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAg                                                                    >  5:14032/1‑149 (MQ=255)
                                                                                 cgcgCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGt                                                                   <  5:343874/149‑1 (MQ=255)
                                                                                    gccgccACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGTgg                                                                 >  2:49926/1‑148 (MQ=255)
                                                                                        ccACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGTGGCAAGa                                                            >  7:284929/1‑149 (MQ=255)
                                                                                              gCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGTGGCAAGACCGAg                                                       <  1:123957/148‑1 (MQ=255)
                                                                                              gCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACAAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGTGGCAAGACCGAg                                                       <  5:327553/148‑1 (MQ=255)
                                                                                              gCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACAAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGTGGCAAGACCGAg                                                       <  5:333475/148‑1 (MQ=255)
                                                                                                  gCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGTGGCAAGACCGAGGTCTa                                                  >  2:168453/1‑149 (MQ=255)
                                                                                                   cACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGACACCGGCAGTGGCAAGACCGAGGTCTAc                                                 <  1:211762/149‑1 (MQ=255)
                                                                                                      tACGTCATGAGCACTGGCTGGCACAGCCGGAACTGGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTg                                                                                        >  1:247152/1‑108 (MQ=255)
                                                                                                      tACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTg                                                                                        <  2:247152/107‑1 (MQ=255)
                                                                                                        cGTCATGAGCACTGGCTGGCACAGCCGTAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGTGGCAAGACCGAGGTCTACCTGc                                             <  2:292856/148‑1 (MQ=255)
                                                                                                         gTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGTCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGTGGCAAGACCGAGGTCTACCTGCaa                                           <  1:168453/149‑1 (MQ=255)
                                                                                                         gTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGTGGCAAGACCGAGGTCTACCTGCaa                                           <  2:86127/149‑1 (MQ=255)
                                                                                                                gCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGTGGCAAGACCGAGGTCTACCTGCAACTGATcc                                    <  7:134927/149‑1 (MQ=255)
                                                                                                                gCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGTGGCAAGACCGAGGTCTACCTGCAACTGATcc                                    <  2:124295/149‑1 (MQ=255)
                                                                                                                  aCTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCAGTGTGCGAGGGTTTTTGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGTGGCAAGACCGAGGTCTACCTGCAACTGATCCGt                                  <  4:133197/149‑1 (MQ=255)
                                                                                                                      gctggcACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGTGGCAAGACCGAGGTCTACCTGCAACTGATCCGTGaaa                              >  6:357867/1‑149 (MQ=255)
                                                                                                                          gcACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGTGGCAAGACCGAGGTCTACCTGCATCTGATCCGTGAAACCCt                          <  6:6063/149‑1 (MQ=255)
                                                                                                                           cacaGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGTGGCAAGACCGAGGTCTACCTGCAACTGATCCGTGAAACCCTg                         >  4:370383/1‑149 (MQ=255)
                                                                                                                            acacCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGTGGCAAGACCGAGGTCTACCTGCAACTGATCCGTGAAACCCTgg                        <  6:356250/149‑1 (MQ=255)
                                                                                                                            acaGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGTGGCAAGACCGAGGTCTACCTGCAACTGATCCGTGAAACCCTg                         >  2:211242/1‑148 (MQ=255)
                                                                                                                              aGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGTGGCAAGACCGAGGTCTACCTGCAACTGATCCGTGAAACCCTGGaa                      <  7:30760/149‑1 (MQ=255)
                                                                                                                                   gAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGTGGCAAGACCGAGGTCTACCTGCAACTGATCCGTGAAACCCTGGAAGCCgg                 <  6:183556/149‑1 (MQ=255)
                                                                                                                                   gAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGTGGCAAGACCGAGGTCTACCTGCAACTGATCCGTGAAACCCTGGAAGCCgg                 >  4:104694/1‑149 (MQ=255)
                                                                                                                                         gtcGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGTGGCAAGACCGAGGTCTACCTGCAACTGATCCGTGAAACCCTGGAAGCCGCCAAg             >  4:253352/3‑148 (MQ=255)
                                                                                                                                          ccGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGTGGCAAGACCGAGGTCTACCTGCAACTGATCCGTGAAACCCTGGAAGCCGGCAAg             >  6:327809/1‑147 (MQ=255)
                                                                                                                                           cGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGTGGCAAGACCGAGGTCTACCTGCAACTGAGCCGTGAAACCCTGGAAGCCGGCAAGCAg          <  5:357867/149‑1 (MQ=255)
                                                                                                                                               cAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGTGGCAAGACCGAGGTCTACCTGCAACTGATCCGTGAAACCCTGGAAGCCGGCAAGCAGgcgc      >  6:26995/1‑149 (MQ=255)
                                                                                                                                                  cGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGTGGCAAGACCGAGGTCTACCTGCAACTGATCCGTGAAACCCTGGAAGCCGGCAAGCAGGCGCTgg   <  3:157622/149‑1 (MQ=255)
                                                                                                                                                   gAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGTGGCAAGACCGAGGTCTACCTGCAACTGATCCGTGAAACCCTGGAAGCCGGCAAGCAGGCGCTGGt  >  6:240365/1‑149 (MQ=255)
                                                                                                                                                   gAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGTGGCAAGACCGAGGTCTACCTGCAACAGATCCGTCAAACCCTGGAAGCCGGCAACCAGGCGCTGGt  >  4:100028/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCCACAGCCTGCTGGCCAAACTCAACCTGAACAAGGACAGCCTCGACCTGCTGCTGGCCAAGGACCTGGTACAGATCGAGACGCGCCGCCACTTGCCGGCACTACGTCATGAGCACTGGCTGGCACAGCCGGAACTGCCGCTCAACGAAGAGCAACGCGAAGCCTTCGACGCCGTGTGCGAGGGTTTTGGCGGTTTCGGCGCGTTCTTGCTGGCCGGCGTCACCGGCAGTGGCAAGACCGAGGTCTACCTGCAACTGATCCGTGAAACCCTGGAAGCCGGCAAGCAGGCGCTGGT  >  NC_002947/5811139‑5811433

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: