Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 886,084 G→A 5.3% W665* (TGG→TGA) ‡ PP_0767 → LuxR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947886,0840GA 5.3% 230.2 / 2.2 75W665* (TGG→TGA) ‡PP_0767LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (37/34);  new base A (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 = 9.76e-01

GCGCGCGGGGCTGGGCGAGGCGCGAGCGTGCCGGGACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGAGCAGCCGGCGGCGGCACCG  >  NC_002947/885937‑886228
                                                                                                                                                   |                                                                                                                                                
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                                                                                                                                    >  7:171439/1‑134 (MQ=255)
                     gcgAGCGTGCCGGGACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCCCAGATGTGcatcaagc                                                                                                                                    <  8:171439/141‑8 (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                                                                                                                         <  2:76531/139‑1 (MQ=255)
                                  gACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTgg                                                                                                                         >  1:76531/1‑139 (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                                                                                                                                               <  7:50934/101‑1 (MQ=255)
                                                  tCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGAt                                                                                                                                               >  8:50934/1‑101 (MQ=255)
                                                  tCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGaa                                                                                               <  6:279705/149‑1 (MQ=255)
                                                     gTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCaaa                                                                                                     >  5:235422/1‑140 (MQ=255)
                                                     gTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCaaa                                                                                                     <  6:235422/140‑1 (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                                        >  5:327492/1‑141 (MQ=255)
                                                                                                                 aGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTgg                                        <  6:327492/141‑1 (MQ=255)
                                                                                                                      gccgaggccgagCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATg                            <  7:7827/148‑1 (MQ=255)
                                                                                                                         gaggccgagCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATggcgg                        >  5:338260/1‑149 (MQ=255)
                                                                                                                          aggccgagCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATggcggc                       >  3:78340/1‑149 (MQ=255)
                                                                                                                            gccgagCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGa                     <  6:374000/149‑1 (MQ=255)
                                                                                                                            gccgagCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGa                     <  8:97413/149‑1 (MQ=255)
                                                                                                                               gagCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGagca                  <  3:226896/149‑1 (MQ=255)
                                                                                                                               gagCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGagca                  <  2:105584/149‑1 (MQ=255)
                                                                                                                                       ggggCACATCCCAGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCCGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTgg                                        <  5:390117/116‑1 (MQ=255)
                                                                                                                                       ggggCACATCCCAGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCCGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTgg                                        >  6:390117/4‑119 (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)
                                                                                                                                                   ggATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGAGCAGCCGGCGGCGGCACCg  <  5:188867/145‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GCGCGCGGGGCTGGGCGAGGCGCGAGCGTGCCGGGACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGAGCAGCCGGCGGCGGCACCG  >  NC_002947/885937‑886228

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