Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,488,151 G→A 92.0% E494K (GAG→AAG)  PP_4824 → two‑component system sensor histidine kinase/response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,488,1510GA92.0% 215.1 / 7.3 75E494K (GAG→AAG) PP_4824two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base G (3/3);  new base A (35/34);  total (38/37)
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.79e-01

CTGCTGACGCTGATCAACGAAATTCTGGACATTTCCAAGCTGGAATCCCGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCGAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTGAAAAACACCGACCAGGGCGAGATCCTGCTGGTGGTGGCGCTGGAT  >  NC_002947/5488007‑5488294
                                                                                                                                                |                                                                                                                                                
ctgctgACGCTGATCAACGAAATTCTGGACATTTCCAAGCTGGAATCCCGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCt                                                                                                                                              >  3:277551/1‑149 (MQ=255)
   ctgACGCTGATCAACGAAATTCTGGACATTTCCAAGCTGGAATCCCGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGa                                                                                                                                            >  3:67559/1‑148 (MQ=255)
   ctgACGCTGATCAACGAAATTCTGGACATTTCCAAGCTGGAATCCCGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGAt                                                                                                                                           <  7:211835/149‑1 (MQ=255)
      aCGCTGATCAACGAAATTCTGGACATTTCCAAGCTGGAATCCCGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCTAGCAACAGAACATCAAGCTGATCa                                                                                                                                         <  7:252465/148‑1 (MQ=255)
      aCGCTGATCAACGAAATTCTGGACATTTCCAAGCTGGAATCCCGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAg                                                                                                                                                >  2:74297/1‑141 (MQ=255)
      aCGCTGATCAACGAAATTCTGGACATTTCCAAGCTGGAATCCCGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCa                                                                                                                                         <  4:67559/148‑1 (MQ=255)
      aCGCTGATCAACGAAATTCTGGACATTTCCAAGCTGGAATCCCGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAAGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAg                                                                                                                                                <  1:74297/141‑1 (MQ=255)
           gATCAACGAAATTCTGGACATTTCCAAGCTGGAATCCCGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCa                                                                                                                                   >  5:18525/1‑149 (MQ=255)
           gATCAACGAAATTCTGGACATTTCCAAGCTGGAATACCGGCAGATCGAGCTGGATGACGTGCAGTTCGAACTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCa                                                                                                                                   >  5:14192/1‑149 (MQ=255)
            aTCAACGAAATTCTGGACATTTCCAAGCTGGAATCCCGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTc                                                                                                                                    <  6:250980/147‑1 (MQ=255)
               aaCGAAATTCTGGACATTTCCAAGCTGGAATCCCGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCAcc                                                                                                                                 >  6:234883/1‑147 (MQ=255)
               aaCGAAATTCTGGACATTTCCAAGCTGGAATCCCGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCa                                                                                                                               >  4:297758/1‑149 (MQ=255)
                aCGAAATTCTGGACATTTCCAAGCTGGAATCCCGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCGAGCTGATCAGCTTCAccccg                                                                                                                              >  8:158623/1‑147 (MQ=255)
                aCGAAATTCTGGACATTTCCAAGCTGGAATCCCGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCGAGCTGATCAGCTTCACCCAg                                                                                                                              >  2:202082/1‑149 (MQ=255)
                 cGAAATTCTGGACATTTCCAAGCTGGAATCCCGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGc                                                                                                                                               >  5:58151/1‑131 (MQ=255)
                 cGAAATTCTGGACATTTCCAAGCTGGAATCCCGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGc                                                                                                                                               <  6:58151/131‑1 (MQ=255)
                  gAAATTCTGGACATTTCCAAGCTGGAATCCCGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGc                                                                                                                             >  4:274579/1‑148 (MQ=255)
                    aaTTCTGGACATTTCCAAGCTGGAATCCCGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCcat                                                                                                                          <  3:274579/149‑2 (MQ=255)
                      ttCTGGACATTTCCAAGCTGGAATCCCGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACaa                                                                                                                        >  8:291753/1‑149 (MQ=255)
                             cATTTCCAAGCTGGAATCCCGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACcgc                                                                                                                 <  4:277551/149‑1 (MQ=255)
                                  ccAAGCTGGAATCCCGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGt                                                                                                             >  6:14649/1‑148 (MQ=255)
                                  ccAAGCTGGAATCCCGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGt                                                                                                             <  3:297758/148‑1 (MQ=255)
                                   cAAGCTGGAATCCCGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGtca                                                                                                           >  4:106879/1‑149 (MQ=255)
                                      gCTGGAATCCCGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGtcatca                                                                                                        <  5:14649/149‑1 (MQ=255)
                                           aaTCCCGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGga                                                                                                   <  6:52192/149‑1 (MQ=255)
                                               ccGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGga                                                                                                   >  5:73143/1‑145 (MQ=255)
                                               ccGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGga                                                                                                   <  6:73143/145‑1 (MQ=255)
                                                cGGCAGATCTAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAAAAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGAccc                                                                                              <  6:262289/149‑1 (MQ=255)
                                                cGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGAcc                                                                                               <  4:212367/148‑1 (MQ=255)
                                                cGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAAAGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGAcc                                                                                               >  3:212367/1‑148 (MQ=255)
                                                  gCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGgaga                                                                                                 >  1:196462/1‑144 (MQ=255)
                                                  gCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGgaga                                                                                                 <  2:196462/144‑1 (MQ=255)
                                                       tCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCgg                                                                                                               <  4:22656/125‑1 (MQ=255)
                                                       tCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCgg                                                                                                               >  3:22656/1‑125 (MQ=255)
                                                               gATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTgcgc                                                                                 <  2:278414/147‑1 (MQ=255)
                                                                    cGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCgg                                                                                                    <  1:309969/123‑1 (MQ=255)
                                                                    cGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCgg                                                                                                    >  2:309969/1‑123 (MQ=255)
                                                                     gTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCtt                                                                         >  2:155790/1‑149 (MQ=255)
                                                                      tGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTg                                                                        >  6:31892/1‑149 (MQ=255)
                                                                       gCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGt                                                                       >  4:252567/1‑149 (MQ=255)
                                                                        cAGTTCGAACTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTc                                                                      <  6:14192/149‑1 (MQ=255)
                                                                         aGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTc                                                                      >  8:180456/1‑148 (MQ=255)
                                                                         aGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTc                                                                      <  6:18525/148‑1 (MQ=255)
                                                                              gACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCct                                                                >  2:70042/1‑149 (MQ=255)
                                                                                 cTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCGAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCctgct                                                             >  2:136862/1‑149 (MQ=255)
                                                                                    aaCGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCctgctg                                                            <  1:124711/147‑1 (MQ=255)
                                                                                     aCGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGACATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGa                                                          >  5:113237/1‑148 (MQ=255)
                                                                                       gcgcTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGaaaa                                                       <  5:31892/149‑1 (MQ=255)
                                                                                             aTCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGccc                                                  <  1:155790/148‑1 (MQ=255)
                                                                                             aTCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGccc                                                  <  5:234883/148‑1 (MQ=255)
                                                                                             aTCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGACATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCt                                                 <  6:113237/149‑1 (MQ=255)
                                                                                                     ttGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGTGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTGAAAAaca                                         <  1:61994/149‑1 (MQ=255)
                                                                                                        ccagAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGAc                                                                                           >  7:263664/4‑96 (MQ=255)
                                                                                                        ccTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGAc                                                                                           <  8:263664/96‑1 (MQ=255)
                                                                                                        ccTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTGAAAAACAcc                                       <  7:94011/148‑1 (MQ=255)
                                                                                                         cTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTGAAAAACACCGa                                    >  5:244426/1‑149 (MQ=255)
                                                                                                                tCTTCCGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTGAAAAACACCGACCAGGGc                             >  5:30238/1‑149 (MQ=255)
                                                                                                                    ccGCGCCAAGGCCGAGCAACAGAACATCGAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTGAAAAACACCGACCAGGGCg                            <  1:136862/146‑1 (MQ=255)
                                                                                                                    ccGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTGAAAAACACCGACCAGGGc                             <  7:291753/145‑1 (MQ=255)
                                                                                                                    ccGCGCCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTGAAAAACACCGACCAGGGc                             >  5:251723/1‑145 (MQ=255)
                                                                                                                      gcgcCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAACCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGcc                                                   <  3:112704/122‑1 (MQ=255)
                                                                                                                      gcgcCAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAACCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGcc                                                   >  4:112704/1‑122 (MQ=255)
                                                                                                                         ccAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTGAAAAACACCGACCAGGGCgag                          <  6:281756/143‑1 (MQ=255)
                                                                                                                         ccAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTGAAAAACACCGACCAGGGCgag                          >  5:281756/1‑143 (MQ=255)
                                                                                                                         ccAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCATGTCGAGCCTGCTGGAAAACGCCCTGAAAAACACCGACCAGGGCgag                          >  5:276797/1‑143 (MQ=255)
                                                                                                                          cAAGGCCGAGCATCAGAACATCAAGCTGATGAGCTTCACACAGCCAGAAGTACAGCGGGTCATCAGCGGAGAACCTACTCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGTAAAACGCCCTGAAAAACACGGACCAGGGCGAGATCctgc                   <  6:251723/149‑1 (MQ=255)
                                                                                                                          cAAGGCCGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTGAAAAACACCGACCAGGGCGAGATCctgc                   >  6:308949/1‑149 (MQ=255)
                                                                                                                               ccGAGCAACAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTGAAAAACACCGACCAGGGCGAGATCCTGCtggtg              >  1:93889/1‑149 (MQ=255)
                                                                                                                                      aCAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTGAAAAACACCGACCAGGGCGAGATCCTGCTGGTGGTGgcgc       >  5:27803/1‑149 (MQ=255)
                                                                                                                                       cAGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTGAAAAACACCGACCAGGGCGAGATCCTGCTGGTGGTGGCGCt      <  6:244426/149‑1 (MQ=255)
                                                                                                                                        aGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTGAAAAACACCGACCAGGGCGAGATCCTGCTGGTGGTGGCGCTg     >  2:292327/1‑149 (MQ=255)
                                                                                                                                        aGAACATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTGAAAAACACCGACCAGGGCGAGATCCTGCTGGTGGTGGCGCTg     <  2:93889/149‑1 (MQ=255)
                                                                                                                                           aCATCGAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTGAAAAACACCGACCAGGGCGAGATCCTGCTGGTGGTGGCGCTGGAt  <  7:158623/149‑1 (MQ=255)
                                                                                                                                            cATCGAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTGAAAAACACCGAACCAGGGCGAGATCCTGCTGGTGGTGGCGCTGGAt  <  1:202082/149‑1 (MQ=255)
                                                                                                                                            cATCAAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTGAAAAACACCGACCAGGGCGAGATCCTGCTGGTGGTGGCGCTGGAt  >  3:169836/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CTGCTGACGCTGATCAACGAAATTCTGGACATTTCCAAGCTGGAATCCCGGCAGATCGAGCTGGATGACGTGCAGTTCGACCTTAACGCGCTGATCGAGGATTGCCTGAACATCTTCCGCGCCAAGGCCGAGCAACAGAACATCGAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTGAAAAACACCGACCAGGGCGAGATCCTGCTGGTGGTGGCGCTGGAT  >  NC_002947/5488007‑5488294

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: