Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,517,377 T→C 5.5% N396N (AAT→AAC)  PP_4849 → molecular chaperone DnaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,517,3770TC 5.5% 215.0 / 2.1 72N396N (AAT→AAC) PP_4849molecular chaperone DnaK
Reads supporting (aligned to +/- strand):  ref base T (30/38);  new base C (2/2);  total (32/40)
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.76e-01

AAGGCCTGCTGGAGCGCGTGCGCGGTAGCGTGACCGAACTGCTGGCCAAAGCTGGGGTAAGCGAAGGCCAGGTCGATACGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCGTTTAGCCTGTACCCGCGAAAGGGCCGGTACAGGCCACTACATGATCAGACCAGCTCC  >  NC_002947/5517229‑5517525
                                                                                                                                                    |                                                                                                                                                    
aaGGCCTGCTGGAGCGCGTGCGCGGTAGCGTGACCGAACTGCTGGCCAAAGCTGGGGTAAGCGAAGGCCAGGTCGATACGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAAt                                                                                                                                                      >  1:301834/1‑149 (MQ=255)
    ccTGCTGGAGCGCGTGCGCGGTAGCGTGACCGAACTGCTGGCCAAAGCTGGGGTAAGCGAAGGCCAGGTCGATACGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCAc                                                                                                                                                  <  3:155694/149‑1 (MQ=255)
     ctgctgGAGCGCGTGCGCGGTAGCGTGACCGAACTGCTGGCCAAAGCTGGGGTAAGCGAAGGCCAGGTCGATACGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACg                                                                                                                                                 >  8:135496/1‑149 (MQ=255)
        ctgGAGCGCGTGCGCGGTAGCGTGACCGAACTGCTGGCCAAAGCT‑GGGTAAGCGAAGGCCAGGTCGATACGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCAt                                                                                                                                             >  8:255396/1‑149 (MQ=255)
           gAGCGCGTGCGCGGTAGCGTGACCGAACTGCTGGCCAAAGCTGGGGTAAGCGAAGGCCAGGTCGATACGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCAtgt                                                                                                                                           >  5:200650/1‑149 (MQ=255)
           gAGCGCGTGCGCGGTAGCGTGACCGAACTGCTGGCCAAAGCTGGGGTAAGCGAAGGCCAGGTCGATACGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCAt                                                                                                                                             <  3:227894/147‑1 (MQ=255)
                 gTGCGCGGTAGCGTGACCGAACTGCTGGCCAAAGCTGGGGTAAGCGAAGGCCAGGTCGATACGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAgg                                                                                                                                     >  2:50479/1‑149 (MQ=255)
                  tGCGCGGTAGCGTGACCGAACTGCTGGCCAAAGCTGGGGTAAGCGAAGGCCAGGTCGATACGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGc                                                                                                                                    <  7:8536/149‑1 (MQ=255)
                     gcgGTAGCGTGACCGAACTGCTGGCCAAAGCTGGGGTAAGCGAAGGCCAGGTCGATACGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAat                                                                                                                                 <  4:89547/149‑1 (MQ=255)
                       ggTAGCGTGACCGAACTGCTGGCCAAAGCTGGGGTAAGCGAAGGCCAGGTCGATACGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAACGCACGGCATGTGGAAGGCAatat                                                                                                                               >  1:63263/1‑149 (MQ=255)
                       ggTAGCGTGACCGAACTGCTGGCCAAAGCTGGGGTAAGCGAAGGCAAGGTCGATACGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAACGCACGGCATGTGGAAGGCAatat                                                                                                                               >  1:67527/1‑149 (MQ=255)
                             gTGACCGAACTTCTGGCCAAATCTTGTGTAAGCTAAGGCCAGGTCGATACGGTGTTCTTCACCGGTGGCTCCAGCGTCATTCCGGCGCTGCGCAACAGTGTTTCGGCTATGCTGCCGAATGCACGGTATGTGGAAGGCAATATCTTTgg                                                                                                                         <  7:193753/149‑1 (MQ=255)
                                  cGAACTGCTGGCCAAAGCTGGGGTAAGCGAAGGCCAGGTCGATACGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGgcagca                                                                                                                    >  3:243745/1‑149 (MQ=255)
                                  cGAACTGCTGGCCAAAGCTGGGGTAAGCGAAGGCCAGGTCGATACGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGgcagca                                                                                                                    <  2:151760/149‑1 (MQ=255)
                                     aCTGCTGGCCAAAGCT‑GGGTAAGCGAAGGCCAGGTCGATACGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCGGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTgg                                                                                                                <  7:255396/149‑1 (MQ=255)
                                       tgctgGCCAAAGCTGGGGTAAGCGAAGGCCAGGTCGATACGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGt                                                                                                               >  4:37886/1‑149 (MQ=255)
                                           ggCCAAAGCTGGGGTAAGCGAAGGCCAGGTCGATACGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCg                                                                                                           >  1:84680/1‑149 (MQ=255)
                                            gCCAAAGCTGGGGTAAGCGAAGGCCAGGTCGATACGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCgg                                                                                                          <  4:49555/149‑1 (MQ=255)
                                               aaaGCTGGGGTAAGCGAAGGCCAGGTCGATACGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCt                                                                                                       <  7:135496/149‑1 (MQ=255)
                                                aaGCTGGGGTAAGCGAAGGCCAGGTCGATACGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTg                                                                                                      <  5:145421/149‑1 (MQ=255)
                                                aaGCTGGGGTAAGCGAAGGCCAGGTCGATACGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTg                                                                                                      <  5:311208/149‑1 (MQ=255)
                                                      gggTAAGCGAAGGCCAGGTCGATACGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCGGCATTGGTAGCGGGCTGGCCATc                                                                                                >  2:61388/1‑149 (MQ=255)
                                                           aGCGAAGGCCAGGTCGATACGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGgc                                                                                           >  8:183767/1‑149 (MQ=255)
                                                                    cAGGTCGATACGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAACGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCg                                                                                  <  2:63263/149‑1 (MQ=255)
                                                                    cAGGTCGATACGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAACGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCg                                                                                  <  2:67527/149‑1 (MQ=255)
                                                                       gTCGATACGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTa                                                                               <  4:243745/149‑1 (MQ=255)
                                                                             aCGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGgcgc                                                                         >  3:20573/1‑149 (MQ=255)
                                                                             aCGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGgcgc                                                                         >  3:20572/1‑149 (MQ=255)
                                                                              cGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGcgcg                                                                        <  3:37886/149‑1 (MQ=255)
                                                                               ggTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCgg                                                                       >  2:182595/1‑149 (MQ=255)
                                                                                 tgtTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGcc                                                                     >  3:5516/1‑149 (MQ=255)
                                                                                   ttcttcACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGgcgc                                                                         >  2:250700/1‑143 (MQ=255)
                                                                                   ttcttcACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCt                                                                    >  2:178323/1‑148 (MQ=255)
                                                                                   ttcttcACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTATCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGgcgc                                                                         <  1:250700/143‑1 (MQ=255)
                                                                                     cttcACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGa                                                                  <  8:258216/148‑1 (MQ=255)
                                                                                     cttcACCGGTGGCTACAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGa                                                                  <  7:183767/148‑1 (MQ=255)
                                                                                        cACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGAttt                                                               <  7:32786/148‑1 (MQ=255)
                                                                                        cACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGAt                                                                 <  2:84680/146‑1 (MQ=255)
                                                                                          ccGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTc                                                            <  2:301834/149‑1 (MQ=255)
                                                                                           cGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCg                                                           >  5:231478/1‑149 (MQ=255)
                                                                                           cGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCg                                                           >  1:101635/1‑149 (MQ=255)
                                                                                                 cTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCGTTTAGc                                                     <  4:8149/149‑1 (MQ=255)
                                                                                                 cTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCGTTTAGc                                                     <  7:254173/149‑1 (MQ=255)
                                                                                                  tCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCGTTTAGcc                                                    >  1:151233/1‑149 (MQ=255)
                                                                                                  tCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCATTTAGcc                                                    >  5:194046/1‑149 (MQ=255)
                                                                                                  tCCAGCGGCATTCCGGCGCTGCGCAACAGAGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCGTTTAGcc                                                    >  3:279160/1‑149 (MQ=255)
                                                                                                   ccAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCt                                                                    >  4:176198/1‑132 (MQ=255)
                                                                                                   ccAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCt                                                                    <  3:176198/132‑1 (MQ=255)
                                                                                                     aGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGAATTGTCGTTTAGCCTg                                                  >  2:98575/1‑148 (MQ=255)
                                                                                                       cGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCGTTTAGCCTGTAc                                               >  7:106766/1‑149 (MQ=255)
                                                                                                           aTTCCGGCGCTGCGCAACAGTGTTTCGGTGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCGTTTAGCCTGTACCcgc                                           <  6:174820/149‑1 (MQ=255)
                                                                                                                    cTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCGTTTAGCCTGTACCCGCGAAAGGGcc                                  <  2:151233/149‑1 (MQ=255)
                                                                                                                     tGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCGTTTAGCCTGTAccc                                             >  5:176351/1‑137 (MQ=255)
                                                                                                                     tGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCGTTTAGCCTGTAccc                                             <  6:176351/137‑1 (MQ=255)
                                                                                                                       cgcAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCGTTTAGCCTGTACCCGCGAAAGGGCCgg                                <  1:178323/148‑1 (MQ=255)
                                                                                                                         cAACAGTGTTTCGGCGATGCTGCCGAATGCACTGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCGTTTAGCCTGTACCcgcg                                          <  2:37945/136‑1 (MQ=255)
                                                                                                                         cAACAGTGTTTCGGCGATGCTGCCGAATGCACCGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCGTTTAGCCTGTACCcgcg                                          >  1:37945/1‑136 (MQ=255)
                                                                                                                          aaCAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCGTTTAGCCTGTACCCGCGAAAGGGCCGGTa                              >  8:197815/1‑147 (MQ=255)
                                                                                                                          aaCAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCGTTTAGCCTGTACCCGCGAAAGGGCCGGTa                              >  8:199540/1‑147 (MQ=255)
                                                                                                                          aaCAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCG‑TACGGCGCGGCCTGATTTGTCGTTTAGCCTGTACCCGCGAAAGGGCCGGTACAg                           >  4:275778/1‑149 (MQ=255)
                                                                                                                          aaCAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCG‑TACGGCGCGGCCTGATTTGTCGTTTAGCCTGTACCCGCGAAAGGGCCGGTACAg                           >  4:271311/1‑149 (MQ=255)
                                                                                                                          aaCAGAGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCGTTTAGCCTGTACCCGCGAAAGGGCCGGTACa                            <  4:279160/149‑1 (MQ=255)
                                                                                                                           tcAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCGTTTAGCCTGTACCCGCGAAAGGGCCGGTACAg                           <  6:138261/148‑1 (MQ=255)
                                                                                                                           aCAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCGTTTAGCCTGTACCCGCGAAAGGGCCGGTACAg                           <  2:101635/149‑1 (MQ=255)
                                                                                                                            cAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCGGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCGTTTAGCCTGTACCCGCGAAAGGGCCGGTACAgg                          <  1:61388/149‑1 (MQ=255)
                                                                                                                            cAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCGTTTAGCCTGTACCCGCGAAAGGGCCGGTACAgg                          <  7:197815/149‑1 (MQ=255)
                                                                                                                             aGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCGTTTAGCCTGTACCCGCGAAAGGGCCGGTACAgga                         <  7:199540/149‑2 (MQ=255)
                                                                                                                             aGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCGTTTAGCCTGTACCCGCGAAAGGGCCGGTACAGGc                         <  8:45986/149‑1 (MQ=255)
                                                                                                                               tgtTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCGTTTAGCCTGTACCCGCGAAAGGGCCGGTACAGGCCa                       >  6:170968/1‑149 (MQ=255)
                                                                                                                                     ggCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCGTTTAGCCTGTACCCGCGAAAGGGCCGGTACAGGCCACTACAt                 <  8:106766/149‑1 (MQ=255)
                                                                                                                                        gATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCGTTTAGCCTGTACCCGCGAAAGGGCCGGTACAGGCCACTACATGAt              <  6:231478/149‑1 (MQ=255)
                                                                                                                                               ccGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCTTTACGGCGCGGCCTGATTTGTCGTTTAGCCTGTACCCGCGAAAGTGCCGGTACAGGCCACTACATGATCAGAc         <  5:170968/147‑1 (MQ=255)
                                                                                                                                               ccGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCGTTTAGCCTGTACCCGCGAAAGGGCCGGTACAGGCCACTACATGATCAGACCa       <  1:182595/149‑1 (MQ=255)
                                                                                                                                                    tGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCGTTTAGCCTGTACCCGCGAAAGGGCCGGTACAGGCCACTACATGATCAGACCAGCTcc  <  6:200650/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAGGCCTGCTGGAGCGCGTGCGCGGTAGCGTGACCGAACTGCTGGCCAAAGCTGGGGTAAGCGAAGGCCAGGTCGATACGGTGTTCTTCACCGGTGGCTCCAGCGGCATTCCGGCGCTGCGCAACAGTGTTTCGGCGATGCTGCCGAATGCACGGCATGTGGAAGGCAATATCTTTGGCAGCATTGGTAGCGGGCTGGCCATCGAGGCGCGCAAGCGTTACGGCGCGGCCTGATTTGTCGTTTAGCCTGTACCCGCGAAAGGGCCGGTACAGGCCACTACATGATCAGACCAGCTCC  >  NC_002947/5517229‑5517525

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