Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,761,718 T→C 8.2% F356L (TTC→CTC)  PP_3325 → ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,761,7180TC 8.2% 204.8 / 6.5 73F356L (TTC→CTC) PP_3325ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base T (34/33);  new base C (4/2);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

TGGCCAACCACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCC  >  NC_002947/3761575‑3761864
                                                                                                                                               |                                                                                                                                                   
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)
                                                 ccAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGAAGCATCAGACCCTGGTGGGGCTGGAGTACGaaaaa                                                                                                >  8:241568/1‑147 (MQ=255)
                                                 ccAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCAGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAAc                                                                                                >  5:274731/1‑148 (MQ=255)
                                                   aGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAc                                                                                             <  4:195563/149‑1 (MQ=255)
                                                   aGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAc                                                                                             <  5:19188/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                       >  6:156068/1‑149 (MQ=255)
                                                                                                                        cccAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGta                       >  5:252352/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)
                                                                                                                                           gCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAagccag    >  1:208392/1‑149 (MQ=255)
                                                                                                                                           gCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGAGATCTACAACCCGGTATACGGCAAGCCGCagcccg    >  8:335300/1‑147 (MQ=255)
                                                                                                                                             aGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAagccagcc  <  7:73579/149‑1 (MQ=255)
                                                                                                                                               ttCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGc             >  5:155813/1‑136 (MQ=255)
                                                                                                                                               ttCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGc             <  6:155813/136‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TGGCCAACCACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAC‑CGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCC  >  NC_002947/3761575‑3761864

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