Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,761,713 G→T 8.3% G354V (GGG→GTG)  PP_3325 → ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,761,7130GT 8.3% 202.5 / 7.3 71G354V (GGG→GTG) PP_3325ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base G (32/33);  new base T (4/2);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

AGCTGCGCCTGGCCAACCACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGC  >  NC_002947/3761566‑3761859
                                                                                                                                                   |                                                                                                                                                   
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                                                                                                                                        >  5:205704/1‑133 (MQ=38)
                       acCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGCTCGAACTGCCCGTGCAGCTCggcctg                                                                                                                                        <  6:205704/138‑6 (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                                                                                                                <  7:198670/140‑1 (MQ=255)
                                             aCAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCtg                                                                                                                >  8:198670/1‑140 (MQ=255)
                                               aGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGa                                                                                                     <  5:229406/149‑1 (MQ=255)
                                               aGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGa                                                                                                     <  6:380474/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                                                                                                                       <  3:322312/99‑1 (MQ=255)
                                                                               caTCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCa                                                                                                                       >  4:322312/1‑99 (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)
                                                                                                                                                   ggCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAagc  >  4:351600/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGCTGCGCCTGGCCAACCACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGC  >  NC_002947/3761566‑3761859

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