Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,291,680 G→C 5.6% A350G (GCC→GGC)  PP_3758 ← two component system response regulator/sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,291,6800GC 5.6% 912.4 / 31.4 269A350G (GCC→GGC) PP_3758two component system response regulator/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (139/115);  new base C (5/10);  total (144/125)
Fisher's exact test for biased strand distribution p-value = 1.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 282 total aligned reads displayed.

TTTGTTCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTGCGCTCGCCGTAACGCATTGCATTG  >  NC_002947/4291531‑4291829
                                                                                                                                                     |                                                                                                                                                     
   gTTCTTTGAATGGTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAg                                                                                                                                                    <  2:3193899/150‑1 (MQ=255)
    ttCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGATGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGc                                                                                                                                                   <  1:3105626/150‑1 (MQ=255)
     tCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCg                                                                                                                                                  <  2:2904491/150‑1 (MQ=255)
     tCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGCCCAGCAGGCCAGCg                                                                                                                                                  >  2:5064975/1‑150 (MQ=255)
       tttGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGAc                                                                                                                                                >  1:1068518/1‑150 (MQ=255)
       tttGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGAc                                                                                                                                                >  2:5528162/1‑150 (MQ=255)
       tttGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGAc                                                                                                                                                >  1:6777485/1‑150 (MQ=255)
       tttGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGAc                                                                                                                                                <  1:3618496/150‑1 (MQ=255)
       tttGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGCCCAGCAGGCCAGCgga                                                                                                                                                <  2:5215169/150‑3 (MQ=255)
       tttGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGCCCAGCAGGCCAGCgga                                                                                                                                                <  1:1856547/150‑3 (MQ=255)
       tttGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGCCCAGCAGGCCAGCgga                                                                                                                                                <  1:941145/150‑3 (MQ=255)
       tttGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGCCCAGCAGGCCAGCgga                                                                                                                                                <  1:5153188/150‑3 (MQ=255)
       tttGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCGCTGCCCAGCAGGCCAGCggag                                                                                                                                               <  2:2346230/150‑4 (MQ=255)
       tttGAATGGTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTGGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGCCCAGCAGGCCAGCgga                                                                                                                                                <  2:5930253/150‑3 (MQ=255)
        ttGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGCCCAGCAGGCCAGCggag                                                                                                                                               <  2:1895094/150‑4 (MQ=255)
         tGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGAccc                                                                                                                                              >  1:2495053/1‑150 (MQ=255)
          gAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCg                                                                                                                                             >  2:3594849/1‑150 (MQ=255)
           aaTGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGc                                                                                                                                            <  2:4911336/150‑1 (MQ=255)
            aTGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCt                                                                                                                                           >  1:5714997/1‑150 (MQ=255)
                 gCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAAc                                                                                                                                      <  2:6874527/150‑1 (MQ=255)
                  ccAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAAct                                                                                                                                     >  2:1942137/1‑150 (MQ=255)
                   cAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAActg                                                                                                                                    >  2:3678804/1‑150 (MQ=255)
                    aGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAActgc                                                                                                                                   <  2:6227524/150‑1 (MQ=255)
                    aGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAActgc                                                                                                                                   >  2:4822046/1‑150 (MQ=255)
                         cAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGaa                                                                                                                              <  1:3451369/150‑1 (MQ=255)
                               cAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCg                                                                                                                        >  1:1036007/1‑150 (MQ=255)
                                aGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCgc                                                                                                                       <  2:91306/150‑1 (MQ=255)
                                  ggTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCgctg                                                                                                                     <  2:2412399/150‑1 (MQ=255)
                                  ggTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCgctg                                                                                                                     >  1:6310402/1‑150 (MQ=255)
                                   gTGAATGTGGTACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCgctgc                                                                                                                    >  2:3121848/1‑150 (MQ=255)
                                   gTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCgctgc                                                                                                                    >  1:558435/1‑150 (MQ=255)
                                   gTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCgctgc                                                                                                                    >  1:2039771/1‑150 (MQ=255)
                                   gTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCgctgc                                                                                                                    <  1:423006/150‑1 (MQ=255)
                                       aTGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGt                                                                                                                <  2:2795528/150‑1 (MQ=255)
                                       aTGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGt                                                                                                                <  1:3765525/150‑1 (MQ=255)
                                         gtgGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTg                                                                                                              <  1:4154458/150‑1 (MQ=255)
                                          tgGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTgg                                                                                                             >  1:4000927/1‑150 (MQ=255)
                                            gCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTggcg                                                                                                           >  1:919531/1‑150 (MQ=255)
                                                ccTTTCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTTCATCCCCAAGCCCAGTCCTCCGCTGGCCTGCTGGGCAGCGTCCCGCTCGAACTGCTTGAAGATCCGCTGCTGGTTggcggctc                                                                                                       <  1:6452636/150‑3 (MQ=255)
                                                ccTAGCCCTCATCGATTACAACCAGTATACGGCCACCTAGAGCCTGTAAGATATGCTCTGAGATGTACATAACCAAGAACAGACCTCCCCTGGCCTGCTGGGCAGCGACACGATCGAACAGCTGGAAGATCCGCTGCTGGTTGGCGgcgc                                                                                                       >  2:4564824/1‑150 (MQ=255)
                                                  ttGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATc                                                                                                                          >  2:5808176/1‑129 (MQ=255)
                                                  ttGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATc                                                                                                                          <  1:5808176/129‑1 (MQ=255)
                                                   tGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCGCCGCTGGCCTGCTGGGCAGCGACCCGATCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTga                                                                                                    <  2:4189930/150‑1 (MQ=255)
                                                    gCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGgcggcg                                                                                                        <  1:4630145/145‑1 (MQ=255)
                                                      cctcctCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTgatac                                                                                                 >  2:2150311/1‑150 (MQ=255)
                                                      cctcctCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTgatac                                                                                                 >  1:4044546/1‑150 (MQ=255)
                                                        tcctcGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTgataccg                                                                                               >  1:3403961/1‑150 (MQ=255)
                                                           tcGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTgataccgata                                                                                            >  1:2393490/1‑150 (MQ=255)
                                                           tcGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTgataccgata                                                                                            >  2:4570769/1‑150 (MQ=255)
                                                               tttCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCTGCTGCCCAGCAGGCCAg                                                                                                                                                    <  1:2344061/90‑1 (MQ=255)
                                                                 tCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGa                                                                                      <  1:5389607/150‑1 (MQ=255)
                                                                   aaCCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATc                                                                                    <  1:1218142/150‑1 (MQ=255)
                                                                    aCCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATcc                                                                                   >  2:202892/1‑150 (MQ=255)
                                                                         cgtCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATcctgcac                                                                              <  1:1560314/148‑1 (MQ=255)
                                                                           tCCGGCCACCGTGCGCCTGGTCGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTTGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATcctgcacct                                                                            <  2:2299690/150‑1 (MQ=255)
                                                                           tCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATcctgcacct                                                                            >  2:4344431/1‑150 (MQ=255)
                                                                           tCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATcctgcacct                                                                            >  2:3008934/1‑150 (MQ=255)
                                                                                  aCCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCa                                                                     >  1:1793064/1‑150 (MQ=255)
                                                                                   ccGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAg                                                                    <  1:3260447/150‑1 (MQ=255)
                                                                                   ccGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAg                                                                    <  2:5588246/150‑1 (MQ=255)
                                                                                   ccGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAg                                                                    <  2:341006/150‑1 (MQ=255)
                                                                                   ccGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAg                                                                    <  2:5612770/150‑1 (MQ=255)
                                                                                         gcCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATg                                                              >  1:5599506/1‑150 (MQ=255)
                                                                                             ggACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATgccgc                                                          <  2:2562696/150‑1 (MQ=255)
                                                                                             ggACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATgccgc                                                          <  2:5293288/150‑1 (MQ=255)
                                                                                              gACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGgcgc                                                                                                       <  2:4055517/104‑1 (MQ=255)
                                                                                              gACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATgccgcc                                                         <  2:6659009/150‑1 (MQ=255)
                                                                                              gACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATgccgcc                                                         <  2:39807/150‑1 (MQ=255)
                                                                                                 gATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTg                                                      >  1:2796523/1‑150 (MQ=255)
                                                                                                  aTCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGc                                                     <  1:4460597/150‑1 (MQ=255)
                                                                                                  aTCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGc                                                     <  1:1674411/150‑1 (MQ=255)
                                                                                                        tCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAac                                               >  1:2110837/1‑150 (MQ=255)
                                                                                                        tCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAac                                               >  2:161321/1‑150 (MQ=255)
                                                                                                           gagaTGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACAcc                                            <  2:809454/150‑1 (MQ=255)
                                                                                                            agaTGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACAccc                                           <  2:1754245/150‑1 (MQ=255)
                                                                                                             gaTGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCg                                          <  2:2369655/150‑1 (MQ=255)
                                                                                                                   cAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTg                                    >  1:751175/1‑150 (MQ=255)
                                                                                                                      ccccAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCAc                                 <  2:3282004/150‑1 (MQ=255)
                                                                                                                      ccccAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCAc                                 >  2:3948026/1‑150 (MQ=255)
                                                                                                                        ccAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCg                               >  2:4746455/1‑150 (MQ=255)
                                                                                                                        ccAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCg                               >  1:4381477/1‑150 (MQ=255)
                                                                                                                         cAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCgg                              >  2:1487591/1‑150 (MQ=255)
                                                                                                                              ccAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTgc                         <  2:4431296/150‑1 (MQ=255)
                                                                                                                               cAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTgcg                        >  2:1078866/1‑150 (MQ=255)
                                                                                                                               cAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTgcg                        >  2:6425007/1‑150 (MQ=255)
                                                                                                                                aGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTgcgc                       >  2:5038538/1‑150 (MQ=255)
                                                                                                                                aGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTgcgc                       >  1:6317300/1‑150 (MQ=255)
                                                                                                                                    cTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCGCCTGCACCGGATTGCGCTCGc                   >  1:3808813/1‑150 (MQ=255)
                                                                                                                                    cTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTGCGCTCGc                   >  1:597669/1‑150 (MQ=255)
                                                                                                                                     tCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTGCGCTCGcc                  >  1:6827954/1‑150 (MQ=255)
                                                                                                                                      tcGCTGCCCAGCAGGCCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTGCGCTCGCCg                 >  1:2712476/2‑150 (MQ=255)
                                                                                                                                      ccGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTGCGCTCGCCg                 >  2:5444404/1‑150 (MQ=255)
                                                                                                                                               tgctgGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTGCGCTCGCCGTAACgcatt        <  2:2704987/150‑1 (MQ=255)
                                                                                                                                               tgctgGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTGCGCTCGCCGTAACgcatt        >  1:4562077/1‑150 (MQ=255)
                                                                                                                                               tgctgGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTGCGCTCGCCGTAACgcatt        <  2:4183290/150‑1 (MQ=255)
                                                                                                                                               tgctgGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTGCGCTCGCCGTAACgcatt        <  1:3050576/150‑1 (MQ=255)
                                                                                                                                                    ggCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTGCGCTCGCCGTAACgcattgcatt   <  2:6317300/150‑1 (MQ=255)
                                                                                                                                                    ggCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTGCGCTCGCCGTAACgcattgcatt   <  2:6688469/150‑1 (MQ=255)
                                                                                                                                                     gCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTGCGCTCGCCGTAACGcattgcattg  >  2:5363828/1‑150 (MQ=255)
                                                                                                                                                     gCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTGCGCTCGCCGTAACGcattgcattg  >  1:5062620/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTGTTCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTGCGCTCGCCGTAACGCATTGCATTG  >  NC_002947/4291531‑4291829

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