Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,336,185 C→G 5.5% A872G (GCC→GGC)  PP_4695 → two‑component system sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,336,1850CG 5.5% 237.8 / 3.3 73A872G (GCC→GGC) PP_4695two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base C (32/37);  new base G (2/2);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

ACCAGAACAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAACATCATGGA  >  NC_002947/5336038‑5336329
                                                                                                                                                   |                                                                                                                                                 
accaGAACAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGcc                                                                                                                                                  >  5:22468/1‑149 (MQ=255)
 ccaGAACAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGccc                                                                                                                                                 <  1:138190/149‑1 (MQ=255)
 ccaGAACAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGcc                                                                                                                                                  >  8:292598/1‑148 (MQ=255)
  caGAACAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGcc                                                                                                                                                  >  5:77717/1‑147 (MQ=255)
    gAACAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAgg                                                                                                                                              <  8:177504/149‑1 (MQ=255)
    gAACAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAgg                                                                                                                                              <  1:67072/149‑1 (MQ=255)
       cAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGCTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGc                                                                                                                                           >  5:29368/1‑149 (MQ=255)
        aGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCt                                                                                                                                          <  5:361502/149‑1 (MQ=255)
          cGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGa                                                                                                                                        <  1:185197/149‑1 (MQ=255)
          cGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGa                                                                                                                                        <  7:328561/149‑1 (MQ=255)
          cGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGa                                                                                                                                        <  1:375613/149‑1 (MQ=255)
          cGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGa                                                                                                                                        <  2:303468/149‑1 (MQ=255)
          cGAGGAGCCCGGGTGCCTGGCTGAAGTGGCACAAGATGCCAGCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCACCTGGGCCAGGCGCTgc                                                                                                                                        <  2:225908/149‑2 (MQ=255)
             ggagCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTATAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATc                                                                                                                                      >  4:23199/1‑148 (MQ=255)
               agCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAAc                                                                                                                                   <  4:110021/149‑1 (MQ=255)
                gcacGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAATGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAAcc                                                                                                                                  <  3:315283/146‑1 (MQ=255)
                   cGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGc                                                                                                                               >  4:262278/1‑149 (MQ=255)
                         ccTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCACCTGGGCCAGGCGCTgc                                                                                                                                        <  7:239324/134‑2 (MQ=255)
                         ccTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCACCTGGGCCAGGCGCTgc                                                                                                                                        >  8:239324/1‑133 (MQ=255)
                             gCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTTGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGc                                                                                                                     <  3:178070/149‑1 (MQ=255)
                               tGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGccc                                                                                                                   >  5:305423/1‑149 (MQ=255)
                                   gTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTATAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGa                                                                                                               <  3:23199/149‑1 (MQ=255)
                                   gTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGa                                                                                                               <  1:144977/149‑1 (MQ=255)
                                    tGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGAt                                                                                                              <  1:336383/149‑1 (MQ=255)
                                      gCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGc                                                                                                            >  3:346308/1‑149 (MQ=255)
                                        acaAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCACCGCCCGCGATGCCt                                                                                                          >  8:72140/1‑149 (MQ=255)
                                        acaAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCt                                                                                                          >  5:358210/1‑149 (MQ=255)
                                           aGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTcgc                                                                                                       <  7:23156/149‑1 (MQ=255)
                                               gCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCgccgcc                                                                                                   >  5:90677/1‑149 (MQ=255)
                                                ccATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGccc                                                                                                  <  7:122436/149‑1 (MQ=255)
                                                                     aCCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCt                                                                                                                                 >  4:102840/1‑97 (MQ=255)
                                                                     aCCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCt                                                                                                                                 <  3:102840/97‑1 (MQ=255)
                                                                       cGGCGCAATTTCGAAGTACAGTTCATCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAg                                                                          >  7:51738/1‑149 (MQ=255)
                                                                         gcgcAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGc                                                                         >  1:74172/1‑148 (MQ=255)
                                                                         gcgcAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCg                                                                        <  7:292598/149‑1 (MQ=255)
                                                                          cgcAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGc                                                                         >  3:233661/1‑147 (MQ=255)
                                                                           gcAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCCTTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAg                                                                      <  7:168847/149‑1 (MQ=255)
                                                                           gcAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAg                                                                      >  4:110932/1‑149 (MQ=255)
                                                                           gcAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAg                                                                      <  7:72140/149‑1 (MQ=255)
                                                                              aTTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGt                                                                    >  3:352476/1‑148 (MQ=255)
                                                                                ttCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGGCCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTaa                                                                 >  5:85368/1‑149 (MQ=255)
                                                                                 tCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAAc                                                                >  7:314511/1‑149 (MQ=255)
                                                                                   gAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGa                                                              <  1:338627/149‑1 (MQ=255)
                                                                                   gAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGa                                                              <  4:346308/149‑1 (MQ=255)
                                                                                           gTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGcaca                                                         >  4:53467/1‑146 (MQ=255)
                                                                                           gTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACAcc                                                       >  2:287779/1‑148 (MQ=255)
                                                                                           gTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACAcc                                                       <  1:287779/148‑1 (MQ=255)
                                                                                               ttcAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGa                                                  >  5:286014/1‑149 (MQ=255)
                                                                                               ttcAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGa                                                  >  1:78823/1‑149 (MQ=255)
                                                                                                    ccTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGGTCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGa                                             <  1:263312/149‑1 (MQ=255)
                                                                                                        tGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGAAGCGAGGTTAACGAGCACACCGTCGACCTGAtc                                           <  6:209786/147‑1 (MQ=255)
                                                                                                          cGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGAtcgtcg                                       <  3:53467/149‑1 (MQ=255)
                                                                                                            aCCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGAtcgtcg                                       >  2:107104/1‑147 (MQ=255)
                                                                                                            aCCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGaa                                     <  4:84801/149‑1 (MQ=255)
                                                                                                            aCCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGaa                                     <  6:29368/149‑1 (MQ=255)
                                                                                                            aCCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGaa                                     <  6:305423/149‑1 (MQ=255)
                                                                                                              ccGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGa                                   <  6:90677/149‑1 (MQ=255)
                                                                                                                  aCCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGaa                               <  4:352476/149‑1 (MQ=255)
                                                                                                                        aggCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGc                         >  2:14479/2‑149 (MQ=255)
                                                                                                                         ggCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCg                        <  6:358210/149‑1 (MQ=255)
                                                                                                                            cGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGa                     >  6:103944/1‑149 (MQ=255)
                                                                                                                             gAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCCGCGGGAt                    >  7:212/1‑149 (MQ=255)
                                                                                                                             gAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGAt                    <  6:286014/149‑1 (MQ=255)
                                                                                                                              aaGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGAtt                   >  6:46715/1‑149 (MQ=255)
                                                                                                                                   gttgatcAGCACCTGGGCCAGGCGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCg                >  1:246289/7‑147 (MQ=255)
                                                                                                                                   gACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCg                >  7:312723/1‑147 (MQ=255)
                                                                                                                                      ccGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCgaagaa           >  7:147629/1‑149 (MQ=255)
                                                                                                                                        gCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAAca         <  6:77717/149‑1 (MQ=255)
                                                                                                                                        gCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAAca         >  5:161237/1‑149 (MQ=255)
                                                                                                                                        gCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAAc          >  3:300374/1‑148 (MQ=255)
                                                                                                                                         cAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAAcat        <  2:78823/149‑1 (MQ=255)
                                                                                                                                          aGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAAcatc       <  5:103944/149‑1 (MQ=255)
                                                                                                                                          aGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAAcat        >  5:239172/1‑148 (MQ=255)
                                                                                                                                             gcatGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAAcatcat     >  1:298893/4‑149 (MQ=255)
                                                                                                                                             gcCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAACATCATg    <  6:22468/149‑1 (MQ=255)
                                                                                                                                              ccTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGCAGAAGGCAGCGGGATTCCCAAGAAaatcagg    >  4:10215/1‑146 (MQ=255)
                                                                                                                                               cTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAACATCATGGa  <  3:262278/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ACCAGAACAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAACATCATGGA  >  NC_002947/5336038‑5336329

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: