Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,846,286 T→A 5.1% W193R (TGG→AGG)  pfeS‑II → histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,846,2860TA 5.1% 239.9 / 2.5 79W193R (TGG→AGG) pfeS‑IIhistidine kinase
Reads supporting (aligned to +/- strand):  ref base T (39/36);  new base A (2/2);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

GCAGCGCTTCAATCCGGGCCAGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGTC  >  NC_002947/1846138‑1846434
                                                                                                                                                     |                                                                                                                                                     
gCAGCGCTTCAATCCGGGCCAGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCt                                                                                                                                                       >  2:224846/1‑149 (MQ=255)
gCAGCGCTTCAATCCGGGCCAGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCt                                                                                                                                                       >  2:219496/1‑149 (MQ=255)
gCAGCGCTTCAATCCGGGCCAGCACCGCTTGCTCTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTAGCTGTATTGTGTAGGCCTGTACCGCATGCTGATCGTGCGGCTGACCCAACTTCGCGAGCAGGCCAATGCCt                                                                                                                                                       >  8:175594/1‑149 (MQ=255)
   gcgcTTCAATCCGGGCCAGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCTTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGc                                                                                                                                                    >  4:204415/1‑149 (MQ=255)
   gcgcTTCAATCCGGGCCAGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGc                                                                                                                                                    <  1:219496/149‑1 (MQ=255)
   gcgcTTCAATCCGGGCCAGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGc                                                                                                                                                    <  1:224846/149‑1 (MQ=255)
       ttCAATCCGGGCCAGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGc                                                                                                                                                >  1:275886/1‑149 (MQ=255)
          aaTCCGGGCCAGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGcc                                                                                                                                               >  5:137794/1‑147 (MQ=255)
          aaTCCGGGCCAGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGcc                                                                                                                                               <  6:137794/147‑1 (MQ=255)
                     gCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCgg                                                                                                                                 >  6:76822/1‑149 (MQ=255)
                        ccGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCTTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGccc                                                                                                                              <  3:204415/149‑1 (MQ=255)
                         cGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCg                                                                                                                             >  8:128827/1‑149 (MQ=255)
                                tgtgGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGAc                                                                                                                      >  3:183776/1‑149 (MQ=255)
                                tgtgGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGAc                                                                                                                      >  7:168122/1‑149 (MQ=255)
                                  tgGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTc                                                                                                                    <  5:76822/149‑1 (MQ=255)
                                    gcgcATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTcgc                                                                                                                  <  8:168122/149‑1 (MQ=255)
                                    gcgcATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTcgc                                                                                                                  <  7:128827/149‑1 (MQ=255)
                                        aTTGTCACCAACGGCATCATCCCGGGGCATGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCa                                                                                                               >  1:22388/1‑149 (MQ=255)
                                            tCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTg                                                                                                                         <  4:194885/134‑1 (MQ=255)
                                            tCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTg                                                                                                                         >  3:194885/1‑134 (MQ=255)
                                            tCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCa                                                                                                               >  7:228344/1‑144 (MQ=255)
                                            tCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCa                                                                                                               <  8:228344/144‑1 (MQ=255)
                                            tCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATc                                                                                                          >  4:83928/1‑149 (MQ=255)
                                              aCCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGc                                                                                                        >  8:176755/1‑149 (MQ=255)
                                                cAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGccc                                                                                                      >  3:53550/1‑149 (MQ=255)
                                                 aaCGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCAgg                                                                                                                                                     >  6:123838/1‑102 (MQ=255)
                                                 aaCGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCAgg                                                                                                                                                     <  5:123838/102‑1 (MQ=255)
                                                      catcatCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCAt                                                                                                 >  5:89575/1‑148 (MQ=255)
                                                       atcatcCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGa                                                                                               <  4:53550/149‑1 (MQ=255)
                                                        tcatcCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGAc                                                                                              <  4:139952/149‑1 (MQ=255)
                                                        tcatcCCGGGGCGGTTCACGTTGCGCTTCTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGAc                                                                                              <  4:134379/149‑1 (MQ=255)
                                                         catcCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGAc                                                                                              >  1:51786/1‑148 (MQ=255)
                                                             ccGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAAt                                                                                          <  2:275886/148‑1 (MQ=255)
                                                               ggggCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTgg                                                                                       >  2:197447/1‑149 (MQ=255)
                                                                      gTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCAGGCGGGCCGACAGCTGATCGGCCCGc                                                                                                                            >  3:274782/1‑106 (MQ=38)
                                                                      gTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCAGGCGGGCCGACAGCTGATCGGCCCGc                                                                                                                            <  4:274782/106‑1 (MQ=38)
                                                                            cTTGCTTTTGTGTGTAGGCCTGTAGCGAATGCTGATCGTGCCGCTGAACCAACTTCGCTAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCGCCATCGCCCTGCATGACGAATTGGGCGAACTGgc                                                                             <  2:233735/147‑1 (MQ=255)
                                                                            cTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCGGGCATGACGAATTGGGCGAACt                                                                                <  7:57455/144‑1 (MQ=255)
                                                                            cTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACt                                                                                >  8:57455/1‑144 (MQ=255)
                                                                            cTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTACGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCCCGCGCCATCGCCAGGCATGACGAATTGGGCGAACTGgc                                                                             >  1:233735/1‑147 (MQ=255)
                                                                              tGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGt                                                                         >  5:98362/1‑149 (MQ=255)
                                                                                 tgttgtGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTg                                                                        >  6:28477/1‑147 (MQ=255)
                                                                                 tgttgtGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTg                                                                        <  5:28477/147‑1 (MQ=255)
                                                                                 tgttgtGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGc                                                                       >  6:314025/1‑148 (MQ=255)
                                                                                  gttgtGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGTAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCt                                                                     >  6:373755/1‑149 (MQ=255)
                                                                                   ttgtGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGc                                                                       <  5:314025/146‑1 (MQ=255)
                                                                                     gtgtgtAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGTAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCttt                                                                   <  5:373755/148‑1 (MQ=255)
                                                                                      tgtgtAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGa                                                                 >  4:71027/1‑149 (MQ=255)
                                                                                        tgtAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGgcg                                                                            >  1:87743/1‑136 (MQ=255)
                                                                                        tgtAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGgcg                                                                            <  2:87743/136‑1 (MQ=255)
                                                                                        tgtAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGAc                                                                <  1:197447/148‑1 (MQ=255)
                                                                                          tAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTgg                                                                                       >  5:310860/1‑123 (MQ=255)
                                                                                          tAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTgg                                                                                       <  6:310860/123‑1 (MQ=255)
                                                                                          tAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCa                                                              <  3:83928/148‑1 (MQ=255)
                                                                                               cTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGCCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGc                                                        >  2:59283/1‑149 (MQ=255)
                                                                                               cTGAACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGc                                                        <  4:310601/149‑1 (MQ=255)
                                                                                                 gTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATg                                                          >  6:109258/1‑145 (MQ=255)
                                                                                                 gTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATg                                                          <  5:109258/145‑1 (MQ=255)
                                                                                                    ccGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGc                                                   <  3:71027/149‑1 (MQ=255)
                                                                                                        aTGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCt                                               <  2:22388/149‑1 (MQ=255)
                                                                                                          gCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTg                                              >  6:322988/1‑148 (MQ=255)
                                                                                                            tGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGACGAGAGGATGc                                             >  4:209577/1‑147 (MQ=255)
                                                                                                             gATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAg                                          >  5:88601/1‑149 (MQ=255)
                                                                                                              aTCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCaa                                           <  4:183776/147‑1 (MQ=255)
                                                                                                              aTCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCaa                                           <  7:176755/147‑1 (MQ=255)
                                                                                                                  tGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGCCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACg                                     <  1:59283/149‑1 (MQ=255)
                                                                                                                     cGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTg                                  <  1:238616/149‑1 (MQ=255)
                                                                                                                       cTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCAc                                      >  5:270828/1‑143 (MQ=255)
                                                                                                                       cTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCAc                                      <  6:270828/143‑1 (MQ=255)
                                                                                                                         gAACCAACTTCGCGAGCAGGCCAATGCCTCGCGGGCCGATCCGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTa                              >  5:186450/1‑149 (MQ=255)
                                                                                                                          aaCCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTAt                             >  4:265949/1‑149 (MQ=255)
                                                                                                                            ccAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGAc                                                                <  2:145183/112‑1 (MQ=255)
                                                                                                                            ccAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGAc                                                                >  1:145183/1‑112 (MQ=255)
                                                                                                                            ccAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATg                            >  2:347186/1‑148 (MQ=255)
                                                                                                                              aaCTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGCCGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATgcag                         >  6:281485/1‑149 (MQ=255)
                                                                                                                                 ttCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATg                            <  7:18531/143‑1 (MQ=255)
                                                                                                                                 ttCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATg                            >  8:18531/1‑143 (MQ=255)
                                                                                                                                        gCAGGCCAATGCCTGGCGGGCCGATCCGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTtgc               >  4:299968/1‑149 (MQ=255)
                                                                                                                                         cAGGCCAATGCCTGGCGGGCCGATCCGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTtgct              <  3:299968/149‑1 (MQ=255)
                                                                                                                                          aGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTtgctg             <  5:322988/149‑1 (MQ=255)
                                                                                                                                            gCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCGGGACTCGCGCACCAACGCCGGGCATGACGAATTGGGCGAACTGGCGCGTGCATTTGACCAGATGGCCGAGCGGCGGCAAGGCACGGTGGCTATGCAGCAACtgctgctgca           <  7:175594/149‑2 (MQ=255)
                                                                                                                                                  gCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCCGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGCCTCCAAGGCACGATTGCTATGCAGCAGCAGTTGCATCGCGACCt     >  8:346779/1‑149 (MQ=255)
                                                                                                                                                     tGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGTc  >  4:11269/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAGCGCTTCAATCCGGGCCAGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGTC  >  NC_002947/1846138‑1846434

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: