Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,825,280 T→A 6.4% T256S (ACC→TCC)  ptxS ← 2‑ketogluconate utilization repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,825,2800TA 6.4% 179.8 / 2.6 63T256S (ACC→TCC) ptxS2‑ketogluconate utilization repressor
Reads supporting (aligned to +/- strand):  ref base T (31/28);  new base A (2/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.53e-01

GCCACGCCGATACGCTCGGTGGGCTGGGCCAGTGCGGTGATGCCTTTGCCGACCAAGGGGTACCAGTCGAGGTCGTCGAGGGCGATCAGCCCGACATCGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCATGCCTGGCCGTCGGCTGATCGACGCGCCGAATGCCTG  >  NC_002947/3825132‑3825427
                                                                                                                                                    |                                                                                                                                                   
gccACGCCGATACGCTCGGTGGGCTGGGCCAGTGCGGTGATGCCTTTGCCGACCAAGGGGTACCAGTCGAGGTCGTCGAGGGCGATCAGCCCGACATCGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGt                                                                                                                                                     >  7:330952/1‑149 (MQ=255)
  cACGCCGATACGCTCGCTGGGCTGGGCCAGTGCGGTGATGCCTTTGCCGACCAAGGGGTACCAGTCGAGGTCGTCGAGGGCGATCAGCCCGACATCGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTgg                                                                                                                                                   <  4:382210/149‑1 (MQ=255)
            ctcTCGGTGGGCTGGGGCAGTGCGGTGATGCCTTTGCCGACCAAGGGGTACCAGTCGAGGTCGTCGAGGGGGATCAGCCCGACATCGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGtt                                                                                                                                         <  1:233616/147‑1 (MQ=255)
            cGCTCGGTGGGCTGGGCCAGTGCGGTGATGCCTTTGCCGACCAAGGGGTACCAGTCGAGGTCGTCGAGGGCGATCAGCCCGACATCGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGt                                                                                                                                          <  7:381068/148‑1 (MQ=255)
                              aGTGCGGTGATGCCTTTGCCGACCAAGGGGTACCAGTCGAGGTCGTCGAGGGCGATCAGCCCGACATCGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTg                                                                                                                       >  2:158980/1‑149 (MQ=255)
                                     tGATGCCTTTGCCGACCAAGGGGTACCAGTCGAGGTCGTCGAGGGCGATCAGCCCGACATCGGCAACGAGGTTGCGCCCGGCCTCGAGCAGCGCCCTGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCa                                                                                                                <  7:393092/149‑1 (MQ=255)
                                           cTTTGCCGACCAAGGGGTACCAGTCGAGGTCGTCGAGGGCGATCAGCCCGACATCGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCa                                                                                                          >  2:325435/1‑149 (MQ=255)
                                                ccGACCAAGGGGTACCAGTCGAGGTCGTCGAGGGCGATCAGCCCGACATCGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACtgt                                                                                                      <  4:297163/148‑1 (MQ=255)
                                                 cGACCAAGGGGTACCAGTCGAGGTCGTCGAGGGCGATCAGCCCGACATCGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACtgtg                                                                                                     <  8:330952/148‑1 (MQ=255)
                                                  gACCAAGGGGTACCAGTCGAGGTCGTCGAGGGCGATCAGCCCGACATCGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTggg                                                                                                   <  1:191924/149‑1 (MQ=255)
                                                  gACCAAGGGGTACCAGTCGAGGTCGTCGAGGGCGATCAGCCCGACATCGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTggg                                                                                                   >  7:192076/1‑149 (MQ=255)
                                                     cAAGGGGTACCGGTCGAGGTCGTCGAGGGCGATCAGCCCGACATCGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGAGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGc                                                                                                  >  2:334659/1‑147 (MQ=255)
                                                       aGGGGTACCAGTCGAGGTCGTCGAGGGCGATCAGCCCGACATCGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGa                                                                                              <  7:158336/149‑1 (MQ=255)
                                                         gggTACCAGTCGAGGTCGTCGAGGGCGATCAGCCCGACATCGGCAACCAGGTTGCGCCCGGCCTCAAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGaaa                                                                                            >  5:219265/1‑149 (MQ=255)
                                                              ccggTCGAGGTCGTCGAGGGCGATCAGCCCGACATCGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGAGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGaa                                                                                        <  1:334659/145‑1 (MQ=255)
                                                               cAGTCGAGGTCGTCGAGGGCGATCAGCCCGACATCCGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGc                                                                                      >  3:155475/1‑149 (MQ=255)
                                                                     aGGTCGTCGAGGGCGATCAGCCCGACATCGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGCCCACTGTGGGCAAGAAACGAAGCCAGATg                                                                                >  7:395633/1‑149 (MQ=255)
                                                                      ggTCGTCGAGGGCGATCAGCCCGACATCGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTgg                                                                               >  7:311193/1‑149 (MQ=255)
                                                                          gtcgAGGGCGATCAGCCCGACATCGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGaa                                                                                        >  4:188990/1‑136 (MQ=255)
                                                                          gtcgAGGGCGATCAGCCCGACATCGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGaa                                                                                        <  3:188990/136‑1 (MQ=255)
                                                                          gtcgAGGGCGATCAGCCCGACATCGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCt                                                                           <  2:308640/149‑1 (MQ=255)
                                                                          gtcgAGGGCGATCAGCCCGACATCGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCt                                                                           <  6:198183/149‑1 (MQ=255)
                                                                             gAGGGCGATCAGCCCGACATCGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCa                                                                        >  2:197659/1‑149 (MQ=255)
                                                                             gAGGGCGATCAGCCCGACATCGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGAGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTg                                                                          >  7:70175/1‑147 (MQ=255)
                                                                              aGGGCGATCAGCGCGACATCGGCAACCTGGTTGCGGCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAg                                                                       <  1:187402/149‑1 (MQ=255)
                                                                                      cAGCCCGACATCGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGgagc                                                               >  6:123165/1‑149 (MQ=255)
                                                                                       aGCCCGACATCGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGgcgc                                                               <  1:197659/148‑1 (MQ=255)
                                                                                           cGACATCGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAtt                                                          >  8:216894/1‑149 (MQ=255)
                                                                                           cGACATCGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAtt                                                          >  8:98889/1‑149 (MQ=255)
                                                                                             aCATCGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACTACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTc                                                        <  2:118385/149‑1 (MQ=255)
                                                                                               aTCGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTc                                                        >  7:365650/1‑147 (MQ=255)
                                                                                                 cGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTcc                                                       <  5:148415/146‑1 (MQ=255)
                                                                                                 cGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTcc                                                       >  6:148415/1‑146 (MQ=255)
                                                                                                 cGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAg                                                     >  7:9119/1‑148 (MQ=255)
                                                                                                   gCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTa                                                   <  4:410331/148‑1 (MQ=255)
                                                                                                   gCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTa                                                   >  3:410331/1‑148 (MQ=255)
                                                                                                   gCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTAc                                                  >  1:394301/1‑149 (MQ=255)
                                                                                                   gCAACCAGGGTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTAc                                                  <  4:116250/149‑1 (MQ=255)
                                                                                                    cAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTAcc                                                 >  5:27350/1‑149 (MQ=255)
                                                                                                     aaCCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTAc                                                  <  1:253854/147‑1 (MQ=255)
                                                                                                     aaCCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTAc                                                  >  2:253854/1‑147 (MQ=255)
                                                                                                          ggTTGCGCCCGGTCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGa                                           >  7:251394/1‑149 (MQ=255)
                                                                                                           gtggcgcCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGAc                                          <  2:407711/146‑1 (MQ=255)
                                                                                                               cgcCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCAt                                      >  6:226638/1‑149 (MQ=255)
                                                                                                                tcccGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCATg                                     >  3:50795/2‑149 (MQ=255)
                                                                                                                gcCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACctgct                                             <  7:53601/141‑1 (MQ=255)
                                                                                                                gcCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACctgct                                             >  8:53601/1‑141 (MQ=255)
                                                                                                                   cGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCATGCCt                                  <  1:325435/149‑1 (MQ=255)
                                                                                                                     gCCTCAAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCATGCCTgg                                <  6:219265/149‑1 (MQ=255)
                                                                                                                          gAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCATGCCTGGCCGt                            <  2:394301/148‑1 (MQ=255)
                                                                                                                            gcagcGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCATGCCTGGCCGTCgg                         <  8:251394/149‑1 (MQ=255)
                                                                                                                               gcgcCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCATGCCTGGCCGTCGGCt                       >  1:72811/1‑148 (MQ=255)
                                                                                                                                 gcCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCATGCCTGGCCGTCGGCTGAt                    <  8:159074/149‑1 (MQ=255)
                                                                                                                                       gTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCATGCCTGGc                               <  2:425643/132‑1 (MQ=255)
                                                                                                                                       gTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCATGCCTGGc                               >  1:425643/1‑132 (MQ=255)
                                                                                                                                       gTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCATGCCTGGCCGTCGGCTGATCGAcgc              <  7:335187/149‑1 (MQ=255)
                                                                                                                                       gTGACCGCCAAGGAGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCATGCCTGGCCGTCGGCTGATCGAcgc              <  8:70175/149‑1 (MQ=255)
                                                                                                                                        tGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCATGCCTGGCCGTCGGCTGATCGAcgcg             >  2:12833/1‑149 (MQ=255)
                                                                                                                                            cGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCATGCCTGGCCGTCGGCTGATCGACGCGCCg          >  5:292663/1‑148 (MQ=255)
                                                                                                                                            cGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCATGCCTGGCCGTCGGCTGATCGACGCGCCg          >  3:95635/1‑148 (MQ=255)
                                                                                                                                              ccAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCATGCCTGGCCGTCGGCTGATCGACGCGCCGAAt       >  3:354311/1‑149 (MQ=255)
                                                                                                                                               cAAGGTGGCAACGCCGTTGAAGGTGTAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCATGCCTGGCCGTCGGCTGAt                    <  5:226695/135‑1 (MQ=255)
                                                                                                                                               cAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCATGCCTGGCCGTCGGCTGAt                    >  6:226695/1‑135 (MQ=255)
                                                                                                                                                  ggTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCATGCCTGGCCGTCGGCTGATCGACGCGCCGAATGCCt   >  4:218407/1‑149 (MQ=255)
                                                                                                                                                    tGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCATGCCTGGCCGTCGGCTGATCGACGCGCCGAATGCCTg  >  5:433633/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCCACGCCGATACGCTCGGTGGGCTGGGCCAGTGCGGTGATGCCTTTGCCGACCAAGGGGTACCAGTCGAGGTCGTCGAGGGCGATCAGCCCGACATCGGCAACCAGGTTGCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCATGCCTGGCCGTCGGCTGATCGACGCGCCGAATGCCTG  >  NC_002947/3825132‑3825427

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