Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,944,532 G→A 8.0% intergenic (+33/+229) actP‑I → / ← PP_1744 acetate permease/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,944,5320GA 8.0% 289.4 / 9.1 99intergenic (+33/+229)actP‑I/PP_1744acetate permease/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (42/49);  new base A (4/4);  total (46/53)
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.11e-01

GTTCTTCTCGGTCACCGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCATACTGTTTAATCAC  >  NC_002947/1944385‑1944674
                                                                                                                                                    |                                                                                                                                              
gTTCTTCTCGGTCACCGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGt                                                                                                                                               >  3:119669/1‑149 (MQ=255)
    ttctCGGTCACCGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCggg                                                                                                                                           >  6:318033/1‑149 (MQ=255)
     tctcGGTCACCGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCgggg                                                                                                                                          >  2:31739/1‑149 (MQ=255)
      ctcGGTCACCGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCgggg                                                                                                                                          <  7:308067/148‑1 (MQ=255)
      ctcGGTCACCGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCgg                                                                                                                                            <  2:44525/146‑1 (MQ=255)
         ggTCACCGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGc                                                                                                                                         <  5:318033/146‑1 (MQ=255)
         ggTCACCGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGc                                                                                                                                         <  7:137641/146‑1 (MQ=255)
         ggTCACCGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGtt                                                                                                                                      <  4:41114/149‑1 (MQ=255)
          gTCACCGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGttt                                                                                                                                     <  7:18159/149‑1 (MQ=255)
             aCCGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGACAAGCGGGGCGttttc                                                                                                                                   <  3:144579/148‑2 (MQ=255)
              ccGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTa                                                                                                                                 <  5:358590/149‑1 (MQ=255)
               cGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCACAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTAt                                                                                                                                >  1:71929/1‑149 (MQ=255)
                  cAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATCTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTAttt                                                                                                                              <  1:53004/149‑1 (MQ=255)
                  cAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTc                                                                                                                             <  7:274726/149‑1 (MQ=255)
                   aaGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCa                                                                                                                            >  8:188024/1‑149 (MQ=255)
                     gTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATg                                                                                                                          <  7:188024/149‑1 (MQ=255)
                        cAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAg                                                                                                                       >  2:313364/1‑149 (MQ=255)
                         aaGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCGTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTAt                                                                                                                                <  8:127328/139‑1 (MQ=255)
                         aaGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTAt                                                                                                                                >  7:127328/1‑139 (MQ=255)
                         aaGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCa                                                                                                                        <  3:21683/147‑1 (MQ=255)
                         aaGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCa                                                                                                                        >  4:21683/1‑147 (MQ=255)
                              ggcTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGtt                                                                                                                                      <  7:61451/128‑1 (MQ=255)
                              ggcTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGtt                                                                                                                                      >  8:61451/1‑128 (MQ=255)
                               gcTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTGGGGGCGTTTTTTATTTCATGCAGTCGGCGa                                                                                                                <  1:31739/149‑1 (MQ=255)
                                 tGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATg                                                                                                              >  6:301766/1‑149 (MQ=255)
                                   aGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATg                                                                                                              <  4:119669/147‑1 (MQ=255)
                                         aGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATg                                                                                                                          <  6:107689/129‑1 (MQ=255)
                                         aGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCa                                                                                                       <  3:54520/148‑1 (MQ=255)
                                                gcTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCgg                                                                                                                   >  3:217788/1‑129 (MQ=255)
                                                gcTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCgg                                                                                                                   <  4:217788/129‑1 (MQ=255)
                                                gcTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCATGGGTGCCTGCGGGGAGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGACACGGTTGTCGGGGCGTTTTTTATTTCATGCAGACGGCGATGGCCTCCAGACAgg                                                                                                 >  5:218135/1‑147 (MQ=255)
                                                gcTGTTCTTACCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAAAGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAgg                                                                                                 <  6:218135/147‑1 (MQ=255)
                                                        tcCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGt                                                                                            <  2:276625/144‑1 (MQ=255)
                                                        tcCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGt                                                                                            >  1:276625/1‑144 (MQ=255)
                                                          ccGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAgg                                                                                                 <  7:311483/137‑1 (MQ=255)
                                                          ccGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAgg                                                                                                 >  8:311483/1‑137 (MQ=255)
                                                          ccGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGGCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGc                                                                                     <  5:301766/149‑1 (MQ=255)
                                                            gCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCAt                                                                                                                           >  3:137479/1‑109 (MQ=255)
                                                            gCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCAt                                                                                                                           <  4:137479/109‑1 (MQ=255)
                                                                   gTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGAc                                                                                                    >  4:227017/1‑125 (MQ=255)
                                                                   gTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGAc                                                                                                    <  3:227017/125‑1 (MQ=255)
                                                                    tGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTccc                                                                           >  4:44581/1‑149 (MQ=255)
                                                                    tGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGACAAGCGGGGCGtttt                                                                                                                                    <  4:187681/92‑1 (MQ=38)
                                                                    tGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGACAAGCGGGGCGtttt                                                                                                                                    <  4:192246/92‑1 (MQ=38)
                                                                    tGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGACAAGCGGGGCGtttt                                                                                                                                    >  3:187681/1‑92 (MQ=38)
                                                                    tGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGACAAGCGGGGCGtttt                                                                                                                                    >  3:192246/1‑92 (MQ=38)
                                                                            cAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGccc                                                                                   >  5:208375/1‑133 (MQ=255)
                                                                            cAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGccc                                                                                   <  6:208375/133‑1 (MQ=255)
                                                                             aGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCg                                                                  >  1:167106/1‑149 (MQ=255)
                                                                              gACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGt                                                                 <  1:313364/149‑1 (MQ=255)
                                                                               aCTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGt                                                                 <  2:71929/148‑1 (MQ=255)
                                                                                   gCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAgcg                                                            >  5:77128/1‑149 (MQ=255)
                                                                                      tGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCt                                                         >  5:74300/1‑149 (MQ=255)
                                                                                       gggTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCt                                                         >  3:62434/1‑148 (MQ=255)
                                                                                        ggTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCtt                                                        >  4:85950/1‑148 (MQ=255)
                                                                                           gCCAGCGGTGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCa                                                                                  >  6:134157/1‑119 (MQ=255)
                                                                                           gCCAGCGGTGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCa                                                                                  <  5:134157/119‑1 (MQ=255)
                                                                                             cAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTccc                                                                           <  6:312000/124‑1 (MQ=255)
                                                                                             cAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTccc                                                                           >  5:312000/1‑124 (MQ=255)
                                                                                             cAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATgga                                                  <  6:74300/149‑1 (MQ=255)
                                                                                              aGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATggag                                                 >  6:151035/1‑149 (MQ=255)
                                                                                              aGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGACCATGCTCACCTCATGTCGTGAGCGCCTTCATgtag                                                 >  4:203036/1‑149 (MQ=255)
                                                                                                  gggCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACg                                             >  8:215722/1‑149 (MQ=255)
                                                                                                    gCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGc                                           <  6:77128/149‑1 (MQ=255)
                                                                                                       gTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTcc                                                                                                        <  2:222737/85‑1 (MQ=255)
                                                                                                       gTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTcc                                                                                                        >  1:222737/1‑85 (MQ=255)
                                                                                                          tcGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACAttt                                     >  7:43310/1‑149 (MQ=255)
                                                                                                            gCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTgg                                   <  5:151035/149‑1 (MQ=255)
                                                                                                             ctcTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTgg                                   <  3:105546/146‑1 (MQ=255)
                                                                                                              aCTAACATCAATCCCCTGCAT‑TGAGAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTgggg                                 >  1:299001/1‑149 (MQ=255)
                                                                                                                tAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTgg                                   >  4:105546/1‑145 (MQ=255)
                                                                                                                  aCATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTgg                                   >  5:319102/1‑143 (MQ=255)
                                                                                                                  aCATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTgg                                   <  6:319102/143‑1 (MQ=255)
                                                                                                                  aCATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTg                             <  8:43310/149‑1 (MQ=255)
                                                                                                                   cATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTgg                            <  4:62434/149‑1 (MQ=255)
                                                                                                                            ccTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAg                   <  7:215722/149‑1 (MQ=255)
                                                                                                                               gCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCa                >  1:175938/1‑149 (MQ=255)
                                                                                                                               gCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCa                >  2:250948/1‑149 (MQ=255)
                                                                                                                                cAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCa                >  3:122582/1‑148 (MQ=255)
                                                                                                                                 aataaaaaaCGCCCCGACAAGCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTgg                                   >  2:131639/6‑129 (MQ=38)
                                                                                                                                 aT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGTTGCTCTGACCATGCTCCCCTCATGTCGTGAGCGCCTTCATGTAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCATa              <  3:203036/149‑1 (MQ=255)
                                                                                                                                 aT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTgg                                   <  3:120272/128‑1 (MQ=255)
                                                                                                                                 aT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTgg                                   >  4:120272/1‑128 (MQ=255)
                                                                                                                                  t‑tGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATggagga                                               >  4:310960/1‑115 (MQ=255)
                                                                                                                                  t‑tGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATggagga                                               <  3:310960/115‑1 (MQ=255)
                                                                                                                                  ataaaaaaCGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCATa              <  5:329046/145‑1 (MQ=255)
                                                                                                                                   taaaaaaCGCCCCGACAAGCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGc                  <  4:60102/141‑1 (MQ=255)
                                                                                                                                   taaaaaaCGCCCCGACAAGCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGc                  >  3:60102/4‑144 (MQ=255)
                                                                                                                                     gAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCATACTg           >  8:8195/1‑149 (MQ=255)
                                                                                                                                      aaaaCGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTg                                                                                              <  4:174879/65‑1 (MQ=255)
                                                                                                                                      aaaaCGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCATACTGt          <  7:8195/149‑1 (MQ=255)
                                                                                                                                       aaaCGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCATACTGt          <  3:85950/148‑1 (MQ=255)
                                                                                                                                       aaaCGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCATACTGt          <  2:167106/148‑1 (MQ=255)
                                                                                                                                           gCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCATACTGTTTaa      >  2:286818/1‑148 (MQ=255)
                                                                                                                                           gCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCATACTGTTTaa      <  1:286818/148‑1 (MQ=255)
                                                                                                                                              ccGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCATACTGTTTAATCAc  >  6:5659/1‑149 (MQ=255)
                                                                                                                                              ccGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCATACTGTTTAATCAc  <  1:250948/149‑1 (MQ=255)
                                                                                                                                                  ttGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGc                  <  6:25379/129‑1 (MQ=255)
                                                                                                                                                  ttGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGc                  >  5:25379/1‑129 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GTTCTTCTCGGTCACCGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCAT‑TGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCATACTGTTTAATCAC  >  NC_002947/1944385‑1944674

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