Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,672,934 T→G 5.3% Q25H (CAA→CAC PP_4136 ← LuxR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,672,9340TG 5.3% 245.9 / 4.0 76Q25H (CAA→CACPP_4136LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (40/32);  new base G (2/2);  total (42/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.90e-01

CTTGATCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAGGC  >  NC_002947/4672793‑4673076
                                                                                                                                             |                                                                                                                                                 
cTTGATCACAGCAGCCACCCGGGTAGTCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGt                                                                                                                                             >  1:198921/1‑148 (MQ=255)
 ttGATCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTc                                                                                                                                           >  5:87747/1‑149 (MQ=255)
    aTCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCAAGTGCTGttgtg                                                                                                                                         <  6:158509/148‑4 (MQ=255)
    aTCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCAAGTGCTGttgtg                                                                                                                                         <  7:322765/148‑4 (MQ=255)
    aTCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGAACTTGCTGTTCACg                                                                                                                                        >  8:20345/1‑149 (MQ=255)
     tCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATTACAGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGt                                                                                                                                       <  1:123572/149‑1 (MQ=255)
       acaGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTg                                                                                                                                      <  5:82236/148‑1 (MQ=255)
         agcagcCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCg                                                                                                                                   >  1:100139/1‑149 (MQ=255)
          gcaggCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGACCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGt                                                                                                                                  <  5:122342/149‑1 (MQ=255)
             gcCACCCGGGTAGGCACTTCAAATTTGAGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAg                                                                                                                               >  6:306920/1‑149 (MQ=255)
                aCCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGc                                                                                                                             >  5:176490/1‑148 (MQ=255)
                  ccGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGaa                                                                                                                          <  3:151860/149‑1 (MQ=255)
                       tAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCggg                                                                                                                     <  7:14140/149‑1 (MQ=255)
                       tAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCggg                                                                                                                     <  5:306920/149‑1 (MQ=255)
                       tAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCggg                                                                                                                     <  2:8604/149‑1 (MQ=255)
                          gCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCg                                                                                                                  >  1:222096/1‑149 (MQ=255)
                          gCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCg                                                                                                                  >  3:239143/1‑149 (MQ=255)
                              ttaaaaTTTGCGCAGGATGTTGGATACGTGAAAATTCAACGTCGACTCCTTCCAATCAAGGATCTGGCCGATCTACCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACTTGGCGTCAGGCTAAACGGGGCGTCCg                                                                                                              >  2:86570/4‑149 (MQ=255)
                              ttCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCg                                                                                                              >  8:214801/1‑149 (MQ=255)
                              ttCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGACGATCTCGCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGACCg                                                                                                              >  2:81592/1‑149 (MQ=255)
                                cAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGc                                                                                                             >  1:179978/1‑148 (MQ=255)
                                    tttGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCt                                                                                                                                               <  4:92024/110‑1 (MQ=255)
                                    tttGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCt                                                                                                                                               >  3:92024/1‑110 (MQ=255)
                                        cgcAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGtt                                                                                                     <  2:100139/148‑1 (MQ=255)
                                         gcAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTtgc                                                                                                   >  7:289089/1‑149 (MQ=255)
                                            ggATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTtgctgc                                                                                                <  6:176490/149‑1 (MQ=255)
                                              aTGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGa                                                                                              <  2:222096/149‑1 (MQ=255)
                                              aTGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGa                                                                                              <  7:103863/149‑1 (MQ=255)
                                               tGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAg                                                                                             >  2:180464/1‑149 (MQ=255)
                                                 ttGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCCCGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAgtg                                                                                           >  4:141608/1‑149 (MQ=255)
                                                 ttGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAgtg                                                                                           <  1:328345/149‑1 (MQ=255)
                                                     aTACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGa                                                                                       <  4:239143/149‑1 (MQ=255)
                                                      tACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAg                                                                                      <  7:214801/149‑1 (MQ=255)
                                                         gTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCCTCCGCCAGGGGTTGCTGCGAGTGTCGGGCGt                                                                                   >  8:351426/1‑149 (MQ=255)
                                                              aaTTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTgtca                                                                              <  3:124740/149‑1 (MQ=255)
                                                                        cGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTgt                                                                                >  8:23549/1‑137 (MQ=255)
                                                                        cGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTgt                                                                                <  7:23549/137‑1 (MQ=255)
                                                                          aCTCCTTGCATTCAAGGATCTGGCCGTTCTGACACGATCTCTTGCCATAAGCACACGAAAACAGGAGTTGCTGTTCACGTGGCGTCAGGCGATCCGGGGCGACCGGCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGaa                                                                  <  1:86570/149‑1 (MQ=255)
                                                                           cTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACGACAACAGGAGTTGCTGTTCACGTGGCGTCAGGCGAAGCGGGGCGACCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGaa                                                                  <  1:81592/148‑1 (MQ=255)
                                                                            tCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑t                                                              >  4:201977/1‑149 (MQ=255)
                                                                              cTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑Tc                                                            <  1:180464/148‑1 (MQ=255)
                                                                                tGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCAc                                                         >  5:205397/1‑149 (MQ=255)
                                                                                   aTTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACggg                                                      >  1:216751/1‑149 (MQ=255)
                                                                                    ttCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCTCTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGt                                                     <  4:158897/149‑1 (MQ=255)
                                                                                    ttCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGt                                                     <  2:86864/149‑1 (MQ=255)
                                                                                    ttCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGt                                                     >  7:63302/1‑149 (MQ=255)
                                                                                     tCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGt                                                     >  4:262240/1‑148 (MQ=255)
                                                                                       aaGGATCTGGCCGATCTCCGACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGGGTACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGTGAt                                                  <  7:351426/149‑1 (MQ=255)
                                                                                       aaGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCCGGCGAACCGGGGCGTCCGCGAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGTGAt                                                  >  6:216650/1‑149 (MQ=255)
                                                                                       aaGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGTGAt                                                  >  7:313242/1‑149 (MQ=255)
                                                                                                 ccGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGTGATTGAGCGGGGc                                        <  8:313242/149‑1 (MQ=255)
                                                                                                   gATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGTGATTGAGCGGGGCAt                                      >  7:154153/1‑149 (MQ=255)
                                                                                                         ccACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGTGATTGAGCGGGGCATtgctgc                                >  6:300719/1‑149 (MQ=255)
                                                                                                         ccACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGTGATTGAGCGGGGCATtgctgc                                >  6:317814/1‑149 (MQ=255)
                                                                                                         ccACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGGCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGTGATTGAGCGGGGCATtgct                                  <  6:87747/147‑1 (MQ=255)
                                                                                                          cACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGTGATTGAGCGGGGCATTGCTGcc                               >  1:211864/1‑149 (MQ=255)
                                                                                                           aCGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGTGATTGAGCGGGGCATTGCTGcc                               >  6:38265/1‑148 (MQ=255)
                                                                                                            cGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGTGATTGAGCGGGGCATTGCTGCCTa                             >  1:116289/1‑149 (MQ=255)
                                                                                                               gCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCt                          >  6:94353/1‑149 (MQ=255)
                                                                                                                ctctTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTg                         <  3:262240/149‑1 (MQ=255)
                                                                                                                      ccATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCa                   >  5:286490/1‑149 (MQ=255)
                                                                                                                       cATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCAc                  >  1:42636/1‑149 (MQ=255)
                                                                                                                         tAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCt                <  6:286490/149‑1 (MQ=255)
                                                                                                                         tAAGCACACCACAACAGCACTTGCTGTTCACGTGCCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGGTGTCAGTCAGGTTTGAAC‑‑TCCACGGGAGATTGAGCGCGGCATTGCTGCCTTCCTGCAGTCCCCt                >  8:278216/1‑149 (MQ=255)
                                                                                                                               caccacAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGaa          <  8:63302/149‑1 (MQ=255)
                                                                                                                               caccacAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGaa          <  5:300719/149‑1 (MQ=255)
                                                                                                                                accacAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAg         >  4:79802/1‑149 (MQ=255)
                                                                                                                                accacAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAg         >  2:86126/1‑149 (MQ=255)
                                                                                                                                accacAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAg         >  6:99748/1‑149 (MQ=255)
                                                                                                                                  gtgaacaGCAAGTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACAGTGAGGCGGGTGATTGAGCGGGGAATTGCAGCCTACCTGCAGTCACCTTACGaa          >  1:303884/4‑149 (MQ=255)
                                                                                                                                  cacaACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGTGATTGAGCGGGGCATtgctgc                                >  4:145898/1‑124 (MQ=255)
                                                                                                                                  cacaACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGTGATTGAGCGGGGCATtgctgc                                <  3:145898/124‑1 (MQ=255)
                                                                                                                                   tgaacaGCAAGTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGca       >  2:204764/3‑148 (MQ=255)
                                                                                                                                     aacaGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAg    <  3:209527/149‑1 (MQ=255)
                                                                                                                                     aacaGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAg    <  2:179978/149‑1 (MQ=255)
                                                                                                                                       caGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAGGc  >  5:147548/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
CTTGATCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAC‑‑TCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAGGC  >  NC_002947/4672793‑4673076

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: