Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,761,710 A→C 8.1% H353P (CAC→CCC)  PP_3325 → ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,761,7100AC 8.1% 215.9 / 8.5 73H353P (CAC→CCC) PP_3325ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base A (33/34);  new base C (4/2);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

TGGAAGCTGCGCCTGGCCAACCACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCC  >  NC_002947/3761562‑3761854
                                                                                                                                                    |                                                                                                                                                 
tGGAAGCTGCGCCTGGCCAACCACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCa                                                                                                                                                   >  4:77688/1‑149 (MQ=255)
    aGCTGCGCCTGGGCAACCACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACggg                                                                                                                                               <  6:102191/149‑1 (MQ=255)
     gCTGCGCCTGGCCAACCACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCACGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCCAGCTGCACggg                                                                                                                                               >  5:89709/1‑148 (MQ=255)
      cTGCGCCTGGCCAACCACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCa                                                                                                                                             <  6:91216/149‑1 (MQ=255)
             tGGCCAACCACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAACTGCCCGTGCAGCTCggc                                                                                                                                      >  4:103017/1‑147 (MQ=255)
                 cAACCACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTgg                                                                                                                                  >  8:253235/1‑149 (MQ=255)
                 cAACCACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAACTGCCCGTGCAGCTCGGCCTgg                                                                                                                                  <  3:103017/149‑1 (MQ=255)
                    ccACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGtg                                                                                                                               >  6:256973/1‑149 (MQ=255)
                    ccACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGAGAGCACCACCATCACGCGCTTCTATCGCGAACCCGACTTCGAATGGAACGACAACATCACCCACGCCGAGCTGCACGGGCACTGCGAGCTGGGtg                                                                                                                               >  4:83827/1‑149 (MQ=255)
                           acCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGCTCGAACTGCCCGTGCAGCTCggcctg                                                                                                                                   <  6:205704/138‑6 (MQ=38)
                           acCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGCTCGAACTGCCCGTGCAGCTCggcctg                                                                                                                                   >  5:205704/1‑133 (MQ=38)
                               aGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCAt                                                                                                                    <  4:373105/149‑1 (MQ=255)
                                gggCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTTCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCAt                                                                                                                    >  4:226885/1‑148 (MQ=255)
                                 ggCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCa                                                                                                                  <  3:190425/149‑1 (MQ=255)
                                   cACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGTATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGa                                                                                                                >  1:56660/1‑149 (MQ=255)
                                    aCGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGAc                                                                                                               >  6:19188/1‑149 (MQ=255)
                                      gCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGAccc                                                                                                             >  5:8781/1‑149 (MQ=255)
                                        tCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGAGACCGCGACTTCGAATGGAACGACAACATCACCAAGGCCGAGCGGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCtg                                                                                                           >  7:160064/1‑149 (MQ=255)
                                            gggCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGGTGGCAGCATCAGACCCtggtgg                                                                                                       >  4:332593/1‑149 (MQ=255)
                                              gCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTgggg                                                                                                     <  4:102719/149‑1 (MQ=255)
                                                 acagtgccGAGCCCCAGGCCCTGGTGAGCACCACCATCACGGGCTTCTATCGCGAACGCGACTTCGAATGGACCGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTg                                                                                                  <  8:160064/143‑1 (MQ=255)
                                                 aCAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCtg                                                                                                           >  8:198670/1‑140 (MQ=255)
                                                 aCAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCtg                                                                                                           <  7:198670/140‑1 (MQ=255)
                                                   aGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGa                                                                                                <  6:380474/149‑1 (MQ=255)
                                                   aGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGa                                                                                                <  5:229406/149‑1 (MQ=255)
                                                     cTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAg                                                                                               >  6:118454/1‑148 (MQ=255)
                                                     cTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGt                                                                                              <  5:355056/149‑1 (MQ=255)
                                                      tCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTa                                                                                             >  6:138167/1‑149 (MQ=255)
                                                        cGAGCACCAGGCCATGGTGAGCACCACCATCACCCGCTTCTATCGCGAACGCGACTTCGAATGGAACTACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTAc                                                                                            <  8:174499/148‑1 (MQ=255)
                                                         gAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGc                                                                                                                                              >  8:98637/1‑97 (MQ=255)
                                                              ccAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGAAGCATCAGACCCTGGTGGGGCTGGAGTACGaaaaa                                                                                      >  8:241568/1‑147 (MQ=255)
                                                              ccAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCAGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAAc                                                                                      >  5:274731/1‑148 (MQ=255)
                                                                aGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAc                                                                                   <  5:19188/149‑1 (MQ=255)
                                                                aGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAc                                                                                   <  4:195563/149‑1 (MQ=255)
                                                                  gCCCTGGTGAGCACCACCATCACGCGCTTCTATCTCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑Cg                                                                                <  8:12556/149‑1 (MQ=255)
                                                                  gCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑Cg                                                                                >  7:77327/1‑149 (MQ=255)
                                                                    ccTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGc                                                                               >  3:221120/1‑148 (MQ=255)
                                                                           aGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCa                                                                       <  6:31837/149‑1 (MQ=255)
                                                                                caccaTCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAg                                                                   >  5:144673/1‑148 (MQ=255)
                                                                                 accaTCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTa                                                                 <  6:274731/149‑1 (MQ=255)
                                                                                   caTCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCa                                                                                                                  >  4:322312/1‑99 (MQ=255)
                                                                                   caTCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCa                                                                                                                  <  3:322312/99‑1 (MQ=255)
                                                                                       aCGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGc                                                            >  5:313868/1‑148 (MQ=255)
                                                                                           gcTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGGTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAgc                                                       <  3:332593/149‑1 (MQ=255)
                                                                                           gcTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAgc                                                       <  1:155119/149‑1 (MQ=255)
                                                                                            cTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAgcg                                                      <  5:256973/149‑1 (MQ=255)
                                                                                               cTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAgcgc                                                     <  6:8781/147‑1 (MQ=255)
                                                                                                  tCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCAc                                                  <  6:805/147‑1 (MQ=255)
                                                                                                  tCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGt                                                >  5:8928/1‑149 (MQ=255)
                                                                                                     cgAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCg                                             >  7:302088/1‑149 (MQ=255)
                                                                                                       aaCGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCa                                           <  6:144673/149‑1 (MQ=255)
                                                                                                       aaCGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCa                                           <  5:118454/149‑1 (MQ=255)
                                                                                                               ttCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATgg                                   >  2:356908/1‑149 (MQ=255)
                                                                                                                 cGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATggg                                  <  5:138167/148‑1 (MQ=255)
                                                                                                                  gAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACACGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGc                                 <  1:286455/149‑1 (MQ=255)
                                                                                                                    aTGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTgg                              <  8:77327/149‑1 (MQ=255)
                                                                                                                     tGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGa                             <  6:89709/149‑1 (MQ=255)
                                                                                                                       gAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTg                               >  2:311212/1‑145 (MQ=255)
                                                                                                                       gAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGat                            <  3:77688/148‑1 (MQ=255)
                                                                                                                        aaCGACAACATCACCCAGGCCGAGCTTCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGatat                          <  3:83827/149‑1 (MQ=255)
                                                                                                                         aCGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATc                         <  1:28611/149‑1 (MQ=255)
                                                                                                                           gACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTa                       >  1:144518/1‑149 (MQ=255)
                                                                                                                               cagccACCCAGCTCGAACTGCCCGTGCAGCTCGGCCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAAc                   >  2:102765/5‑149 (MQ=37)
                                                                                                                                 gccACCCAGCTCGAACTGCCCGTGCAGCTCGGCCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAAccc                 >  3:115254/3‑149 (MQ=38)
                                                                                                                                  tCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCg                <  1:33945/149‑1 (MQ=255)
                                                                                                                                  tCACCCAGCTCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATgg                                   <  3:287120/130‑1 (MQ=255)
                                                                                                                                  tCACCCAGCTCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATgg                                   >  4:287120/1‑130 (MQ=255)
                                                                                                                                     cccAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGta             >  5:252352/1‑149 (MQ=255)
                                                                                                                                     cccAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGta             >  6:156068/1‑149 (MQ=255)
                                                                                                                                      ccAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTAccc                                                              >  7:343958/1‑99 (MQ=255)
                                                                                                                                      ccAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTAccc                                                              <  8:343958/99‑1 (MQ=255)
                                                                                                                                      ccAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGta             >  8:73579/1‑148 (MQ=255)
                                                                                                                                       cAGGCCGAGCTGCACGGGCAGTTAGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCGGGAGGACGAAAACTAC‑CGCACCAGCAAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGTCTGGATATATACAACACGGtata           >  7:356276/1‑149 (MQ=255)
                                                                                                                                                 tGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGcc  >  8:35655/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TGGAAGCTGCGCCTGGCCAACCACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCC  >  NC_002947/3761562‑3761854

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