Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,761,701 C→T 8.0% A350V (GCC→GTC) ‡ PP_3325 → ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,761,7010CT 8.0% 218.9 / 6.4 75A350V (GCC→GTC) ‡PP_3325ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base C (34/35);  new base T (4/2);  total (38/37)
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.93e-01

TGGAAGCTGCGCCTGGCCAACCACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATA  >  NC_002947/3761562‑3761845
                                                                                                                                           |                                                                                                                                                 
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                                                                                                                          >  5:205704/1‑133 (MQ=38)
                           acCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGCTCGAACTGCCCGTGCAGCTCggcctg                                                                                                                          <  6:205704/138‑6 (MQ=38)
                             ccAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGc                                                                                                                                              >  1:267778/1‑116 (MQ=255)
                             ccAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGc                                                                                                                                              <  2:267778/116‑1 (MQ=255)
                               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                                                                                       <  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)
                                                                  gCCCTGGTGAGCACCACCATCACGCGCTTCTATCTCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCg                                                                       <  8:12556/149‑1 (MQ=255)
                                                                  gCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCg                                                                       >  7:77327/1‑149 (MQ=255)
                                                                    ccTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGc                                                                      >  3:221120/1‑148 (MQ=255)
                                                                           aGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCa                                                              <  6:31837/149‑1 (MQ=255)
                                                                                caccaTCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAg                                                          >  5:144673/1‑148 (MQ=255)
                                                                                 accaTCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTa                                                        <  6:274731/149‑1 (MQ=255)
                                                                                   caTCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCa                                                                                                         <  3:322312/99‑1 (MQ=255)
                                                                                   caTCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCa                                                                                                         >  4:322312/1‑99 (MQ=255)
                                                                                       aCGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGc                                                   >  5:313868/1‑148 (MQ=255)
                                                                                           gcTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGGTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAgc                                              <  3:332593/149‑1 (MQ=255)
                                                                                           gcTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAgc                                              <  1:155119/149‑1 (MQ=255)
                                                                                            cTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAgcg                                             <  5:256973/149‑1 (MQ=255)
                                                                                               cTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAgcgc                                            <  6:8781/147‑1 (MQ=255)
                                                                                                  tCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCAc                                         <  6:805/147‑1 (MQ=255)
                                                                                                  tCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGt                                       >  5:8928/1‑149 (MQ=255)
                                                                                                     cgAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCg                                    >  7:302088/1‑149 (MQ=255)
                                                                                                       aaCGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCa                                  <  6:144673/149‑1 (MQ=255)
                                                                                                       aaCGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCa                                  <  5:118454/149‑1 (MQ=255)
                                                                                                               ttCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATgg                          >  2:356908/1‑149 (MQ=255)
                                                                                                                 cGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATggg                         <  5:138167/148‑1 (MQ=255)
                                                                                                                  gAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACACGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGc                        <  1:286455/149‑1 (MQ=255)
                                                                                                                    aTGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTgg                     <  8:77327/149‑1 (MQ=255)
                                                                                                                     tGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGa                    <  6:89709/149‑1 (MQ=255)
                                                                                                                       gAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTg                      >  2:311212/1‑145 (MQ=255)
                                                                                                                       gAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGat                   <  3:77688/148‑1 (MQ=255)
                                                                                                                        aaCGACAACATCACCCAGGCCGAGCTTCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGatat                 <  3:83827/149‑1 (MQ=255)
                                                                                                                         aCGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATc                <  1:28611/149‑1 (MQ=255)
                                                                                                                           gACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTa              >  1:144518/1‑149 (MQ=255)
                                                                                                                               cagccACCCAGCTCGAACTGCCCGTGCAGCTCGGCCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAAc          >  2:102765/5‑149 (MQ=37)
                                                                                                                                 gccACCCAGCTCGAACTGCCCGTGCAGCTCGGCCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAAccc        >  3:115254/3‑149 (MQ=38)
                                                                                                                                  tCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCg       <  1:33945/149‑1 (MQ=255)
                                                                                                                                  tCACCCAGCTCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATgg                          <  3:287120/130‑1 (MQ=255)
                                                                                                                                  tCACCCAGCTCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATgg                          >  4:287120/1‑130 (MQ=255)
                                                                                                                                     cccAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGta    >  5:252352/1‑149 (MQ=255)
                                                                                                                                     cccAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGta    >  6:156068/1‑149 (MQ=255)
                                                                                                                                      ccAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTAccc                                                     >  7:343958/1‑99 (MQ=255)
                                                                                                                                      ccAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTAccc                                                     <  8:343958/99‑1 (MQ=255)
                                                                                                                                      ccAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGta    >  8:73579/1‑148 (MQ=255)
                                                                                                                                       cAGGCCGAGCTGCACGGGCAGTTAGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCGGGAGGACGAAAACTACCGCACCAGCAAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGTCTGGATATATACAACACGGtata  >  7:356276/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                 
TGGAAGCTGCGCCTGGCCAACCACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATA  >  NC_002947/3761562‑3761845

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: