Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,863,777 G→C 5.2% P900P (CCG→CCC 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,7770GC 5.2% 316.3 / 5.9 96P900P (CCG→CCCPP_3413two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base G (46/45);  new base C (3/2);  total (49/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

TACCCGGCCAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTG  >  NC_002947/3863637‑3863924
                                                                                                                                            |                                                                                                                                                    
tACCCGGCCAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGa                                                                                                                                              >  4:41378/1‑149 (MQ=255)
tACCCGGCCAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGa                                                                                                                                              >  3:359475/1‑149 (MQ=255)
tACCCGGCCAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCCGTccatga                                                                                                                                              >  7:299725/1‑144 (MQ=255)
   ccGGCCAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGAcag                                                                                                                                           >  6:280923/1‑147 (MQ=255)
   ccGGCCAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACgg                                                                                                                                           >  3:49782/1‑149 (MQ=255)
     ggCCAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACAGGc                                                                                                                                         <  5:280923/149‑1 (MQ=255)
       ccAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCAt                                                                                                                                       <  8:313414/149‑1 (MQ=255)
         aaCCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGa                                                                                                                                     <  5:36513/149‑1 (MQ=255)
          aCCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCCGTCCATGACCGGCATgct                                                                                                                                    <  8:299725/149‑3 (MQ=255)
          aCCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCCGTCCATGACCGGCATgct                                                                                                                                    <  5:68359/149‑3 (MQ=255)
            cTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCAATTGCAGCATGCCGGTCATGGACGGGCATGAACt                                                                                                                                  >  7:128756/1‑149 (MQ=255)
              gCTTTTGCTGGAACGGCAACTGCAAACCCTGGAGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGc                                                                                                                                         <  6:317118/140‑1 (MQ=255)
              gCTTTTGCTGGAACGGCAACTGCAAACCCTGGAGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACAGATTGCAGCATGCCGGTCATGGACGGGc                                                                                                                                         >  5:317118/1‑140 (MQ=255)
                    gCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATgg                                                                                                                                               >  5:306989/1‑128 (MQ=255)
                    gCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATgg                                                                                                                                               <  6:306989/128‑1 (MQ=255)
                     cTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGgcg                                                                                                                        >  5:7657/1‑149 (MQ=255)
                      tGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCACGAAGCGCGCTACGACCTGGTCATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGgcg                                                                                                                        >  7:107708/1‑148 (MQ=255)
                       ggAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGAACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCa                                                                                                                      >  8:334076/1‑149 (MQ=255)
                        gAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCAt                                                                                                                     >  4:118859/1‑149 (MQ=255)
                         aaCCGCAACTGCAAACCCTGCGGCACCATGTGACACTGCCGGAGAATGGCGAGATCGCCCTACCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGAAGCATGACGGTCATGGACGGGCATGAACTCACCCGGCGCAtt                                                                                                                    >  4:348517/1‑149 (MQ=255)
                             gCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATg                                                                                                                                      <  6:366881/128‑1 (MQ=255)
                             gCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATg                                                                                                                                      >  5:366881/1‑128 (MQ=255)
                              cAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCAt                                                                                                                     <  4:184113/143‑1 (MQ=255)
                              cAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCAt                                                                                                                     >  3:184113/1‑143 (MQ=255)
                                aCTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGcc                                                                                                             <  2:306615/149‑1 (MQ=255)
                                 cTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCAAGAACTCACCCGGCGCATTCGCAGCCt                                                                                                            <  4:359475/149‑1 (MQ=255)
                                    cAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCTCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTg                                                                                                           >  6:307597/1‑147 (MQ=255)
                                    cAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTg                                                                                                           >  3:370654/1‑147 (MQ=255)
                                    cAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGa                                                                                                         <  3:67797/149‑1 (MQ=255)
                                         ccTGGGGCACCATGTGACACTGGCGGAGACTGGCGAGATCGCCCTTGCACGCTGGCCGGAAGCGCGCTTCGACCTGGTGATAACCGATTGCAGCATGCCGTTCATGGCCGGGCATGAACTCACCCCGCGCATTAGCAGCCTGGAGGGc                                                                                                     >  6:100703/1‑148 (MQ=255)
                                         ccTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCGGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGc                                                                                                     >  6:10441/1‑148 (MQ=255)
                                             gggCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAgcg                                                                                                >  6:36021/1‑149 (MQ=255)
                                               gccccATGGGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGcgcg                                                                                              <  6:372090/146‑1 (MQ=255)
                                               gCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGAATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGcgcg                                                                                              >  5:117291/1‑149 (MQ=255)
                                                 aCCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGc                                                                                            <  6:7657/149‑1 (MQ=255)
                                                     tgtgACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGc                                                                                        >  4:90367/1‑149 (MQ=255)
                                                       tgACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCca                                                                                      <  2:164433/149‑1 (MQ=255)
                                                         acacTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCcagc                                                                                    <  1:363737/149‑1 (MQ=255)
                                                         acacTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCcag                                                                                     >  6:33077/1‑148 (MQ=255)
                                                         acacTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCcag                                                                                     <  7:48332/148‑1 (MQ=255)
                                                           acTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCcagcag                                                                                  >  7:359752/1‑149 (MQ=255)
                                                               gCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGcgcg                                                                                              <  4:60594/133‑1 (MQ=255)
                                                               gCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGcgcg                                                                                              >  3:60594/1‑133 (MQ=255)
                                                               gCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGcc                                                                               <  7:338534/148‑1 (MQ=255)
                                                               gCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCg                                                                              >  8:91000/1‑149 (MQ=255)
                                                               gCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCg                                                                              >  6:241403/1‑149 (MQ=255)
                                                                 ggAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTa                                                                            <  3:90367/149‑1 (MQ=255)
                                                                   agaATGGCGAGATCGCTCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTAtt                                                                          <  5:307597/149‑1 (MQ=255)
                                                                   agaATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTAtt                                                                          <  5:241403/149‑1 (MQ=255)
                                                                   agaATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTAt                                                                           >  8:269488/1‑148 (MQ=255)
                                                                    gaATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTc                                                                         >  8:6231/1‑149 (MQ=255)
                                                                         gCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAgg                                                                                 <  5:312120/136‑1 (MQ=255)
                                                                         gCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAgg                                                                                 >  6:312120/1‑136 (MQ=255)
                                                                           gagaTCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGt                                                                  <  4:144474/149‑1 (MQ=255)
                                                                            agaTCGCCCTTGCACGCTGGCAGGAAGCGCGCTTAGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGt                                                                  >  4:308098/1‑148 (MQ=255)
                                                                               tCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGAcc                                                              >  6:187768/1‑149 (MQ=255)
                                                                                cGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGcc                                                                               >  7:185485/1‑131 (MQ=255)
                                                                                cGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGcc                                                                               <  8:185485/131‑1 (MQ=255)
                                                                                 gCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGc                                                            >  2:83135/1‑149 (MQ=255)
                                                                                  cccTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGcc                                                           <  7:91000/149‑1 (MQ=255)
                                                                                    cTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCaa                                                         <  5:36021/149‑1 (MQ=255)
                                                                                    cTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCaa                                                         <  3:41378/149‑1 (MQ=255)
                                                                                     atGCACTCTGGCAGGAATCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGTGCATGAACTCACCCGGCGCAGGCGCAGCCTGGAGGGCGAGCTCGGCCTGCCAGCAGGCCTTCTTCTTGGCGTGACCGCCAAt                                                        <  5:100703/148‑1 (MQ=255)
                                                                                      tGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGc                                                            <  4:61285/144‑1 (MQ=255)
                                                                                      tGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCagat                                                       >  3:61285/1‑146 (MQ=255)
                                                                                       gCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATg                                                       <  5:33077/148‑1 (MQ=255)
                                                                                         aCGCTGGCAGGAAGCGCACTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGcc                                                     >  1:330731/1‑148 (MQ=255)
                                                                                          cGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTCACCGCCAATGCCCa                                                   <  7:334076/149‑1 (MQ=255)
                                                                                            cTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAgg                                                 >  1:341457/1‑149 (MQ=255)
                                                                                               gCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGc                                                <  7:6231/147‑1 (MQ=255)
                                                                                               gCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGc                                                <  1:83135/147‑1 (MQ=255)
                                                                                                cAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCga                                             <  1:15310/149‑1 (MQ=255)
                                                                                                 aGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCgag                                            >  7:40661/1‑149 (MQ=255)
                                                                                                   gAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGTCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCATGCCgagga                                          <  3:348517/149‑1 (MQ=255)
                                                                                                     aGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGc                                        <  7:269488/149‑1 (MQ=255)
                                                                                                      gcgcgcTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCg                                       <  2:341457/149‑1 (MQ=255)
                                                                                                        gcacTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTa                                     <  2:330731/146‑1 (MQ=255)
                                                                                                                aCCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGc                             <  4:49782/149‑1 (MQ=255)
                                                                                                                  cTGGTGATCACCGATTGCAGCATGCCGGTCATGTTCTGGCATGAACTCACCCGGCGCATTCGCAGCCTGGTGTGCGAGCGCTGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCt                           <  8:359752/149‑1 (MQ=255)
                                                                                                                   tGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTg                          >  6:191326/1‑149 (MQ=255)
                                                                                                                     gTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGc                        >  3:134683/1‑149 (MQ=255)
                                                                                                                       gATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGc                        >  5:58699/1‑147 (MQ=255)
                                                                                                                         tCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGc                    >  4:113470/1‑149 (MQ=255)
                                                                                                                          cACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATg                                                       >  8:139627/1‑113 (MQ=255)
                                                                                                                          cACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATg                                                       <  7:139627/113‑1 (MQ=255)
                                                                                                                           aCCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTg                              >  3:359703/1‑137 (MQ=255)
                                                                                                                           aCCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTg                              <  4:359703/137‑1 (MQ=255)
                                                                                                                           aCCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCgg                  >  2:26120/1‑149 (MQ=255)
                                                                                                                           aCCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCgg                  >  4:287702/1‑149 (MQ=255)
                                                                                                                                   gttcaTGCCCGTCCATGACCGGCATGAACTCACCCGGACGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGAc           >  1:267386/4‑149 (MQ=255)
                                                                                                                                   gttcaTGCCCGTCCATGACCGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGAc           >  7:440/4‑148 (MQ=255)
                                                                                                                                   cagcaTGCCGGTCATGGACGGGCATGACCTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGAc           <  5:10441/148‑1 (MQ=255)
                                                                                                                                    agcaTGCCGGTCATGGACGGGCATGAACTCACCCGGTGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACg          >  6:24816/1‑148 (MQ=255)
                                                                                                                                    agcaTGCAGGTCATGGAGGGGCTTGAACTCACCCGGCGCATTGGCAGGCTGGAGGGCGAGCGCGGCGTGGCAGCAGGCCGTATTCTTGGCGTGACCGCGAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGa         <  8:107708/149‑1 (MQ=255)
                                                                                                                                        tGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGc     <  3:113470/149‑1 (MQ=255)
                                                                                                                                           cGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTg  <  3:118859/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
TACCCGGCCAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTG  >  NC_002947/3863637‑3863924

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: