Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,311,322 C→G 5.8% A71P (GCA→CCA)  PP_1144 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,311,3220CG 5.8% 254.9 / 3.4 80A71P (GCA→CCA) PP_1144membrane protein
Reads supporting (aligned to +/- strand):  ref base C (40/34);  new base G (3/2);  total (44/36)
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.80e-01

CTGCGTTCGAGGCTGTGCATGGCCAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGCAGAACCAAGAGCGGTTCATGACGGCA  >  NC_002947/1311177‑1311469
                                                                                                                                                 |                                                                                                                                                   
ctgcGTTCGAGGCTGTGCATGGCCAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAg                                                                                                                                                  >  4:348740/1‑149 (MQ=255)
 tgcGTTCGAGGCTGTGCATGGCCAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAgg                                                                                                                                                 >  8:100001/1‑149 (MQ=255)
     ttCGAGGCTGTGCATGGCCAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGgcagc                                                                                                                                             <  3:45548/149‑1 (MQ=255)
              gtgCATGGCCAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAg                                                                                                                                                  <  6:256429/135‑1 (MQ=255)
              gtgCATGGCCAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAg                                                                                                                                                  >  5:256429/1‑135 (MQ=255)
               tgCAGGGCCAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCaa                                                                                                                                    <  2:157637/148‑1 (MQ=255)
                 cATGGCCAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAg                                                                                                                                   >  4:247265/1‑147 (MQ=255)
                  aTGGCCAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCAt                                                                                                                                <  6:53468/149‑1 (MQ=255)
                     gCCAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGtt                                                                                                                             >  7:197481/1‑149 (MQ=255)
                      ccAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTc                                                                                                                            <  8:281565/149‑1 (MQ=255)
                      ccAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTc                                                                                                                            >  4:321327/1‑149 (MQ=255)
                      ccAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTc                                                                                                                            >  7:124321/1‑149 (MQ=255)
                           cAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTg                                                                                                                       >  8:124496/1‑149 (MQ=255)
                           cAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGCCCAGCATGCTGATGAAAAGCATGGCCAATATGCAGCCTGCCAGGCAGCAGGCTACAAGCATGTTCCATTg                                                                                                                       >  8:135836/1‑149 (MQ=255)
                               ccACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTcgg                                                                                                                   <  1:230692/149‑1 (MQ=255)
                                cACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGgcagca                                                                                                                                            <  4:266413/123‑1 (MQ=255)
                                cACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTcggt                                                                                                                  <  6:44055/149‑1 (MQ=255)
                                     ccAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGACGTAGTCCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTcggtgggcg                                                                                                             >  7:304028/1‑149 (MQ=255)
                                        gTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTcggtgggcggtg                                                                                                          <  7:135836/149‑1 (MQ=255)
                                        gTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTcggtgggcggtg                                                                                                          <  7:22877/149‑1 (MQ=255)
                                        gTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTcggtgggcggtg                                                                                                          <  8:124321/149‑1 (MQ=255)
                                              tgcgaaGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCAAAAAGGAAAGCATGATGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGgtcg                                                                                                    >  5:284192/6‑149 (MQ=255)
                                               gAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGACAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGgtcgt                                                                                                   >  1:151872/1‑149 (MQ=255)
                                                 gCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGtcgtcg                                                                                                 <  8:22370/149‑1 (MQ=255)
                                                 gCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGtcgtcg                                                                                                 >  8:81845/1‑149 (MQ=255)
                                                  cAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCa                                                                                                                          >  6:389106/1‑123 (MQ=255)
                                                  cAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCa                                                                                                                          <  5:389106/123‑1 (MQ=255)
                                                     tGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCATGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCtgt                                                                                             >  5:12999/1‑149 (MQ=255)
                                                       ataGAATGGTAATGGATTAGTCGTAGTGCACGTCGACGGGGTGCCGATAGGCTAGGATGCTGATGAAATGCATGGCCCATATGCCGCCTGGCAGGCTGCAGGCTCCTAGCATGTTCGATTGTCGGTGTGCGGTGGGGTCGTCGCTGTgg                                                                                           <  8:304028/147‑1 (MQ=255)
                                                       aaaGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGACAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCCCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGTGCGGTGGGGTCGTCGCTGTgg                                                                                           >  5:232183/1‑149 (MQ=255)
                                                       aaaGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTgg                                                                                           <  7:100001/149‑1 (MQ=255)
                                                                tCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCtgtg                                                                                            <  4:74385/139‑1 (MQ=255)
                                                                tCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCtgtg                                                                                            >  3:74385/1‑139 (MQ=255)
                                                                    ggATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCgg                                                                              <  2:339640/149‑1 (MQ=255)
                                                                       tACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCa                                                                           >  6:313714/1‑149 (MQ=255)
                                                                        aCGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCAt                                                                          >  2:167487/1‑149 (MQ=255)
                                                                          gtcgtAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGt                                                                        >  7:262429/1‑149 (MQ=255)
                                                                          gtcgtAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGCCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGt                                                                        >  2:75313/1‑149 (MQ=255)
                                                                             gtAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGa                                                                     <  5:313714/149‑1 (MQ=255)
                                                                                  gCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTg                                                                >  2:87/1‑149 (MQ=255)
                                                                                    aCGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTc                                                                       <  1:379617/140‑1 (MQ=255)
                                                                                    aCGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTc                                                                       >  2:379617/1‑140 (MQ=255)
                                                                                                    gAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCgcaa                                              <  8:197481/149‑1 (MQ=255)
                                                                                                    gAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTAACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCgcgc                                                >  8:111144/1‑147 (MQ=255)
                                                                                                     aaaGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCgcaag                                             >  3:171420/1‑149 (MQ=255)
                                                                                                       aGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGaaa                                                          <  1:98340/134‑1 (MQ=255)
                                                                                                       aGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGaaa                                                          >  2:98340/1‑134 (MQ=255)
                                                                                                           cAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCgcaa                                              >  5:11928/1‑142 (MQ=255)
                                                                                                           cAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCgcaa                                              <  6:11928/142‑1 (MQ=255)
                                                                                                           cAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACg                                        >  2:87369/1‑148 (MQ=255)
                                                                                                             gCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATCTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGaaa                                     >  1:213412/1‑149 (MQ=255)
                                                                                                              cATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGcc                                                    <  1:3101/133‑1 (MQ=255)
                                                                                                              cATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGcc                                                    >  2:3101/1‑133 (MQ=255)
                                                                                                              cATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAAt                                    <  1:167487/149‑1 (MQ=255)
                                                                                                               aTGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGATCGTGGCGAAACAGGCCCCGCAAGCAACCAAAt                                    >  4:81598/1‑148 (MQ=255)
                                                                                                                 gCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAgg                                 >  6:37955/1‑149 (MQ=255)
                                                                                                                  cTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGc                                <  1:87/149‑1 (MQ=255)
                                                                                                                   tgatgaAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGATCCATTGTCGGTGGGCGGTGGGGTCGTCGATGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGc                                >  2:361324/1‑148 (MQ=255)
                                                                                                                    gatgaAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGc                                >  1:157188/1‑147 (MQ=255)
                                                                                                                    gatgaAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTACAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGc                                >  4:117919/1‑147 (MQ=255)
                                                                                                                     atgaAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGa                              <  3:394222/148‑1 (MQ=255)
                                                                                                                     atgaAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGa                              <  8:262429/148‑1 (MQ=255)
                                                                                                                      tgaAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGc                            <  7:88716/149‑1 (MQ=255)
                                                                                                                         aaTGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTg                                                                >  4:345237/1‑110 (MQ=255)
                                                                                                                         aaTGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTg                                                                <  3:345237/110‑1 (MQ=255)
                                                                                                                         aaTGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGCAg                          >  7:237533/1‑148 (MQ=255)
                                                                                                                          aTGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGTGAAACAGGCCGCGCAAGCAACGa                                       >  8:87133/1‑134 (MQ=255)
                                                                                                                          aTGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGTGAAACAGGCCGCGCAAGCAACGa                                       <  7:87133/134‑1 (MQ=255)
                                                                                                                          aTGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGCAGaa                        >  7:86965/1‑149 (MQ=255)
                                                                                                                               tGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGCAGAACCAag                   <  3:117919/149‑1 (MQ=255)
                                                                                                                                   ccATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGCAGAACCAAGAGCg               <  4:171420/149‑1 (MQ=255)
                                                                                                                                       ctgctgcCTGGCAGGCGGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGCAGAACCAAGAGCGGtt            >  2:329054/2‑148 (MQ=255)
                                                                                                                                       ctgctgcCTGGCAGGCGGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGCAGAACCAAGAGCGGTTc           <  1:329054/148‑1 (MQ=255)
                                                                                                                                        tgctgcCTGGCAGGCGGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAATCAGGCCGCgcaa                                              >  2:144286/1‑113 (MQ=255)
                                                                                                                                        tgctgcCTGGCAGGCGGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAATCAGGCCGCgcaa                                              >  2:139480/1‑113 (MQ=255)
                                                                                                                                        tGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGCAGAACCAAGAGCGGTTc           <  1:87369/148‑1 (MQ=255)
                                                                                                                                        tGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGCAGAACCAAGAGCGGTTCa          <  2:213412/149‑1 (MQ=255)
                                                                                                                                           cgccTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGCAGAACCAAGAGCGGTTCATGa       >  2:356409/1‑149 (MQ=255)
                                                                                                                                            gccTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGCAGAACCAAGAGCGGTTCATGAc      >  3:164995/1‑149 (MQ=255)
                                                                                                                                             ccTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGGGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGCAGAACCAAGAGCGGTTCATGACg     <  3:348740/149‑1 (MQ=255)
                                                                                                                                                 ccAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGCAGAACCAAGAGCGGTTCATGACGGCa  <  8:86965/148‑1 (MQ=255)
                                                                                                                                                 ccAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGCAGAACCAAGAGCGGTTCATGACGGCa  <  3:81598/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CTGCGTTCGAGGCTGTGCATGGCCAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGCAGAACCAAGAGCGGTTCATGACGGCA  >  NC_002947/1311177‑1311469

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