Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,657,941 G→C 5.5% P185A (CCA→GCA)  PP_1450 ← TPS family activation/secretion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,657,9410GC 5.5% 236.8 / 3.0 72P185A (CCA→GCA) PP_1450TPS family activation/secretion protein
Reads supporting (aligned to +/- strand):  ref base G (32/36);  new base C (2/2);  total (34/38)
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.85e-01

GTGCCAGCGTGAAGCAACCTTGCCCGGCTGCACGACCACGCGGGTGCCGCCCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGTTCGATGGTTTCGACAAAGCCTTCGACAATGACGATGTCCAACGGCATGCCCTCTGCTGGGGGCTGGCGCAGGTAGGGGCGG  >  NC_002947/1657795‑1658083
                                                                                                                                                  |                                                                                                                                              
gTGCCAGCGTGAAGCAACCTTGCCCGGCTGCACGACCACGCGGGTGCCGCCCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGATCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCt                                                                                                                                              >  4:209783/1‑149 (MQ=255)
        gTGAAGCAACCTTGCCCGGCTGCACGACCACGCGGGTGCCGCCCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGc                                                                                                                                      <  7:103377/149‑1 (MQ=255)
          gAAGCAACCTTGCCCGGCTGCACGACCACGCGGGTGCCGCCCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGCCTTCCCAGGcct                                                                                                                                    >  4:47787/1‑147 (MQ=255)
           aaGCAACCTTGCCCGGCTGCACGACCACGCGGGTGCCGCCCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAgg                                                                                                                                   >  8:116477/1‑149 (MQ=255)
             gCAACCTTGCCCGGCTGCACGACCACGCGGGTGCCGCCCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGcc                                                                                                                                 <  4:112058/149‑1 (MQ=255)
                   ttGCCCGGCTGCACGACCACGCGGGTGCCGCCCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGCCACTGTACCTGcccgaccttg                                                                                                                                   <  2:251151/141‑8 (MQ=255)
                   ttGCCCGGCTGCACGACCACGCGGGTGCCGCCCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGCCACTGTACCTGcccgaccttg                                                                                                                                   >  1:251151/1‑134 (MQ=255)
                     gcaaggttGCACGACCACGCGGGTGCCGCCCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAgg                                                                                                                         <  3:44971/142‑1 (MQ=255)
                     gCCCGGCTGCACGACCACGCGGGTGCCGCCCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGaa                                                                                                                           <  7:116477/147‑1 (MQ=255)
                      cccGGCTGCACGACCACGCGGGTGCCGCCCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCGGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAgg                                                                                                                         <  7:122843/148‑1 (MQ=255)
                       ccGGCTGCACGACCACGCGGGTGCCGCCCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCa                                                                                                                       >  3:88333/1‑149 (MQ=255)
                          gCTGCACGACCACGCGGGTGCCGCCCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTg                                                                                                                    <  1:51710/149‑1 (MQ=255)
                          gCTGCACGACCACGCGGGTGCCGCCCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGATCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTg                                                                                                                    <  3:209783/149‑1 (MQ=255)
                           cTGCACGACCACGCGGGTGCCGCCCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGCCTTCCCAGGCCTGCTTGGGCAGCCACTgt                                                                                                                   <  3:47787/149‑2 (MQ=38)
                              cACGACCACGCGGGTGCCGCCCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGgaga                                                                                                                <  3:47235/149‑1 (MQ=255)
                              cACGACCACGCGGGTGCCGCCCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGga                                                                                                                  <  4:88333/147‑1 (MQ=255)
                              cACGACCACGCGGGTGCCGCCCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGga                                                                                                                  >  8:151977/1‑147 (MQ=255)
                               aCGACCACGCGGGTGCCGCCCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGga                                                                                                                  >  3:79350/1‑146 (MQ=255)
                                 gACCACGCGGGTGCCGCCCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGCGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGgaga                                                                                                                >  7:312699/1‑146 (MQ=255)
                                 gACCACGCGGGTGCCGCCCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGgagagag                                                                                                             <  4:80193/149‑1 (MQ=255)
                                     aCGCGGGTGCCGCCCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAg                                                                                                         <  4:333478/149‑1 (MQ=255)
                                     aCGCGGGTGCCGCCCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAg                                                                                                         >  5:233511/1‑149 (MQ=255)
                                     aCGCGGGTGCCGCCCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAg                                                                                                         <  6:233511/149‑1 (MQ=255)
                                      cgcgGGTGCCGCCCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAgg                                                                                                        <  7:151977/149‑1 (MQ=255)
                                       gcgGGTGCCGCCCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAggg                                                                                                       >  2:75619/1‑149 (MQ=255)
                                           gTGCCGCCCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAg                                                                                                   >  6:195415/1‑149 (MQ=255)
                                              ccgccCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCTGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTc                                                                                                >  7:185610/1‑149 (MQ=255)
                                               cgccCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGgaga                                                                                                                >  6:222748/1‑132 (MQ=255)
                                               cgccCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGgaga                                                                                                                >  6:4588/1‑132 (MQ=255)
                                               cgccCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGgaga                                                                                                                <  5:4588/132‑1 (MQ=255)
                                               cgccCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGgaga                                                                                                                <  5:222748/132‑1 (MQ=255)
                                               cgccCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCg                                                                                               >  4:138621/1‑149 (MQ=255)
                                                 cccTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCggg                                                                                             <  3:138621/149‑1 (MQ=255)
                                                   cTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGgc                                                                                           <  8:185610/149‑1 (MQ=255)
                                                    tGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAgg                                                                                                        <  2:237143/135‑1 (MQ=255)
                                                    tGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAgg                                                                                                        >  1:237143/1‑135 (MQ=255)
                                                            gCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGtt                                                                                  <  5:195415/149‑1 (MQ=255)
                                                                ggCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGTTCGAt                                                                              <  4:79350/149‑1 (MQ=255)
                                                                    gcaAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGgagagaga                                                                                                            >  3:275054/1‑115 (MQ=255)
                                                                    gcaAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGgagagaga                                                                                                            <  4:275054/115‑1 (MQ=255)
                                                                        ggTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGTTCGATGGTTTCGa                                                                      >  8:172148/1‑149 (MQ=255)
                                                                            gcgcCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGTTCGATGGTTTCGACaaa                                                                  <  1:75619/149‑1 (MQ=255)
                                                                            gcgcCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGTTCGATGGTTTCGACaaa                                                                  <  7:325985/149‑1 (MQ=255)
                                                                            gcgcCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGTTCGATGGTTTCGACa                                                                    >  6:150219/1‑147 (MQ=255)
                                                                              gcCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCaa                                                                                                                                        >  8:278988/1‑77 (MQ=255)
                                                                              gcCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCaa                                                                                                                                        <  7:278988/77‑1 (MQ=255)
                                                                               cccAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGTCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGTTCGATGGTTTCGACAAAGcc                                                               >  8:220655/1‑149 (MQ=255)
                                                                                ccAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGTTCGATGGTTTCGACAAAGCCt                                                              >  7:141152/1‑149 (MQ=255)
                                                                                  aGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGTTCGATGGTTTCGACAAAGCCTTc                                                            <  8:312699/149‑1 (MQ=255)
                                                                                      cATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGTCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGTTCGATGGTTTCGACAAAGCCTTCGa                                                          <  7:220655/147‑1 (MQ=255)
                                                                                       aTAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGgccagcca                                                                                     >  5:261471/1‑119 (MQ=255)
                                                                                       aTAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGgccagcca                                                                                     <  6:261471/119‑1 (MQ=255)
                                                                                        tAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGTTCGATGGTTTCGACAAAGCCTTCGACAAt                                                      <  5:150219/149‑1 (MQ=255)
                                                                                                     ggTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGTTCGATGGTTTCGACAAAGCCTTCGACAATGACGAt                                                >  4:238476/1‑142 (MQ=255)
                                                                                                     ggTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGTTCGATGGTTTCGACAAAGCCTTCGACAATGACGAt                                                <  3:238476/142‑1 (MQ=255)
                                                                                                     ggTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGTTCGATGGTTTCGACAAAGCCTTCGACAATGACGATGTCCAAc                                         >  3:173815/1‑149 (MQ=255)
                                                                                                           gCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGTCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCGAGTTCGATGGTTTCGACAAAGCCTTCGACAATGACGATGTCCAACGGc                                      <  4:142208/146‑1 (MQ=255)
                                                                                                           gCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGTTCGATGGTTTCGACAAAGCCTTCGACAATGACGAt                                                <  6:148054/136‑1 (MQ=255)
                                                                                                           gCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGTTCGATGGTTTCGACAAAGCCTTCGACAATGACGAt                                                >  5:148054/1‑136 (MQ=255)
                                                                                                           gCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGTTCGATGGTTTCGACAAAGCCTTCGACAATGACGATGTCCAACGGc                                      >  3:142208/1‑146 (MQ=255)
                                                                                                                     gCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGTTCGATGGTTTCGACAAAGCCTTCGACAATGACGATGTCCAACGGCATGCCCTctgctg                         <  7:172148/149‑1 (MQ=255)
                                                                                                                            tCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGTTCGATGGTTTCGACAAAGCCTTCGACAATGACGATGTCCAACGGCATGCCCTCTGCTgggg                      <  1:11877/145‑1 (MQ=255)
                                                                                                                            tCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGTTCGATGGTTTCGACAAAGCCTTCGACAATGACGATGTCCAACGGCATGCCCTCTGCTgggg                      >  2:11877/1‑145 (MQ=255)
                                                                                                                            tCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGTTCGATGGTTTCGACAAAGCCTTCGACAATGACGATGTCCAACGGCATGCCCTCTGCTgg                        >  7:306828/1‑143 (MQ=255)
                                                                                                                            tCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGTTCGATGGTTTCGACAAAGCCTTCGACAATGACGATGTCCAACGGCATGCCCTCTGCTgg                        <  8:306828/143‑1 (MQ=255)
                                                                                                                            tCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGTTCGATGGTTTCGACAAAGCCTTCGACAATGACGATGTCCAACGGCATGCCCTCTGCTGGGGGCTg                  <  8:141152/149‑1 (MQ=255)
                                                                                                                             cAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGTTCGATGGTTTCGACAAAGCCTTCGACAATGACGATGTCCAACGGCATGCCCTCTGCTGGGGGCTgg                 >  5:138842/1‑149 (MQ=255)
                                                                                                                                ggTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGTTCGATGGTTTCGAc                                                                     <  3:192521/94‑1 (MQ=255)
                                                                                                                                ggTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGTTCGATGGTTTCGAc                                                                     >  4:192521/1‑94 (MQ=255)
                                                                                                                                    gggCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGTTCGATGGTTTCGACAAAGCCTTCGACAATGACGATGTCCAACGGCATGCCCTCTGCTGGGGGCTGGCGCAgg           >  1:61104/1‑148 (MQ=255)
                                                                                                                                          gTACAGTGGCTGCCCAAGCAGGCCTGGGAACGCACTCGACAGACACAGGGGTAGCTCGGGGCCAGCCAGTACGATGGTGTCGACAAAGCCTTCGACAATGCCGATGTCAACCGGCATGCCCTCTGCTTGGGGCTGGCGCAGGTAGGggc    >  2:268767/1‑149 (MQ=255)
                                                                                                                                            aCAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGTTCGATGGTTTCGACAAAGCCTTCGACAATGACGATGTCCAACGGCATACCCTCTGCTGGGGGCTGGCGCAGGTAGGggcgg  >  1:3807/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
GTGCCAGCGTGAAGCAACCTTGCCCGGCTGCACGACCACGCGGGTGCCGCCCTGCAGTTCGCCGGGCAGCAAGGTGGCGCCCAGTTCATAGGCGCGAAGGCGGTTGAGCTGGTCCAGGCCTTGCTCAAGGTCGGGCAGGTACAGTGGCTGCCCAAGCAGGCCTGGGAAGGCACTGGAGAGAGACAGGGGTAGGTCGGGGCCAGCCAGTTCGATGGTTTCGACAAAGCCTTCGACAATGACGATGTCCAACGGCATGCCCTCTGCTGGGGGCTGGCGCAGGTAGGGGCGG  >  NC_002947/1657795‑1658083

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