Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,584,611 C→T 7.1% intergenic (+98/+29) PP_2264 → / ← folD‑II sugar ABC transporter substrate‑binding protein/bifunctional methylenetetrahydrofolate dehydrogenase/methenyltetrahydrofolate cyclohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,584,6110CT 7.1% 245.0 / 6.1 84intergenic (+98/+29)PP_2264/folD‑IIsugar ABC transporter substrate‑binding protein/bifunctional methylenetetrahydrofolate dehydrogenase/methenyltetrahydrofolate cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base C (39/39);  new base T (3/3);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.53e-01

GTCAGCTACAACGAACTGCTGAAGTCCTGGGAGGCGGCCAGGAAGTGATCGGCTGAACGCCGGGGTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCGACGCCGCCCGGCACCGGTGTGATCCAGCCGGCGCGGGGCAGGGCAGTCTC  >  NC_002947/2584466‑2584756
                                                                                                                                                 |                                                                                                                                                  
gTCAGCTACAACGAACTGCTGAAGTCCTGGGAGGCGGCCAGGAAGTGATCGGCTGAACGCCGGGGTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCaaa                                                                                                                                                 <  6:340000/149‑1 (MQ=255)
         aaCGAACTGCTGAAGTCCTGGGAGGCGGCCAGGAAGTGATCGGCTGAACGCCGGGGTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGtt                                                                                                                                        >  7:203186/1‑149 (MQ=255)
          aCGAACTGCTGAAGTCCTGGGAGGCGGCCAGGAAGTGATCGGCTGAACGCCGGGGTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGttt                                                                                                                                       <  6:204498/149‑1 (MQ=255)
            gAACTGCTGAAGTCCTGGGAGGCGGCCAGGAAGTGATCGGCTGAACGCCGGGGTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGCCAAAAAAAACGTCACCTGCAAAGGTGCCGTTTTc                                                                                                                                     >  3:251369/1‑149 (MQ=255)
                tgctgAAGTCCTGGGAGGCGGCCAGGAAGTGATCGGCTGAACGCCGGGGTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGc                                                                                                                                   <  8:124782/147‑1 (MQ=255)
                tgctgAAGTCCTGGGAGGCGGCCAGGAAGTGATCGGCTGAACGCCGGGGTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGc                                                                                                                                   >  7:124782/1‑147 (MQ=255)
                  ctgAAGTCCTGGGAGGCGGCCAGGAAGTGATCGGCTGAACGCCGGGGTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTTTGCAGGTGCCGttttttttg                                                                                                                                 <  3:337027/147‑6 (MQ=255)
                  ctgAAGTCCTGGGAGGCGGCCAGGAAGTGATCGGCTGAACGCCGGGGTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCa                                                                                                                               >  3:74857/1‑149 (MQ=255)
                   tgAAGTCCTGGGAGGCGGCCAGGAAGTGATCGGCTGAACGCCGGGGTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCa                                                                                                                               >  7:86213/1‑148 (MQ=255)
                   tgAAGTCCTGGGAGGCGGCCAGGAAGTGATCGGCTGAACGCCGGGGTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCAc                                                                                                                              >  3:216189/1‑149 (MQ=255)
                    gAAGTCCTGGGAGGCGGCCAGGAAGTGATCGGCTGAACGCCGGGGTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTTTGCAGGTGCCGtt                                                                                                                                        >  1:134653/1‑138 (MQ=255)
                    gAAGTCCTGGGAGGCGGCCAGGAAGTGATCGGCTGAACGCCGGGGTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTTTGCAGGTGCCGtt                                                                                                                                        <  2:134653/138‑1 (MQ=255)
                     aaGTCCTGGGAGGCGGCCAGGAAGTGATCGGCTGAACGCCGGGGTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCa                                                                                                                               <  8:86213/146‑1 (MQ=255)
                     aaGTCCTGGGAGGCGGCCAGGAAGTGATCGGCTGAACGCCGGGGTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCAc                                                                                                                              <  2:202122/147‑1 (MQ=255)
                          cTGGGAGGCGGCCAGGAAGTGATCGGCTGAACGCCGGGGTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTAt                                                                                                                       <  8:22678/149‑1 (MQ=255)
                             ggAGGCGGCCAGGAAGTGATCGGCTGAACGCCGGGGTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTAc                                                                                                                    <  4:216189/149‑1 (MQ=255)
                                 gcggcCAGGAAGTGATCGGCTGAACGCCGGGGTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTtg                                                                                                                 <  4:74857/148‑1 (MQ=255)
                                  cggcCAGGAAGTGATCGGCTGAACGCCGGGGTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTtgtg                                                                                                               <  5:51560/149‑1 (MQ=255)
                                  cggcCAGGAAGTGATCGGCTGAACGCCGGGGTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTtgtg                                                                                                               <  5:372/149‑1 (MQ=255)
                                   ggcCAGGAAGTGATCGGCTGAACGCCGGGGTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACg                                                                                                                             <  4:241325/134‑1 (MQ=255)
                                   ggcCAGGAAGTGATCGGCTGAACGCCGGGGTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACg                                                                                                                             >  3:241325/1‑134 (MQ=255)
                                      cAGGAAGTGATCGGCTGAACGCCGGGGTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGc                                                                                                           >  6:137385/1‑149 (MQ=255)
                                            gTGATCGGCTGAACGCCGGGGTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCtcttct                                                                                                     <  1:273944/149‑1 (MQ=255)
                                                tCGGCTGAACGCCGGGGTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCg                                                                                                 >  6:86785/1‑149 (MQ=255)
                                                  ggCTGAACGCCGGGGTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGc                                                                                                   <  6:43186/145‑1 (MQ=255)
                                                  ggCTGAACGCCGGGGTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGc                                                                                                   >  5:43186/1‑145 (MQ=255)
                                                        aCGCCGGGGTGGTATCACCGACCGCGGGGGCCGTGCAACCGCTCGCGCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGc                                                                                         <  2:311020/149‑1 (MQ=255)
                                                          gCCGGGGTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCgt                                                                                      <  8:203186/149‑1 (MQ=255)
                                                             ggggTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGtt                                                                                   <  6:246849/149‑1 (MQ=255)
                                                              gggTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCTCGCAGTATTACTTGTGCAGCTCTTCTGcc                                                                                                  >  4:51800/1‑134 (MQ=255)
                                                              gggTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCTCGCAGTATTACTTGTGCAGCTCTTCTGcc                                                                                                  <  3:51800/134‑1 (MQ=255)
                                                                   gtTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAATAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCca                                                                              <  2:32588/148‑1 (MQ=255)
                                                                   gtTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCca                                                                              <  3:327350/148‑1 (MQ=255)
                                                                    tttCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCcagc                                                                            >  3:330673/1‑149 (MQ=255)
                                                                    tttCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCcag                                                                             >  4:321615/1‑148 (MQ=255)
                                                                     ttCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCcag                                                                             >  5:248419/1‑147 (MQ=255)
                                                                        accgaccgCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAgg                                                                         >  6:191326/1‑148 (MQ=255)
                                                                        accgaccgCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAgg                                                                         >  1:212326/1‑148 (MQ=255)
                                                                        accgaccgCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAgg                                                                         >  2:248089/1‑148 (MQ=255)
                                                                        accgaccgCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAgg                                                                         >  7:225011/1‑148 (MQ=255)
                                                                          cgaccgCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCg                                                                                                 <  7:101004/123‑1 (MQ=255)
                                                                          cgaccgCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCg                                                                                                 >  8:101004/1‑123 (MQ=255)
                                                                          cgaccgCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCAc                                                                      <  1:326367/149‑1 (MQ=255)
                                                                                       cGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCg                                                         >  1:130052/1‑149 (MQ=255)
                                                                                          gCAACTGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGcc                                                      >  5:143670/1‑149 (MQ=255)
                                                                                           cAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGccc                                                     >  2:285399/1‑149 (MQ=255)
                                                                                            aaCCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCg                                                    >  2:303311/1‑149 (MQ=255)
                                                                                             aCCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAg                                                                                                            <  2:190462/93‑1 (MQ=255)
                                                                                             aCCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAg                                                                                                            >  1:190462/1‑93 (MQ=255)
                                                                                               cgctcgctCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCGAcg                                                 >  8:198322/1‑149 (MQ=255)
                                                                                                gctcgctcGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCGAcgc                                                <  3:321615/149‑1 (MQ=255)
                                                                                                          gAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCGACGCCGCCCGGCAc                                      <  8:225011/149‑1 (MQ=255)
                                                                                                          gAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCGACGCCGCCCGGCAc                                      <  5:191326/149‑1 (MQ=255)
                                                                                                            aGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATAGAGCGTGTTCTCCAGCAGGCACGCCCGGGt                                                              >  3:302773/1‑123 (MQ=255)
                                                                                                            aGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGt                                                              <  4:297682/123‑1 (MQ=255)
                                                                                                            aGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGt                                                              >  3:297682/1‑123 (MQ=255)
                                                                                                            aGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCGACGCCGCCCGGCACCg                                    <  7:198322/149‑1 (MQ=255)
                                                                                                            aGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATAAAGCGTGTTCTCCAGCAGGCACGCCCGGGt                                                              <  4:302773/123‑1 (MQ=255)
                                                                                                                gCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGccc                                                                  >  4:311277/1‑115 (MQ=255)
                                                                                                                  cACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCGACGCCGCCCGGCACCGGTGTGa                              <  2:130052/149‑1 (MQ=255)
                                                                                                                     tAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCGACGCCGCCCGGCACCGGTGTGATc                            >  5:151963/1‑148 (MQ=255)
                                                                                                                     tAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCGACGCCGCCCGGCACCGGTGTGATc                            >  5:45053/1‑148 (MQ=255)
                                                                                                                       gTCACGGGCAAAAAAAACTGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCGACGCCGCCCGGCACCGGTGTGATCCa                          <  4:251369/148‑1 (MQ=255)
                                                                                                                       gTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGCCGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCGACGCCGCCCGGCACCGGTGTGATCCa                          <  1:285399/149‑1 (MQ=255)
                                                                                                                       gTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCcagca                                                                           >  2:26543/1‑99 (MQ=255)
                                                                                                                       gTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCcagca                                                                           <  1:26543/99‑1 (MQ=255)
                                                                                                                        tCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCGACGCCGCCCGGCACCGGTGTGATCCAGc                        >  5:51603/1‑149 (MQ=255)
                                                                                                                          aCGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCGACGCCGCCCGGCACCGGTGTGATCCAGcc                       <  1:248089/148‑1 (MQ=255)
                                                                                                                          aCGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCGACGCCGCCCGGCACCGGTGTGATCCAGCCg                      >  4:333276/1‑149 (MQ=255)
                                                                                                                          aCGGGCAAAAAAAACCGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCGCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATc                                                          <  4:259856/113‑1 (MQ=255)
                                                                                                                               cAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCGACGCCGCCCGGCACCGGTGTGATCCAGCCGGcgcg                 >  6:47943/1‑149 (MQ=255)
                                                                                                                                tgcgaaaaCGGCACCTTTGCAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCGACGCCGCCCGGCACCGGTGTGATCCAGCCGGCGCgg                >  5:163466/5‑149 (MQ=255)
                                                                                                                                aaaaaaaaCGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCGACGCCGCCCGGCACCGGTGTGATCCAGCCGGCGCgg                <  2:146568/149‑1 (MQ=255)
                                                                                                                                aaaaaaaaCGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCGACGCCGCCCGGCACCGGTGTGATCCAGCCGGCGCgg                >  1:146568/1‑149 (MQ=255)
                                                                                                                                  cgaaaaCGGCACCTTTGCAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCGACGCCGCCCGGCACCGGTGTGATCCAGCCGGCGCgggg              <  6:163466/147‑1 (MQ=255)
                                                                                                                                  aaaaaaCGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCGACGCCGCCCGGCACCGGTGTGATCCAGCCGGCGCgggg              >  6:14195/1‑149 (MQ=255)
                                                                                                                                    aaaaCGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCGACGCCGCCCGGCACCGGTGTGATCCAGCCGGCGCgggg              >  1:164408/1‑147 (MQ=255)
                                                                                                                                     aaaCGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCGACGCCGCCCGGCACCGGTGTGATCCAGCCGGCGCGgggcag           <  6:45053/149‑1 (MQ=255)
                                                                                                                                      aaCGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCGGGCACGCCCGGGTCATCGGCCCGACGCCGCCCGGCACCGGTGTGATCCAGCCGGCGCGgggcagg          <  4:330673/149‑1 (MQ=255)
                                                                                                                                       aCGGCACCTTTGCAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCGACGCCGCCCGGCACCGGTGTGATCCAGCCGGCGCGgggcaggg         >  8:143208/1‑149 (MQ=255)
                                                                                                                                       aCGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGTACGCCCGGGTCATCGGCCCGACGCCGCCCGGCACCGGTGTGATCCAGCCGGCGCGgggcaggg         >  1:260423/1‑149 (MQ=255)
                                                                                                                                              cTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCGACGCCGCCCGGCACCGGTGTGATCCAGCCGGCGCGGGGCAGGGCAGtctc  <  5:86785/149‑1 (MQ=255)
                                                                                                                                               tGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCGACGCCGCCCGGCACCGGTGTGATCCAGCCGGCGCGgggcagg          >  2:169015/1‑140 (MQ=255)
                                                                                                                                               tGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCGACGCCGCCCGGCACCGGTGTGATCCAGCCGGCGCGgggcagg          <  1:169015/140‑1 (MQ=255)
                                                                                                                                                 cAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCGACGCCGCCCGGCACCGGTGTGATCCAGCCGgcgc                  <  1:3218/130‑1 (MQ=255)
                                                                                                                                                 cAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCGACGCCGCCCGGCACCGGTGTGATCCAGCCGgcgc                  >  2:3218/1‑130 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GTCAGCTACAACGAACTGCTGAAGTCCTGGGAGGCGGCCAGGAAGTGATCGGCTGAACGCCGGGGTGGTTTCACCGACCGCGGGGGCCGTGCAACCGCTCGCTCGCGAAGTGGCCACTAGTCACGGGCAAAAAAAACGGCACCTGCAAAGGTGCCGTTTTCGCATCACGCAGTATTACTTGTGCAGCTCTTCTGCCGCATACAGCGTGTTCTCCAGCAGGCACGCCCGGGTCATCGGCCCGACGCCGCCCGGCACCGGTGTGATCCAGCCGGCGCGGGGCAGGGCAGTCTC  >  NC_002947/2584466‑2584756

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: