Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,886,876 A→G 5.9% S1133P (TCA→CCA)  PP_1695 ← sodium‑solute symporter/two‑component system sensor histidine kinase/response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,886,8760AG 5.9% 247.6 / 4.3 87S1133P (TCA→CCA) PP_1695sodium‑solute symporter/two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base A (41/39);  new base G (2/3);  total (44/43)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-01

AAGAGCTTGCCGGCAGACCTGCTGGAGCGGGCGCTGTCCGATGAGGCGCTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCAGCCCGGTGCCGGTCTCGCCGTCGTCCAGGTGGTAGTCCACCAATGCCAGGTGC  >  NC_002947/1886729‑1887024
                                                                                                                                                   |                                                                                                                                                    
aaGAGCTTGCCGGCAGACCTGCTGGAGCGGGCGCTGTCCGATGAGGCGCTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGAt                                                                                                                                                     <  5:223508/149‑1 (MQ=255)
     cTTGCCGGCAGACCTGCTGGAGCGGGCGCTGTCCGATGAGGCGCTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATg                                                                                                                                                >  8:106566/1‑149 (MQ=255)
      ttGCCGGCAGACCTGCTGGAGCGGGCGCTGTCCGATGAGGCGCTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGAT‑GCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTa                                                                                                                                              >  4:385386/1‑149 (MQ=255)
        gCCGGCAGACCTGCTGGAGCGGGCGCTGTCCGATGAGGCGCTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTa                                                                                                                                              <  7:106566/148‑1 (MQ=255)
            gCAGACCTGCTGGAGCGGGCGCTGTCCGATGAGGCGCTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGCTGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGCTGCATGTACCAgg                                                                                                                                         >  7:395948/1‑149 (MQ=255)
             cAGACCTGCTGGAGCGGGCGCTGTCCGATGAGGCGCTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAggg                                                                                                                                        <  3:325327/149‑1 (MQ=255)
               gACCTGCTGGAGCGGGCGCTGTCCGATGAGGCGCTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGGTACATGCATCAGGGct                                                                                                                                      <  4:23887/149‑2 (MQ=255)
                aCCTGCTGGAGCGGGCGCTGTCCGATGAGGCGCTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGa                                                                                                                                     >  2:53384/1‑149 (MQ=255)
                 ccTGCTGGAGCGGGCGCTGTCCGATGAGGCGCTGCCCGAGGTAGTGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGAt                                                                                                                                    >  5:246600/1‑149 (MQ=255)
                 ccTGCTGGAGCGGGCGCTGTCCGATGAGGCGCTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGAt                                                                                                                                    <  4:46449/149‑1 (MQ=255)
                     ctgGAGCGGGCGCTGTCCGATGAGGCGCTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGGTACATGCATCAGGGctggac                                                                                                                                  <  1:173933/147‑6 (MQ=255)
                     ctgGAGCGGGCGCTGTCCGATGAGGCGCTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGtt                                                                                                                                >  7:131998/1‑149 (MQ=255)
                       ggAGCGGGCGCTGTCCGATGAGGCGCTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTc                                                                                                                              >  3:63746/1‑149 (MQ=255)
                         aGCGGGCGCTGTCCGATGAGGCGCTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGTTGCATGTACCAGGGCGATGGTTTCCt                                                                                                                            <  7:107382/149‑1 (MQ=255)
                         aGCGGGCGCTGTCCGATGAGGCGCTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCt                                                                                                                            <  7:21893/149‑1 (MQ=255)
                          gCGGGCGCTGTCCGATGAGGCGCTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCtt                                                                                                                           <  8:45024/149‑1 (MQ=255)
                           cGGGCGCTGTCCGATGAGGCGCTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTtg                                                                                                                          >  7:289401/1‑149 (MQ=255)
                           cGGGCGCTGTCCGATGAGGCGCTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTtg                                                                                                                          <  6:134787/149‑1 (MQ=255)
                               cgcTGTCCGATGAGGCGCTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTTGCCAGGTAGTGCAGCCCTGATGCATGTACCAGGGCGATGGTTTcc                                                                                                                             >  5:374563/1‑142 (MQ=255)
                               cgcTGTCCGATGAGGCGCTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTcc                                                                                                                             <  6:374563/142‑1 (MQ=255)
                                gcTGTCCGATGAGGCGCTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGcc                                                                                                                     >  6:439799/1‑149 (MQ=255)
                                gcTGTCCGATGAGGCGCTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGcc                                                                                                                     <  6:199071/149‑1 (MQ=255)
                                  tGTCCGATGAGGCGCTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGt                                                                                                                   >  3:86245/1‑149 (MQ=255)
                                      cGATGAGGCGCTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCgg                                                                                                                <  8:289401/148‑1 (MQ=255)
                                      cGATGAGGCGCTGCCCGAGGTAGCGGTTACTGAGCCAGGCACAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGgc                                                                                                               >  2:341585/1‑149 (MQ=255)
                                           aGGCGCTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGa                                                                                                          <  8:249903/149‑1 (MQ=255)
                                            ggCGCTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGTAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGa                                                                                                          >  7:3649/1‑148 (MQ=255)
                                             gcgcTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGAt                                                                                                         >  6:275406/1‑148 (MQ=255)
                                             gcgcTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATc                                                                                                        <  1:341585/149‑1 (MQ=255)
                                                cTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCAc                                                                                                      >  7:371428/1‑148 (MQ=255)
                                                cTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCAc                                                                                                      >  5:204998/1‑148 (MQ=255)
                                                   cccGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTa                                                                                                    <  1:397050/147‑1 (MQ=255)
                                                     cGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACAcc                                                                                                <  8:222827/149‑1 (MQ=255)
                                                      gAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGTAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCg                                                                                               <  8:3649/149‑1 (MQ=255)
                                                        ggTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCggg                                                                                             <  5:439799/149‑1 (MQ=255)
                                                          tAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTa                                                                                           <  8:7334/149‑1 (MQ=255)
                                                               gTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGc                                                                                      <  1:53384/149‑1 (MQ=255)
                                                               gTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGc                                                                                      <  8:131998/149‑1 (MQ=255)
                                                                   cTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCa                                                                                   <  1:106859/148‑1 (MQ=255)
                                                                   cTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGAGCTGATCACTACACAGGGTCCCGGCTCAc                                                                                  >  5:100045/1‑149 (MQ=255)
                                                                      aGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCAcc                                                                                 >  2:113736/1‑147 (MQ=255)
                                                                        ccAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGc                                                                             >  6:15109/1‑149 (MQ=255)
                                                                         cAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGcc                                                                            <  7:353134/149‑1 (MQ=255)
                                                                          aGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCATCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCg                                                                           <  6:246600/149‑1 (MQ=255)
                                                                              tCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCg                                                                                               >  8:14353/1‑125 (MQ=255)
                                                                              tCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCg                                                                                               <  7:14353/125‑1 (MQ=255)
                                                                                aGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCa                                                                                   <  4:211150/135‑1 (MQ=255)
                                                                                aGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCa                                                                                   >  3:211150/1‑135 (MQ=255)
                                                                                  aTGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGaa                                                                   <  5:15109/149‑1 (MQ=255)
                                                                                  aTGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGaa                                                                   <  1:333662/149‑1 (MQ=255)
                                                                                      cGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAagcca                                                               <  4:418069/149‑1 (MQ=255)
                                                                                       gATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACAcc                                                                                                <  3:291082/115‑1 (MQ=255)
                                                                                       gATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACAcc                                                                                                >  4:291082/1‑115 (MQ=255)
                                                                                        aTTGTGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAagccagc                                                             >  7:120449/1‑149 (MQ=255)
                                                                                             gCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCAt                                                         >  3:7807/1‑148 (MQ=255)
                                                                                             gCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCAt                                                         >  4:253454/1‑148 (MQ=255)
                                                                                             gCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCAt                                                         >  1:55836/1‑148 (MQ=255)
                                                                                              cAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCAt                                                         >  5:398787/1‑147 (MQ=255)
                                                                                              cAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCa                                                       <  1:428393/149‑1 (MQ=255)
                                                                                                 ggCGCTCACGGCCGCAGGCTTGTCAGGTTTTGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCATcc                                                    <  6:100045/149‑1 (MQ=255)
                                                                                                        aaGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCAGCCCGGTGcc                                             >  1:31851/1‑149 (MQ=255)
                                                                                                         aGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCAGCCCGGTGCCg                                            <  3:385386/149‑1 (MQ=255)
                                                                                                           gCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCAGCCCGGTGCCGGt                                          <  8:371428/149‑1 (MQ=255)
                                                                                                            ccGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCAGCCCGGTGCCGGtc                                         <  7:179890/149‑1 (MQ=255)
                                                                                                            ccGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCAGCCCGGTGCCGGtc                                         >  1:145728/1‑149 (MQ=255)
                                                                                                             cGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACACAGCCGTGCTCGAAGCCAGCCCATCAGCCCGGTGACGGtct                                        >  5:45520/1‑149 (MQ=255)
                                                                                                                 ggCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCAGCCCGGTGCCGGTCTCGcc                                    <  6:204998/149‑1 (MQ=255)
                                                                                                                    ttGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCAGCCCGGTGCCGGTCTCGCcgtc                                 >  7:266329/1‑149 (MQ=255)
                                                                                                                        ccGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCAGCCCGGTACCGGTCTCGCCgtcgtc                              >  3:147950/1‑148 (MQ=255)
                                                                                                                          ggTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCAGCCCGGTACCGGTCTCGCCGTCGTCCAg                           <  4:147950/149‑1 (MQ=255)
                                                                                                                              tGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCAGCCCGGTGCCGGTCTCGCCGTCGTCCAggtgg                       >  5:275086/1‑149 (MQ=255)
                                                                                                                               ggCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCAGCCCGGTGCCGGTCTCGCCGTCGTCCAggtgg                       <  4:7807/148‑1 (MQ=255)
                                                                                                                                  cAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCAGCCCGGTGCCGGTCTCGCCGTCGTCCAGGTGGTAGt                   >  7:237477/1‑149 (MQ=255)
                                                                                                                                  cAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCAGCCCGGTGCCGGTCTCGCCGTCGTCCAGGTGGTAGt                   >  5:283451/1‑149 (MQ=255)
                                                                                                                                     gTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCAGCCCGGTGCCGGTCTCGCCGTCGTCCAGGTGGTAGTcca                >  4:88748/1‑149 (MQ=255)
                                                                                                                                      tAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCAGCCCGGTGCCGGTCTCGCCGTCGTCCAGGTGGTAGTccac               <  8:266329/149‑1 (MQ=255)
                                                                                                                                        gTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCAGCCCGGTGCCGGTCTCGCCGTCGTCCAGGTGGTAGTccacc              >  6:413935/1‑148 (MQ=255)
                                                                                                                                        gTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCAGCCCGGTGCCGGTCGCGCCGTCGTCCAGGTGGTAGTccacc              >  6:338416/1‑148 (MQ=255)
                                                                                                                                         tCCACCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCAGTCGGCGCTGATCACTACACAAGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCAGCCCGGTGCCGGTCTCGCCGTCGTCCAGGTGGTAGTCCACCaa            >  1:131740/1‑149 (MQ=255)
                                                                                                                                         catcgCCCTGGTACATGCATCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCAGCCCGGTGCCGGTCTCGCCGTCGTCCAGGTGGTAGTCCACCaa            >  7:375935/5‑149 (MQ=255)
                                                                                                                                          atcgCCCTGGTACATGCATCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCAGccc                                                   <  8:342460/106‑1 (MQ=255)
                                                                                                                                          atcgCCCTGGTACATGCATCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCAGccc                                                   >  7:342460/4‑109 (MQ=255)
                                                                                                                                           cAGCCCTGATGCATGTCCCACGGCGATGCTTTCATTGCTGCCGTCGGCGCTGATACCTTCACCGGGTACCGGCTCACCCTGCCGAGGTCGAAGCCAGCCCATGCGCCCGGTGCCGGGCTCGCCGTCGTCCAGCTGGTAGTccccaaaag          >  2:319046/1‑142 (MQ=255)
                                                                                                                                              cccTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCAGCCCGGTGCCGGTCTCGCcgtcg                                >  2:408309/1‑124 (MQ=255)
                                                                                                                                              cccTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCAGCCCGGTGCCGGTCTCGCcgtcg                                <  1:408309/124‑1 (MQ=255)
                                                                                                                                              cccTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCAGCCCGGTGCCGGTCTCGCCGTCGTCCAGGTGGTAGTCCACCAATGCCa       <  4:86245/149‑1 (MQ=255)
                                                                                                                                              cccAGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACAGGGTACCGGCTCACCCAGCCGTGCTCCAAGCCAGCCCATCAGCCCGGTGCCGGTCTCGCCGTCGTCCAGGTGGTAGTCCACCAATGc         >  8:279280/1‑147 (MQ=255)
                                                                                                                                                   aTGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCAGCCCGGTGCCGGTCTCGCCGTCGTCCAGGTGGTAGTCCACCAATGCCAGGTGc  <  8:237477/149‑1 (MQ=255)
                                                                                                                                                   aTGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCAGCCCGGTGCCGGTCTCGCCGTCGTCCAGGTGGTAGTCCACCAATGCCAGGTGc  >  7:249839/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AAGAGCTTGCCGGCAGACCTGCTGGAGCGGGCGCTGTCCGATGAGGCGCTGCCCGAGGTAGCGGTTACTGAGCCAGGCTCAGATGGCGATTGAGCAGGGCGCGCAAGGCCGCAGGCTTGACCGGTTTGGCCAGGTAGTCCAGCCCTGATGCATGTACCAGGGCGATGGTTTCCTTGCTGCCGTCGGCGCTGATCACTACACCGGGTACCGGCTCACCCAGCCGTGCTCGAAGCCAGCCCATCAGCCCGGTGCCGGTCTCGCCGTCGTCCAGGTGGTAGTCCACCAATGCCAGGTGC  >  NC_002947/1886729‑1887024

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