Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,944,530 T→C 8.2% intergenic (+31/+231) 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,5300TC 8.2% 270.1 / 10.3 97intergenic (+31/+231)actP‑I/PP_1744acetate permease/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (41/48);  new base C (4/4);  total (45/52)
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.94e-01

TGGTTCTTCTCGGTCACCGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCATACTGTTTAATCAC  >  NC_002947/1944383‑1944674
                                                                                                                                                    |                                                                                                                                                
tGGTTCTTCTCGGTCACCGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCtt                                                                                                                                                 >  1:44525/1‑149 (MQ=255)
  gTTCTTCTCGGTCACCGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGt                                                                                                                                               >  3:119669/1‑149 (MQ=255)
      ttctCGGTCACCGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCggg                                                                                                                                           >  6:318033/1‑149 (MQ=255)
       tctcGGTCACCGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCgggg                                                                                                                                          >  2:31739/1‑149 (MQ=255)
        ctcGGTCACCGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCgggg                                                                                                                                          <  7:308067/148‑1 (MQ=255)
        ctcGGTCACCGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCgg                                                                                                                                            <  2:44525/146‑1 (MQ=255)
           ggTCACCGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGc                                                                                                                                         <  7:137641/146‑1 (MQ=255)
           ggTCACCGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGc                                                                                                                                         <  5:318033/146‑1 (MQ=255)
           ggTCACCGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGtt                                                                                                                                      <  4:41114/149‑1 (MQ=255)
            gTCACCGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGttt                                                                                                                                     <  7:18159/149‑1 (MQ=255)
               aCCGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGACAAGCGGGGCGttttc                                                                                                                                   <  3:144579/148‑2 (MQ=255)
                ccGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTa                                                                                                                                 <  5:358590/149‑1 (MQ=255)
                 cGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCACAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTAt                                                                                                                                >  1:71929/1‑149 (MQ=255)
                    cAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATCTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTAttt                                                                                                                              <  1:53004/149‑1 (MQ=255)
                    cAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTc                                                                                                                             <  7:274726/149‑1 (MQ=255)
                     aaGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCa                                                                                                                            >  8:188024/1‑149 (MQ=255)
                       gTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATg                                                                                                                          <  7:188024/149‑1 (MQ=255)
                          cAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAg                                                                                                                       >  2:313364/1‑149 (MQ=255)
                           aaGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCGTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTAt                                                                                                                                <  8:127328/139‑1 (MQ=255)
                           aaGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTAt                                                                                                                                >  7:127328/1‑139 (MQ=255)
                           aaGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCa                                                                                                                        >  4:21683/1‑147 (MQ=255)
                           aaGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCa                                                                                                                        <  3:21683/147‑1 (MQ=255)
                                ggcTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGtt                                                                                                                                      <  7:61451/128‑1 (MQ=255)
                                ggcTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGtt                                                                                                                                      >  8:61451/1‑128 (MQ=255)
                                 gcTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTGGGGGCGTTTTTTATTTCATGCAGTCGGCGa                                                                                                                <  1:31739/149‑1 (MQ=255)
                                   tGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATg                                                                                                              >  6:301766/1‑149 (MQ=255)
                                     aGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATg                                                                                                              <  4:119669/147‑1 (MQ=255)
                                           aGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATg                                                                                                                          <  6:107689/129‑1 (MQ=255)
                                           aGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCa                                                                                                       <  3:54520/148‑1 (MQ=255)
                                                  gcTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCgg                                                                                                                   <  4:217788/129‑1 (MQ=255)
                                                  gcTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCgg                                                                                                                   >  3:217788/1‑129 (MQ=255)
                                                  gcTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCATGGGTGCCTGCGGGGAGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGACACGGTTGTCGGGGCGTTTTTTATTTCATGCAGACGGCGATGGCCTCCAGACAgg                                                                                                 >  5:218135/1‑147 (MQ=255)
                                                  gcTGTTCTTACCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAAAGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAgg                                                                                                 <  6:218135/147‑1 (MQ=255)
                                                          tcCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGt                                                                                            >  1:276625/1‑144 (MQ=255)
                                                          tcCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGt                                                                                            <  2:276625/144‑1 (MQ=255)
                                                            ccGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAgg                                                                                                 <  7:311483/137‑1 (MQ=255)
                                                            ccGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAgg                                                                                                 >  8:311483/1‑137 (MQ=255)
                                                            ccGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGGCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGc                                                                                     <  5:301766/149‑1 (MQ=255)
                                                              gCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCAt                                                                                                                           <  4:137479/109‑1 (MQ=255)
                                                              gCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCAt                                                                                                                           >  3:137479/1‑109 (MQ=255)
                                                                     gTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGAc                                                                                                    >  4:227017/1‑125 (MQ=255)
                                                                     gTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGAc                                                                                                    <  3:227017/125‑1 (MQ=255)
                                                                      tGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTccc                                                                           >  4:44581/1‑149 (MQ=255)
                                                                      tGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGACAAGCGGGGCGtttt                                                                                                                                    <  4:187681/92‑1 (MQ=38)
                                                                      tGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGACAAGCGGGGCGtttt                                                                                                                                    <  4:192246/92‑1 (MQ=38)
                                                                      tGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGACAAGCGGGGCGtttt                                                                                                                                    >  3:187681/1‑92 (MQ=38)
                                                                      tGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGACAAGCGGGGCGtttt                                                                                                                                    >  3:192246/1‑92 (MQ=38)
                                                                              cAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGccc                                                                                   <  6:208375/133‑1 (MQ=255)
                                                                              cAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGccc                                                                                   >  5:208375/1‑133 (MQ=255)
                                                                               aGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCg                                                                  >  1:167106/1‑149 (MQ=255)
                                                                                gACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGt                                                                 <  1:313364/149‑1 (MQ=255)
                                                                                 aCTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGt                                                                 <  2:71929/148‑1 (MQ=255)
                                                                                     gCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAgcg                                                            >  5:77128/1‑149 (MQ=255)
                                                                                        tGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCt                                                         >  5:74300/1‑149 (MQ=255)
                                                                                         gggTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCt                                                         >  3:62434/1‑148 (MQ=255)
                                                                                          ggTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCtt                                                        >  4:85950/1‑148 (MQ=255)
                                                                                             gCCAGCGGTGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCa                                                                                  >  6:134157/1‑119 (MQ=255)
                                                                                             gCCAGCGGTGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCa                                                                                  <  5:134157/119‑1 (MQ=255)
                                                                                               cAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTccc                                                                           >  5:312000/1‑124 (MQ=255)
                                                                                               cAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTccc                                                                           <  6:312000/124‑1 (MQ=255)
                                                                                               cAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATgga                                                  <  6:74300/149‑1 (MQ=255)
                                                                                                aGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATggag                                                 >  6:151035/1‑149 (MQ=255)
                                                                                                aGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGACCATGCTCACCTCATGTCGTGAGCGCCTTCATgtag                                                 >  4:203036/1‑149 (MQ=255)
                                                                                                    gggCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACg                                             >  8:215722/1‑149 (MQ=255)
                                                                                                      gCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGc                                           <  6:77128/149‑1 (MQ=255)
                                                                                                         gTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTcc                                                                                                        >  1:222737/1‑85 (MQ=255)
                                                                                                         gTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTcc                                                                                                        <  2:222737/85‑1 (MQ=255)
                                                                                                            tcGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACAttt                                     >  7:43310/1‑149 (MQ=255)
                                                                                                              gCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTgg                                   <  5:151035/149‑1 (MQ=255)
                                                                                                               ctcTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTgg                                   <  3:105546/146‑1 (MQ=255)
                                                                                                                aCTAACATCAATCCCCTGCATTGAGAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTgggg                                 >  1:299001/1‑149 (MQ=255)
                                                                                                                  tAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTgg                                   >  4:105546/1‑145 (MQ=255)
                                                                                                                    aCATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTgg                                   >  5:319102/1‑143 (MQ=255)
                                                                                                                    aCATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTgg                                   <  6:319102/143‑1 (MQ=255)
                                                                                                                    aCATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTg                             <  8:43310/149‑1 (MQ=255)
                                                                                                                     cATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTgg                            <  4:62434/149‑1 (MQ=255)
                                                                                                                              ccTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAg                   <  7:215722/149‑1 (MQ=255)
                                                                                                                                 gCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCa                >  1:175938/1‑149 (MQ=255)
                                                                                                                                 gCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCa                >  2:250948/1‑149 (MQ=255)
                                                                                                                                  cATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCa                >  3:122582/1‑148 (MQ=255)
                                                                                                                                   aataaaaaaCGCCCCGACAAGCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTgg                                   >  2:131639/6‑129 (MQ=38)
                                                                                                                                   aTTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGTTGCTCTGACCATGCTCCCCTCATGTCGTGAGCGCCTTCATGTAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCATa              <  3:203036/149‑1 (MQ=255)
                                                                                                                                   aTTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTgg                                   >  4:120272/1‑128 (MQ=255)
                                                                                                                                   aTTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTgg                                   <  3:120272/128‑1 (MQ=255)
                                                                                                                                    ttGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATggagga                                               >  4:310960/1‑115 (MQ=255)
                                                                                                                                    ttGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATggagga                                               <  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)
                                                                                                                                                    |                                                                                                                                                
TGGTTCTTCTCGGTCACCGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCATACTGTTTAATCAC  >  NC_002947/1944383‑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: