Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,166,531 A→G 5.9% L450P (CTG→CCG) ‡ hrpA ← ATP‑dependent helicase HrpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,166,5310AG 5.9% 247.1 / 4.5 84L450P (CTG→CCG) ‡hrpAATP‑dependent helicase HrpA
Reads supporting (aligned to +/- strand):  ref base A (41/38);  new base G (2/3);  total (43/41)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

ATCGGCAAGCGCGCCAACTGGCGGCCAATTGGCGTCAGCTGGTTCTCGCGGTTAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGTTGAAGTCCTCTTCGCTGTAGAGACGCACGCAAATGCCCGGCTCGACCCGGCCGCAACGGCCTTTA  >  NC_002947/5166387‑5166677
                                                                                                                                                  |                                                                                                                                                  
aTCGGCAAGCGCGCCAACTGGCGGCCAATTGGCGTCAGCTGGTTCTCGCGGTTAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGgcg                                                                                                                                                <  1:373274/149‑1 (MQ=255)
  cGGCAAGCGCGCCAACTGGCGGCCAATTGGCGTCAGCTGGTTCTCGCGGTTAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCTGGgcgc                                                                                                                                               >  4:69065/1‑148 (MQ=255)
  cGGCAAGCGCGCCAACTGGCGGCCAATTGGCGTCAGCTGGTTCTCGCGGTTAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGgcg                                                                                                                                                >  1:8644/1‑147 (MQ=255)
   ggCAAGCGCGCCAACTGGCGGCCAATTGGCGTCAGCTGGTTCTCGCGGTTAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGgcgc                                                                                                                                               >  1:7747/1‑147 (MQ=255)
   ggCAAGCGCGCCAACTGGCGGCCAATTGGCGTCAGCTGGTTCTCGCGGTTAACCGCCGAAAGCTC‑‑CTGCAACAGGTCGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGgcgc                                                                                                                                               <  2:7747/147‑1 (MQ=255)
    gCAAGCGCGCCAACTGGCGGCCAATTGGCGTCAGCTGGTTCTCGCGGTTAACCGCCGAAAGCTCGACTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGgcg                                                                                                                                                >  1:117853/1‑147 (MQ=255)
        gcgcgcCAACTGGCGGCCAATTGGCGTCAGCTGGTTCTCGCGGTTAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCTGGGCGCaatcgat                                                                                                                                        <  3:69065/149‑7 (MQ=255)
        gcgcgcCAACTGGCGGCCAATTGGCGTCAGCTGGTTCTCGCGGTTAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCAGCTCGATGAACGGGAAGGCATCGATTGCGCCTGGGCGCaatcgat                                                                                                                                        <  4:104780/149‑7 (MQ=255)
            gcCAATTGGCCGCCAGTTGGCGTCAGCTGGTTCTCGCGGTTAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACAt                                                                                                                                    >  2:209091/1‑149 (MQ=255)
              cAACTGGCGGCCAATTGGCGTCAGCTGGTTCTCGCGGTTAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCt                                                                                                                                  >  8:35171/1‑149 (MQ=255)
               aaCTGGCGGCCAATTGGCGTCAGCTGGTTCTCGCGGTTAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTg                                                                                                                                 <  2:118072/149‑1 (MQ=255)
                aCTGGCGTCCAATTGGCGTCAGCTGGTTCTCGCTGTTAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGc                                                                                                                                <  2:253879/149‑1 (MQ=255)
                 cTGGCGGCCAATTGGCGTCAGCTGGTTCTCGCGGTTAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCa                                                                                                                               <  3:131414/149‑1 (MQ=255)
                    gcggcCAATTGGCGTCAGCTGGTTCTCGCGGTTAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAgg                                                                                                                             <  6:52761/148‑1 (MQ=255)
                     cggcCAATTGGCGTCAGCTGGTTCTCGCGGTTAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGAt                                                                                                                           <  2:260419/149‑1 (MQ=255)
                      ggcCAATTGGCGTCAGCTGGTTCTCGCGGTTAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGAt                                                                                                                           >  1:357668/1‑148 (MQ=255)
                      ggcCAATTGGCGTCAGCTGGTTCTCGCGGTTAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGAt                                                                                                                           >  8:142203/1‑148 (MQ=255)
                      ggcCAATTGGCGTCAGCTGGTTCTCGCGGTTAACCGCCGAAAGCTC‑‑ATGCAAAACGTTGAAGCCGTCGCTGATGGCTTTGCCATACGGCGGCTCGATGAACGCCAACGCATCGATTGAGCCCAGGCGCAGGTGAAACAACTGCACGaa                                                                                                                           >  5:235697/1‑147 (MQ=255)
                       gcCAATTGGCGTCAGCTGGTTCTCGCGGTTAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGTCTTTGCCATCCGGCGGCTCGATGAACGGGAATGAATCGTTTGCGCCCAGGCGCAGGTGCAACATCTTCAGGATCa                                                                                                                         <  5:151999/149‑1 (MQ=255)
                        ccAATTGGCGTCAGCTGGTTCTCGCGGTTAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCCCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCAc                                                                                                                        <  5:39848/149‑1 (MQ=255)
                          aaTTGGCGTCAGCTGGTTCTCGCGGTTAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCAccg                                                                                                                      <  8:159966/149‑1 (MQ=255)
                            ttGGCGTCAGCTGGTTCTCGCGGTTAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATAGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCAccgcc                                                                                                                    >  3:257318/1‑149 (MQ=255)
                                cGTCAGCTGGTTCTCGCGGTTAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCTGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCa                                                                                                                <  2:117853/149‑1 (MQ=255)
                                 gTCAGCTGGTTCTCGCGGTTAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAg                                                                                                               <  7:35171/149‑1 (MQ=255)
                                       tGGTTCTCGCGGTTAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGt                                                                                                         <  2:357668/149‑1 (MQ=255)
                                         gTTCTCGCGGTTAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAGCGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGt                                                                                                         <  3:174325/147‑1 (MQ=255)
                                         gTTCTCGCGGTTAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGt                                                                                                         <  1:209091/147‑1 (MQ=255)
                                         gTTCTCGCGGTTAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTg                                                                                                        <  7:151490/148‑1 (MQ=255)
                                             tcGCGGTTAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTgcg                                                                                                      >  3:51051/1‑146 (MQ=255)
                                               gcgGTAAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCa                                                                                                    <  6:305779/146‑1 (MQ=255)
                                                    tAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTcc                                                                                            >  8:115700/1‑149 (MQ=255)
                                                     aaCCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCg                                                                                           >  1:18777/1‑149 (MQ=255)
                                                      aCCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCgg                                                                                          <  4:59197/149‑1 (MQ=255)
                                                       ccgccgAAAGCAC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCggg                                                                                         >  1:194126/1‑149 (MQ=255)
                                                          ccgAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCg                                                                                      <  3:204716/149‑1 (MQ=255)
                                                           cgAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGc                                                                                                                                <  6:5051/106‑1 (MQ=255)
                                                           cgAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGc                                                                                                                                >  5:5051/1‑106 (MQ=255)
                                                           cgAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCgg                                                                                     <  7:142203/149‑1 (MQ=255)
                                                                  tC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGc                                                                              >  8:190767/1‑149 (MQ=255)
                                                                  tC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGACATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGACCAGGCGCAGGTGCATCATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCTGAAACTCCGCGTCGGAAAACGc                                                                              >  8:319401/1‑149 (MQ=255)
                                                                        gCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTg                                                                                                           <  1:186758/116‑1 (MQ=255)
                                                                        gCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTg                                                                                                           >  2:186758/1‑116 (MQ=255)
                                                                           aCAGGTTAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGc                                                                      >  1:345200/1‑149 (MQ=255)
                                                                             aGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCt                                                                     <  7:115700/149‑1 (MQ=255)
                                                                             aGGTGGAAGCCGTCGCTGATTGCTTTGCCATCCGGCGGCTCGATGATCTGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCGGGGTCGGTAAACGCTGGCCGGCt                                                                     <  6:235697/149‑1 (MQ=255)
                                                                               gTTGAAGCCGTCGCTGATGGCTTTGCCATCGGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGt                                                                   >  1:33744/1‑149 (MQ=255)
                                                                               gTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAAtctc                                                                                             <  7:181108/123‑1 (MQ=255)
                                                                               gTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAAtctc                                                                                             >  8:181108/1‑123 (MQ=255)
                                                                               gTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGt                                                                   >  3:335140/1‑149 (MQ=255)
                                                                               gTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGt                                                                   >  3:331171/1‑149 (MQ=255)
                                                                                ttGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGtt                                                                  <  4:51051/149‑1 (MQ=255)
                                                                                    aGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGTTGaa                                                               <  4:304404/148‑1 (MQ=255)
                                                                                     gCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCGAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGTTGAAg                                                              <  6:306608/148‑1 (MQ=255)
                                                                                      ccGTCGCTGATGGCTTTGCCATCTGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGttagggtc                                                            >  3:210159/1‑143 (MQ=255)
                                                                                      ccGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGTTGAAGTc                                                            <  5:285382/149‑1 (MQ=255)
                                                                                      ccGTCGCTGATGGCTTTGCCATCCGGCGGCACGATGAACGGGAAGGCATCGAATGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCACAATCTCCGGGACGGTAAACGCTGGCCGGCTGTTGAAGTc                                                            >  5:58121/1‑149 (MQ=255)
                                                                                             tGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGTTGAAGTCCTCTTCg                                                     >  5:141313/1‑149 (MQ=255)
                                                                                               aTGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGTTGAAGTCCTCTTCGCt                                                   <  7:298650/149‑1 (MQ=255)
                                                                                                  gCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGTTGAAGTCCTCTTCGCTGTa                                                <  2:33744/149‑1 (MQ=255)
                                                                                                       gCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGTTGAAGTCCTCTTCGCTGTAGAGac                                           >  7:312177/1‑149 (MQ=255)
                                                                                                         cATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGTTGAAGTCCTCTTCGCTGTAGAGacgc                                         <  5:173351/149‑1 (MQ=255)
                                                                                                            ccGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGGAGGTGCAACATCTGCAGCATCACCCCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACCCTGGCCGGCTGTTGAAGTCCTCTTCGCTGTAGAGacgcacg                                      >  5:262766/1‑149 (MQ=255)
                                                                                                                    tCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGTTGAAGTCCTCTTCGCTGTAGAGACGCACGCAAATGcc                              >  8:207351/1‑149 (MQ=255)
                                                                                                                    tCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGTTGAAGTCCTCTTCGCTGTAGAGACGCACGCAAATGcc                              >  8:321399/1‑149 (MQ=255)
                                                                                                                        tGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGTTGAAGTCCTCTTCGCTGTAGAGACGCACGCAAATGCCCgg                           >  6:36111/1‑148 (MQ=255)
                                                                                                                            cGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGTTGAAGTCCTCTTCGCTGTAGAGACGCACGCAAATGCCCGGCTCGa                      >  8:117355/1‑149 (MQ=255)
                                                                                                                             gggAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGTTGAAGTCCTCTTCGCTGTAGAGACGCACGCAAATGCCCGGCTCGAc                     >  3:158082/1‑149 (MQ=255)
                                                                                                                                aaGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGTTGAAGTCCTCTTCGCTGTAGAGACGCACGCAAATGCCCGGCTCGACCCg                  <  7:190767/149‑1 (MQ=255)
                                                                                                                                  ggCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGTTGAAGTCCTCTTCGCTGTAGAGACGCACGCAAATGCCCGGCTCGACCCGGc                >  7:77400/1‑149 (MQ=255)
                                                                                                                                  ggCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGTTGAAGTCCTCTTCGCTGTAGAGACGCACGCAAATGCCCGGCTCGACCCGGc                >  3:364858/1‑149 (MQ=255)
                                                                                                                                    cATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGTTGAAGTCCTCTTCGCTGTAGAGACGCACGCAAATGCCCGGCTCGACCCGGCCg              >  5:122687/1‑149 (MQ=255)
                                                                                                                                     aTCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGTTGAAGTCCTCTTCGCTGTAGAGACGCACGCAAATGCCCGGCTCGACCCGGCCg              >  1:282655/1‑148 (MQ=255)
                                                                                                                                     aTCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGTTGAAGTCCTCTTCGCTGTAGAGACGCACGCAAATGCCCGGCTCGACCCGGCCg              >  8:68962/1‑148 (MQ=255)
                                                                                                                                     aTCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGTTGAAGTCCTCTTCGCTGTAGAGACGCACGCAAATGCCCGGCTCGACCCGGCCg              >  2:393434/1‑148 (MQ=255)
                                                                                                                                     aTCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGTTGAAGTCCTCTTCGCTGTAGAGACGCACGCAAATGCCCGGCTCGACCCGGCCGc             <  6:141313/149‑1 (MQ=255)
                                                                                                                                        acctGCGCCTGGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGTTGAAGTCCTCTTCGCTGTa                                                <  6:360667/108‑1 (MQ=255)
                                                                                                                                        acctGCGCCTGGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGTTGAAGTCCTCTTCGCTGTa                                                >  5:360667/4‑111 (MQ=255)
                                                                                                                                         aTTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGCTTGGTTCGCAGAATCTCCGTGTCGGTAAAGGCTGGCCGGCTGTTGAAGTCCTCTTCGCTGTAGAGACGCACGCAAATGCCCGGCTCGACCCGGCCGCAACg         <  4:257318/149‑1 (MQ=255)
                                                                                                                                          ttGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGTTGAAGTCCTCTTCGCTGTAGAGACGCACGCAAATGCCCGGCTCGACCCGGCCGCAACgg        >  2:19227/1‑149 (MQ=255)
                                                                                                                                            gcgcCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGTTGAAGTCCTCTTCGCTGTAGAGACGCACGCAAATGCCCGGCTCGACCCGGCCGCAACGGc       <  6:19656/148‑1 (MQ=255)
                                                                                                                                               cccAGGCGCAGGTGCAACATCTGCAGGATCCCCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGGCGGCTGTAGAAGTCCTCTTCGCTGTAGAGCCGCACGCTAATGCGCGGCTATACCCGGCCGCAACGGccat    >  5:297971/1‑146 (MQ=255)
                                                                                                                                               cccAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGTTGAAGTCCTCTTCGCTGTAGAGACGCACGCAAATGCCCGGCTCGACCCGGCCGCAACGGCCttt   <  7:117355/149‑1 (MQ=255)
                                                                                                                                                ccAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGTTGAAGTCCTCTTCGCTGTAGAGACGCACGCAAATGCCCGGCTCGACCCGGCCGCAACGGCCTTTa  >  7:227238/1‑149 (MQ=255)
                                                                                                                                                ccAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGGCGGTAAACGCTGGCCGGCTGTTGAAGTCCTCTTCGCTGTAGAGACGCACGCAAATTCCCGGCTCGACCCGGCCGCAACGGCCTTTa  <  7:319401/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ATCGGCAAGCGCGCCAACTGGCGGCCAATTGGCGTCAGCTGGTTCTCGCGGTTAACCGCCGAAAGCTC‑‑CTGCAACAGGTTGAAGCCGTCGCTGATGGCTTTGCCATCCGGCGGCTCGATGAACGGGAAGGCATCGATTGCGCCCAGGCGCAGGTGCAACATCTGCAGGATCACCGCCGCCAGGTTGGTGCGCAGAATCTCCGGGTCGGTAAACGCTGGCCGGCTGTTGAAGTCCTCTTCGCTGTAGAGACGCACGCAAATGCCCGGCTCGACCCGGCCGCAACGGCCTTTA  >  NC_002947/5166387‑5166677

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: