Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,963,305 C→A 5.6% G117V (GGC→GTC)  fleS ← two‑component system sensor histidine kinase FleS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,963,3050CA 5.6% 229.0 / 2.4 71G117V (GGC→GTC) fleStwo‑component system sensor histidine kinase FleS
Reads supporting (aligned to +/- strand):  ref base C (34/33);  new base A (2/2);  total (36/35)
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.21e-01

CAGGTCGTTGAGCAGCACCAATTGGCCAGGCTCGGCATCCAGTGAGCGGGTGGCGATCGACAGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCACCCCGCCAGGCAGCAAGTCCAGCAGGTTCTGCAGGCGGTTGGCCAGGCGTTCCTTTTCACTCAGTTCGTC  >  NC_002947/4963157‑4963453
                                                                                                                                                    |                                                                                                                                                    
ctggTCGTGGAGCAGCACCACTTGGCCAGGCTCGGCATCCAGGGAGCGGGTGGCGATCGACAGGCGGCGACCGTCGCGCAGCGAGACTTCGTGTCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGc                                                                                                                                                      <  5:380890/147‑1 (MQ=255)
 agGTCGTTGAGCAGCACCAATTGGCCAGGCTCGGCATCCAGTGAGCGGGTGGCGATCGACAGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGcc                                                                                                                                                     >  2:21091/1‑149 (MQ=255)
 agGTCGTTGAGCAGCACCAATTGGCCAGGCTCGGCATCCAGTGAGCGGGTGGCGATCGACAGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGcc                                                                                                                                                     >  2:190742/1‑149 (MQ=255)
  ggTCGTTGAGCAGCACCAATTGGCCAGGCTCGGCATCCAGTGAGCGGGTGGCGATCGACAGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCg                                                                                                                                                    >  5:310041/1‑149 (MQ=255)
   gTCGTTGAGCAGCACCAATTGGCCAGGCTCGGCATCCAGTGAGCGGGTGGCGATCGACAGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGa                                                                                                                                                   >  7:332190/1‑149 (MQ=255)
    tCGTTGAGCAGCACCAATTGGCCAGGCTCGGCATCCAGTGAGCGGGTGGCGATCGACAGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGAt                                                                                                                                                  <  1:21091/149‑1 (MQ=255)
        tGAGCAGCACCAATTGGCCAGGCTCGGCATCCAGTGAGCGGGTGGCGATCGACAGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCag                                                                                                                                               <  2:336790/148‑1 (MQ=255)
              gcaCCAATTGGCCAGGCTCGGCATCCAGTGAGCGGGTGGCGATCGACAGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAgg                                                                                                                                            <  6:254602/145‑1 (MQ=255)
              gcaCCAATTGGCCAGGCTCGGCATCCAGTGAGCGGGTGGCGATCGACAGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAgg                                                                                                                                            >  5:254602/1‑145 (MQ=255)
               caCCAATTGGCCAGGCTCGGCATCCAGTGAGCGGGTGGCGATCGACAGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGACGATCAGAGGCTCAc                                                                                                                                       >  7:43002/1‑149 (MQ=255)
               caCCAATTGGCCAGGCTCGGCATCCAGTGAGCGGGTGGCGATCGACAGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGACGATCAGAGGCTCAc                                                                                                                                       >  7:37091/1‑149 (MQ=255)
                aCCAATTGGCCAGGCTCGGCATCCAGTGAGCGGGTGGCGATCGACAGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCAcc                                                                                                                                      <  1:190742/149‑1 (MQ=255)
                aCCAATTGGCCAGGCTCGGCATCCAGTGAGCGGGTGGCGATCAACAGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCAcc                                                                                                                                      >  6:340961/1‑149 (MQ=255)
                  cAATTGGCCAGGCTCGGCATCCAGTGAGCGGGTGGCGATCGACAGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGa                                                                                                                                    <  8:112180/149‑1 (MQ=255)
                  cAATTGGCCAGGCTCGGCATCCAGTGAGCGGGTGGCGATCGACAGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGa                                                                                                                                    <  8:332190/149‑1 (MQ=255)
                        gCCAGGCTCGGCATCCAGTGAGCGGGTGGCGATCGACAGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACt                                                                                                                              >  7:290430/1‑149 (MQ=255)
                         ccAGGCTCGGCATCCTGTGAGCGGGTGGCGATCGACAGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGTGCAAAGCTACGGGCGATCACCTGGCGGCACAATTGGCCGATCAGAGGCTCACCGAGCAACTc                                                                                                                             <  7:21872/149‑1 (MQ=255)
                               tCGGCATCCAGTGAGCGGGTGGCGATCGACAGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAggc                                                                                                                       >  6:286267/1‑149 (MQ=255)
                                  gCATCCAGTGAGCGGGTGGCGATCGACAGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAggcggc                                                                                                                    >  4:174563/1‑149 (MQ=255)
                                        agggatcGGGGGGCGATCGACAGGCGGGGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGTGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGtt                                                                                                              <  8:321462/143‑1 (MQ=255)
                                        aGTGAGCGGGTGGCGATCGACAGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGtt                                                                                                              <  8:59457/149‑1 (MQ=255)
                                         gTGAGCGGGTGGCGATCGACAGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACAGCGCAACTCGCAGGCGGCAGGGTTg                                                                                                             >  7:137568/1‑149 (MQ=255)
                                            aGCGGGTGGCGATCGACAGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGc                                                                                                           <  5:86436/148‑1 (MQ=255)
                                            aGCGGGTGGCGATCGACAGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCt                                                                                                          >  1:223864/1‑149 (MQ=255)
                                                      gATCGACAGGCGGCGACCGTCGCGCAGCGAGACTTCGTGTCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCAc                                                                                                  >  1:190812/1‑147 (MQ=255)
                                                          gACAGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGcc                                                                                            <  4:340284/149‑1 (MQ=255)
                                                           aCAGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGcccc                                                                                           >  7:62215/1‑149 (MQ=255)
                                                           aCAGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGccc                                                                                           <  3:114026/149‑1 (MQ=255)
                                                             aGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTg                                                                                                             >  5:179737/1‑129 (MQ=255)
                                                             aGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTg                                                                                                             <  6:179737/129‑1 (MQ=255)
                                                             aGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCt                                                                                          >  2:168193/1‑148 (MQ=255)
                                                              ggcggcGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGa                                                                                        <  1:36071/149‑1 (MQ=255)
                                                                 ggcgACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCa                                                                                     >  7:15393/1‑149 (MQ=255)
                                                                    gACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGCTCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATc                                                                                   >  7:138709/1‑148 (MQ=255)
                                                                    gACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCCCCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATc                                                                                   >  7:86442/1‑148 (MQ=255)
                                                                    gACCGGCGCGCAGCGAGACTTCGTGGCCGTCGGCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATcg                                                                                  <  2:354288/149‑1 (MQ=255)
                                                                     aCCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATcga                                                                                 >  7:172881/1‑149 (MQ=255)
                                                                       cGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCggg                                                                                                                >  6:14103/1‑116 (MQ=255)
                                                                       cGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCggg                                                                                                                <  5:14103/116‑1 (MQ=255)
                                                                           gcgcAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATcgatcacga                                                                           <  5:286267/149‑1 (MQ=255)
                                                                           gcgcAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATcgatcacga                                                                           >  2:6458/1‑149 (MQ=255)
                                                                            cgcAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATcgatcacgat                                                                          <  2:101233/149‑1 (MQ=255)
                                                                                gCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCAcc                                                                      >  7:327210/1‑149 (MQ=255)
                                                                                      cTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCACCCCGCCa                                                                <  7:269254/149‑1 (MQ=255)
                                                                                       ttCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCAcc                                                                                                 >  5:67892/1‑115 (MQ=255)
                                                                                       ttCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCAcc                                                                                                 <  6:67892/115‑1 (MQ=255)
                                                                                              ccGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCACCCCGCCAGGCAGCaa                                                        <  8:15393/149‑1 (MQ=255)
                                                                                              ccGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCACCCCGCCAGGCAGCaa                                                        <  1:6458/149‑1 (MQ=255)
                                                                                              ccGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCACCCCGCCAGGCAGCaa                                                        <  1:168193/149‑1 (MQ=255)
                                                                                                       ctGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCACCCCGCCAGGCAGCAAGTCcagcag                                               <  7:270566/148‑1 (MQ=255)
                                                                                                          cgcgGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCACCCCGCCAGGCAGCAAGTCCAGCAGGtt                                            >  8:249208/1‑149 (MQ=255)
                                                                                                              gggCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCACCCCGCCAGGCAGCAAGTCCAGCAGGTTCTGc                                        <  3:174563/149‑1 (MQ=255)
                                                                                                                  aaaGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCACCCCGCCAGGCAGCAAGTCCAGCAGGTTCTGCAGGc                                    >  7:335788/1‑149 (MQ=255)
                                                                                                                  aaaGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCACCCCGCCAGGCAGCAAGTCCAGCAGGTTCTGCAGGc                                    >  2:315707/1‑149 (MQ=255)
                                                                                                                        aCGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCACCCCGCCAGGCAGCAAGTCCCGCAGGTTCTGCAGGCGGTTgg                              >  7:127928/1‑149 (MQ=255)
                                                                                                                         cGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCACCCCGCCAGGCAGCAAGTCCAGCAGGTTCTGCAGGCGGTTGGc                             >  3:28931/1‑149 (MQ=255)
                                                                                                                         cGGGCGATCACCTGGCGCCACAATTGGACGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCACCCCGCCAGGCAGCAAGTCCAGCAGGTTCTGCAGGCGGTTg                               <  8:43002/147‑1 (MQ=255)
                                                                                                                         cGGGCGATCACCTGGCGCCACAATTGGACGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCACCCCGCCAGGCAGCAAGTCCAGCAGGTTCTGCAGGCGGTTg                               <  8:37091/147‑1 (MQ=255)
                                                                                                                          gggCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCACCCCGCCAGGCAGCAAGTCcagc                                                 <  3:45150/128‑1 (MQ=255)
                                                                                                                          gggCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCACCCCGCCAGGCAGCAAGTCcagc                                                 >  4:45150/1‑128 (MQ=255)
                                                                                                                            gCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCACCCCGCCAGGCAGCAAGTCCAGCAGGTTCTGCAGGCGGTTGGCCAg                          >  2:234749/1‑149 (MQ=255)
                                                                                                                             cGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCACCCCGCCAGGCAGCAAGTCCAGCAGGTTCTGCAGGCGGTTGGCCAgg                         >  4:107042/1‑149 (MQ=255)
                                                                                                                              gATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCACCCCGCCAGGCAGCAAGTCCAGCAGGTTCTGCAGGCGGTTGGCCAGGc                        <  8:327210/149‑1 (MQ=255)
                                                                                                                               aTCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCACCCCGCCAGGCAGCAAGTCCAGCAGGTTCTGCAGGCGGTTGGCCAGGCg                       >  5:204145/1‑149 (MQ=255)
                                                                                                                                tCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCACCCCGCCAAGCAGCAAGTCCAGCAGGTTCTg                                         <  7:232829/130‑1 (MQ=255)
                                                                                                                                tCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCAGACCATCGATAACGAACACCCCGCCAAGCAGCAACTCCAGCAGGTTATg                                         >  8:232829/1‑130 (MQ=255)
                                                                                                                                   ccTGGCGCCACAATTGGCGGATCAGAGGCTCACCGAGCAACTGGCAGGGGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGTTCACCGCGCCAGGCAGCAAGTCCAGCAGGTTCTGCAGGCGGTTGGCCAGGCGTTcc                   <  8:86442/149‑1 (MQ=255)
                                                                                                                                   ccTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCACCCCGCCAGGCAGCAAGTCCAGCAGGTTCTGCAGGCGGTTGGCCAGGCGTTcc                   >  7:295083/1‑149 (MQ=255)
                                                                                                                                       gcgcCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCACCCCGCCAGGCAGCAAGTCCAGCAGGTTCTGCAGGCGGTTGGCCAGGCGTTCCttt                <  7:51664/148‑1 (MQ=255)
                                                                                                                                         gcCACAATTGGCCGATCAGAGGCTCCCCGAGCAACTCGCAGGCGGCCGGGTTGGCGTCGCGCACCAGGCCCTGCCCATCGATCACGATCACCCCGCCAGGCAGCAAGTCCAGCAGGTTCGGCAGGCGGGGGGCCAGGCGTTCCTTTTcc             >  6:299143/1‑148 (MQ=255)
                                                                                                                                         gcCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCACCCCGCCAGGCAGCAAGTCCAGCAGGTTCTGCAGGCGGTTGGCCAGGCGTTCCTTTTCa             <  8:138709/149‑1 (MQ=255)
                                                                                                                                                   gCCGATCAGAGGCTCTCCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCACCCCGCCAGGCAGCAAGTCCAGCAGGTTCTGCAGGCGGTTGGCCAGGCGTTCCTTTTCACTCAGTtcg    <  8:172881/148‑1 (MQ=255)
                                                                                                                                                   gCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCACCCCGCCAGGCAGCAAGTCCAGCAGGTTCTGCAGGCGGTTGGCCAGGCGTTCCTTTTCACTCAGTtcgt   >  4:162629/1‑149 (MQ=255)
                                                                                                                                                   gCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCACCCCGCCAGGCAGCAAGTCCAGCAGGTTCTGCAGGCGGTTGGCCAGGCGTTCCTTTTCACTCAGTtcgt   <  2:223864/149‑1 (MQ=255)
                                                                                                                                                    ccGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCACCCCGCCAGGCAGCAAGTCCAGCAGGTTCTGCAGGCGGTTGGCCAGGCGTTCCTTTTCACTCAGTtcgtc  <  6:204145/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAGGTCGTTGAGCAGCACCAATTGGCCAGGCTCGGCATCCAGTGAGCGGGTGGCGATCGACAGGCGGCGACCGTCGCGCAGCGAGACTTCGTGGCCGTCGTCCTTGCGCGGGGCAAAGCTACGGGCGATCACCTGGCGCCACAATTGGCCGATCAGAGGCTCACCGAGCAACTCGCAGGCGGCCGGGTTGGCTTCGCGCACCAGGCCCTGACCATCGATCACGATCACCCCGCCAGGCAGCAAGTCCAGCAGGTTCTGCAGGCGGTTGGCCAGGCGTTCCTTTTCACTCAGTTCGTC  >  NC_002947/4963157‑4963453

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