Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,863,839 C→T 7.2% P921L (CCA→CTA)  PP_3413 → two‑component system sensor histidine kinase/response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,863,8390CT 7.2% 203.9 / 4.4 69P921L (CCA→CTA) PP_3413two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base C (32/32);  new base T (2/3);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

ACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCACCTGCCCCAGGTGCGCCCTCAGCAGCCG  >  NC_002947/3863696‑3863987
                                                                                                                                               |                                                                                                                                                    
acTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCtggcag                                                                                                                                                 <  3:121276/149‑1 (MQ=255)
acTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCAAGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCcagcag                                                                                                                                                 >  1:147852/1‑149 (MQ=255)
    gCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCg                                                                                                                                             <  2:240103/149‑1 (MQ=255)
    gCGGAGAAGGGCGAGAGCGCCCTTGCACGCTGGCAGGAAGCGCGCGGCGACCTGGTGATCACCGATTGCAGCATGCCGGGCATGGACGGGCATTAACTCACCCGGCGCATTCCCAGCCGGGAGGGCGAGCGCGGCCTGCCAGCAGGcc                                                                                                                                              <  6:242363/148‑1 (MQ=255)
       gagaATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCTGGCAGGCCgcgc                                                                                                                                          <  7:43717/149‑4 (MQ=255)
          aaTGGCGAGCTCGCCCTTGCACGATGGCATGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCt                                                                                                                                       <  6:288974/149‑1 (MQ=255)
            tGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTg                                                                                                                                     >  1:203643/1‑149 (MQ=255)
            tGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTg                                                                                                                                     >  1:225953/1‑149 (MQ=255)
             ggCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTgg                                                                                                                                    >  4:174363/1‑149 (MQ=255)
              gCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTgg                                                                                                                                    <  5:38555/148‑1 (MQ=255)
              gCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGc                                                                                                                                   >  1:264809/1‑149 (MQ=255)
                  gaTCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGa                                                                                                                               >  6:221432/1‑149 (MQ=255)
                   aTCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGAc                                                                                                                              <  3:40556/149‑1 (MQ=255)
                      gCCCTTGCACGCTGGCTGGAAGGGCGCTGCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGc                                                                                                                           <  2:147852/149‑1 (MQ=255)
                      gCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGc                                                                                                                           >  6:193088/1‑149 (MQ=255)
                       cccTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCg                                                                                                                            <  8:257882/147‑1 (MQ=255)
                       cccTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCg                                                                                                                            >  7:257882/1‑147 (MQ=255)
                       cccTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGc                                                                                                                           >  3:186118/1‑148 (MQ=255)
                              aCGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGcc                                                                                                                    >  8:89162/1‑148 (MQ=255)
                               cGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCa                                                                                                                  <  2:225953/149‑1 (MQ=255)
                                    gCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGc                                                                                                                                   <  3:283349/127‑1 (MQ=255)
                                    gCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGc                                                                                                                                   >  4:283349/1‑127 (MQ=255)
                                          aGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGTGc                                                                                                       <  5:44837/149‑1 (MQ=255)
                                          aGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGc                                                                                                       <  4:52962/149‑1 (MQ=255)
                                             gcgcTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTa                                                                                                    >  5:25432/1‑149 (MQ=255)
                                               gcTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACt                                                                                                  <  4:173674/149‑1 (MQ=255)
                                                     aCCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGc                                                                                            <  6:25432/149‑1 (MQ=255)
                                                     aCCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGc                                                                                            <  5:158341/149‑1 (MQ=255)
                                                      ccTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGcc                                                                                           >  4:265711/1‑149 (MQ=255)
                                                      ccTGGTGATCACAGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGcc                                                                                           >  8:87031/1‑149 (MQ=255)
                                                       cTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCt                                                                                          <  1:114860/149‑1 (MQ=255)
                                                        tGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTg                                                                                         >  6:5860/1‑149 (MQ=255)
                                                          gTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTg                                                                                         <  5:221432/147‑1 (MQ=255)
                                                            gATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGc                                                                                       >  3:62086/1‑147 (MQ=255)
                                                                aCCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCgg                                                                                 >  5:302647/1‑149 (MQ=255)
                                                                aCCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCg                                                                                  <  5:193088/148‑1 (MQ=255)
                                                                    aTTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATg                                                                             >  3:58800/1‑149 (MQ=255)
                                                                      tGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCTGGCAGGCCgcg                                                                                                                                           >  4:271245/1‑83 (MQ=255)
                                                                      tGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCTGGCAGGCCgcg                                                                                                                                           <  3:271245/85‑3 (MQ=255)
                                                                        gttcaTGCCCGTCCATGACCGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGAc                                                                          >  1:19347/4‑148 (MQ=255)
                                                                             tGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGc                                                                    <  5:5860/149‑1 (MQ=255)
                                                                              gCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGcc                                                                   >  5:140187/1‑149 (MQ=255)
                                                                               cctgTCATGTACGGGCATGAACTCACCCGGCGCATTGGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGcc                                                                   <  7:78348/145‑1 (MQ=255)
                                                                               ccGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCt                                                                  <  7:89162/149‑1 (MQ=255)
                                                                                   tCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGc                                                                    >  5:178570/1‑143 (MQ=255)
                                                                                   tCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGc                                                                    <  6:178570/143‑1 (MQ=255)
                                                                                   tCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGtt                                                               >  1:159718/1‑148 (MQ=255)
                                                                                      tGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAg                                                           <  6:302647/149‑1 (MQ=255)
                                                                                       ggACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGGGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGc                                                          >  2:136912/1‑149 (MQ=255)
                                                                                        gACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGc                                                          <  2:264809/148‑1 (MQ=255)
                                                                                          cGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGa                                                       >  1:292307/1‑149 (MQ=255)
                                                                                           gggCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGAt                                                      >  5:1644/1‑149 (MQ=255)
                                                                                             gCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTg                                                    <  2:42868/149‑1 (MQ=255)
                                                                                             gCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCATGCCAGCAGCCCGTATTCTTGGCGTGACCGCCAAAGCCCAGGCCGAGGAGCGTACTCGCTGACAGGCCAGCGGCATGGACCAATGCCTGTTCAACCCGATTg                                                    >  1:275101/1‑149 (MQ=255)
                                                                                                   aCTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTgc                                              <  4:62086/149‑1 (MQ=255)
                                                                                                     tCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTgcg                                             >  2:28351/1‑148 (MQ=255)
                                                                                                      cACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCa                                           >  8:224079/1‑149 (MQ=255)
                                                                                                       aCCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCAc                                          >  3:140435/1‑149 (MQ=255)
                                                                                                       aCCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCAc                                          <  2:19347/149‑1 (MQ=255)
                                                                                                       aCCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCAc                                          >  5:66376/1‑149 (MQ=255)
                                                                                                         ccGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGc                                        <  7:224079/149‑1 (MQ=255)
                                                                                                              gcATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCTAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTc                                      >  7:59219/1‑146 (MQ=255)
                                                                                                               cATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGa                                  <  4:140435/149‑1 (MQ=255)
                                                                                                                     cAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCAcc                            >  8:129266/1‑149 (MQ=255)
                                                                                                                        ccTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCt                                                                  >  3:253258/1‑108 (MQ=255)
                                                                                                                        ccTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCt                                                                  <  4:253258/108‑1 (MQ=255)
                                                                                                                            gAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCACCTGcccc                      <  3:174363/148‑1 (MQ=255)
                                                                                                                              gggCGAGCGCGGCCTGCTAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCACCTGCCCCAgg                   <  8:59219/149‑1 (MQ=255)
                                                                                                                                   aGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCACCTGCCCCAGGTgcgc              >  2:304545/1‑149 (MQ=255)
                                                                                                                                     cgcgGCCTGCCAGCAGGCCGTAGTCTTGGCGTGACCGCCAATGCCGAGGCCGAGGAGGGTACTCGCTGCCTTGCCAGCTGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTTCTCACTCTCAAGACCCTCCTGCCCCAGGTtcgtg             <  8:278452/148‑3 (MQ=255)
                                                                                                                                           cTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCACCTGCCCCAGGTGCGCCCTcagca      <  4:58800/149‑1 (MQ=255)
                                                                                                                                              ccAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCACCTGCCCCAGGTGCGCCCTCAGCAGcc   <  6:140187/149‑1 (MQ=255)
                                                                                                                                               cagcagGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCACCTGCCCCAGGTGCGCCCTCAGCAGCCg  >  6:439/1‑149 (MQ=255)
                                                                                                                                               cagcagGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCACCTGCCCCAGGTGCGCCCTCAGCAGCCg  >  1:44887/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
ACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCACCTGCCCCAGGTGCGCCCTCAGCAGCCG  >  NC_002947/3863696‑3863987

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: