Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 886,083 G→C 5.3% W665S (TGG→TCG) ‡ PP_0767 → LuxR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947886,0830GC 5.3% 245.4 / 3.9 75W665S (TGG→TCG) ‡PP_0767LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (37/34);  new base C (2/2);  total (39/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 = 8.05e-01

TGCGCGCGGGGCTGGGCGAGGCGCGAGCGTGCCGGGACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGAGCAGCCGGCGGCGGC  >  NC_002947/885936‑886224
                                                                                                                                                   |                                                                                                                                             
tGCGCGCGGGGGTGGGCGAGGCGCGTGGGTGCCGGGACATCAGTGTGCTCTTCGGTGATTGCGTGTTTGCCACGCTCGATGGCCGGGATGGTCACTTTGCCGAGGCCTTTGCGGTGCTGGCCGAGGCCGAGCGCTTGATGCACATCTg                                                                                                                                               <  7:60124/148‑1 (MQ=255)
 gcgcgcGGGGCTGGGCGAGGCGCGAGCGTGCCGGGACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTgg                                                                                                                                              <  3:283926/148‑1 (MQ=255)
            tGGGCGAGGCGCGAGCGTGCCGGGACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCg                                                                                                                                  <  5:142674/149‑1 (MQ=255)
              ggCGAGGCGCGAGCGTGTCGGGACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGt                                                                                                                                >  5:269594/1‑149 (MQ=255)
              ggCGAGGCGCGAGCGTGCCGGGACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGt                                                                                                                                <  7:146100/149‑1 (MQ=255)
                cGAGGCGCGAGCGTGCCGGGACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGAAGCACATCTGGGATGTGCCCCCg                                                                                                                                  >  6:409660/1‑145 (MQ=255)
                cGAGGCGCGAGCGTGCCGGGACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGAAGCACATCTGGGATGTGCCCCCg                                                                                                                                  <  5:409660/145‑1 (MQ=255)
                  aGGCGCGAGCGTGTCGGGACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTctac                                                                                                                            <  6:269594/149‑1 (MQ=255)
                      gcgAGCGTGCCGGGACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCCCAGATGTGcatcaagc                                                                                                                                <  8:171439/141‑8 (MQ=255)
                      gcgAGCGTGCCGGGACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCCCAGATGTGcatcaagc                                                                                                                                >  7:171439/1‑134 (MQ=255)
                        gAGCGTGCCGGGACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTg                                                                                                                      >  7:290966/1‑149 (MQ=255)
                         aGCGTGCCGGGACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTgg                                                                                                                     <  2:67555/149‑1 (MQ=255)
                                 gggACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGAt                                                                                                              >  8:362649/1‑148 (MQ=255)
                                  ggACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTa                                                                                                            <  3:205335/149‑1 (MQ=255)
                                   gACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTgg                                                                                                                     >  1:76531/1‑139 (MQ=255)
                                   gACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTgg                                                                                                                     <  2:76531/139‑1 (MQ=255)
                                    aCATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTAcc                                                                                                          >  8:19499/1‑149 (MQ=255)
                                      aTCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCt                                                                                                        >  5:374000/1‑149 (MQ=255)
                                         aGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTgg                                                                                                      >  7:97413/1‑148 (MQ=255)
                                          gtgtGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGt                                                                                                     >  5:279705/1‑148 (MQ=255)
                                           tgtgCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTc                                                                                                    >  7:36651/1‑148 (MQ=255)
                                             tgCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCaaa                                                                                                 >  8:60701/1‑149 (MQ=255)
                                              gCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGAtgt                                                                                                                                         >  2:269473/1‑108 (MQ=255)
                                              gCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGAtgt                                                                                                                                         <  1:269473/108‑1 (MQ=255)
                                               cTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAAt                                                                                                <  8:290966/148‑1 (MQ=255)
                                               cTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATg                                                                                               >  4:389557/1‑149 (MQ=255)
                                                tCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGc                                                                                              >  6:28069/1‑149 (MQ=255)
                                                  aTCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGc                                                                                              <  7:362649/147‑1 (MQ=255)
                                                  aTCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCg                                                                                             >  8:7827/1‑148 (MQ=255)
                                                   tCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGAt                                                                                                                                           >  8:50934/1‑101 (MQ=255)
                                                   tCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGAt                                                                                                                                           <  7:50934/101‑1 (MQ=255)
                                                   tCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGaa                                                                                           <  6:279705/149‑1 (MQ=255)
                                                      gTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCaaa                                                                                                 <  6:235422/140‑1 (MQ=255)
                                                      gTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCaaa                                                                                                 >  5:235422/1‑140 (MQ=255)
                                                       tCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATCACCCTGGTCAAATGCGAACtgt                                                                                       <  7:60701/149‑1 (MQ=255)
                                                        cATTGCGTGATAGCCACGATCGATGGCCGCGAAGGTCACTTTGCCGAGCCCTTTGCGGAGCTGGCCGAGGACGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACATGGCCATGATTACCCTGGTCAAATGCGAACtgtg                                                                                      >  4:270402/1‑149 (MQ=255)
                                                         aTTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACtgtg                                                                                      <  6:381909/148‑1 (MQ=255)
                                                          ttGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGtggc                                                                                    >  6:273836/1‑149 (MQ=255)
                                                            gCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGtggctg                                                                                  >  6:278697/1‑149 (MQ=255)
                                                             cGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGtggctgg                                                                                 >  1:105584/1‑149 (MQ=255)
                                                                     ccACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAgggg                                                                         >  8:220893/1‑149 (MQ=255)
                                                                     ccACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAgggg                                                                         >  6:44838/1‑149 (MQ=255)
                                                                     ccACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAggg                                                                          <  5:278697/148‑1 (MQ=255)
                                                                       aCGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCg                                                                       >  5:315156/1‑149 (MQ=255)
                                                                             gATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACg                                                                   <  8:36651/147‑1 (MQ=255)
                                                                             gATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGa                                                                 >  1:228111/1‑149 (MQ=255)
                                                                             gATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGa                                                                 <  2:228111/149‑1 (MQ=255)
                                                                              aTGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGAc                                                                <  5:273836/149‑1 (MQ=255)
                                                                                    gggAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTcgcg                                                          <  6:315156/149‑1 (MQ=255)
                                                                                    gggAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTcgcg                                                          <  7:19499/149‑1 (MQ=255)
                                                                                          gTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTcgcg                                                          >  1:89841/1‑143 (MQ=255)
                                                                                          gTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTcgcg                                                          <  2:89841/143‑1 (MQ=255)
                                                                                          gTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGAACGGACCTCGCGGAGTcc                                                    >  3:205788/1‑149 (MQ=255)
                                                                                              cTTTGCCGCGGCCTTTGCGGCGGGGCCCGTGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGc                                                <  1:206528/149‑1 (MQ=255)
                                                                                                    cGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAgc                                          >  5:386547/1‑149 (MQ=255)
                                                                                                     gAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAgcg                                         >  6:188867/1‑149 (MQ=255)
                                                                                                           tttGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGa                                                                 <  6:338612/119‑1 (MQ=255)
                                                                                                           tttGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGa                                                                 >  5:338612/1‑119 (MQ=255)
                                                                                                               cGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGAACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCa                                <  4:205788/148‑1 (MQ=255)
                                                                                                                  aGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTgg                                    <  6:327492/141‑1 (MQ=255)
                                                                                                                  aGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTgg                                    >  5:327492/1‑141 (MQ=255)
                                                                                                                       gccgaggccgagCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATg                        <  7:7827/148‑1 (MQ=255)
                                                                                                                          gaggccgagCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATggcgg                    >  5:338260/1‑149 (MQ=255)
                                                                                                                           aggccgagCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATggcggc                   >  3:78340/1‑149 (MQ=255)
                                                                                                                             gccgagCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGa                 <  8:97413/149‑1 (MQ=255)
                                                                                                                             gccgagCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGa                 <  6:374000/149‑1 (MQ=255)
                                                                                                                                gagCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGagca              <  3:226896/149‑1 (MQ=255)
                                                                                                                                gagCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGagca              <  2:105584/149‑1 (MQ=255)
                                                                                                                                        ggggCACATCCCAGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCCGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTgg                                    >  6:390117/4‑119 (MQ=255)
                                                                                                                                        ggggCACATCCCAGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCCGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTgg                                    <  5:390117/116‑1 (MQ=255)
                                                                                                                                        gATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGAGCAGCcggcgg      >  8:62466/1‑149 (MQ=255)
                                                                                                                                          tGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGAGCAGCcggcggcg    >  1:170421/1‑149 (MQ=255)
                                                                                                                                           gCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGAGCAGCcggcggcgg   >  8:296021/1‑149 (MQ=255)
                                                                                                                                            cacaTCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGAGCAGCcggcggcgg   <  4:78340/148‑1 (MQ=255)
                                                                                                                                            cacaTCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGAGCAGCCggcggcggc  <  3:389557/149‑1 (MQ=255)
                                                                                                                                             acaTCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGAGCAGCcggcgg      >  2:407020/1‑144 (MQ=255)
                                                                                                                                             acaTCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGAGCAGCcggcgg      <  1:407020/144‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
TGCGCGCGGGGCTGGGCGAGGCGCGAGCGTGCCGGGACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGAGCAGCCGGCGGCGGC  >  NC_002947/885936‑886224

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: