Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,989,452 G→A 18.3% intergenic (+43/+46) PP_5242 → / ← ycaC‑II GAF domain‑containing sensor histidine kinase/putative hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,989,4520GA18.3% 415.7 / 80.7 187intergenic (+43/+46)PP_5242/ycaC‑IIGAF domain‑containing sensor histidine kinase/putative hydrolase
Reads supporting (aligned to +/- strand):  ref base G (77/75);  new base A (13/21);  total (90/97)
Fisher's exact test for biased strand distribution p-value = 2.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-01

Only 100 of 193 total aligned reads displayed.

CGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGT  >  NC_002947/5989303‑5989598
                                                                                                                                                      |                                                                                                                                                  
cGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCgg                                                                                                                                                    <  1:406119/150‑1 (MQ=255)
cGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCgg                                                                                                                                                    <  1:1423288/150‑1 (MQ=255)
cGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCgg                                                                                                                                                    <  1:1424671/150‑1 (MQ=255)
     tGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTa                                                                                                                                               <  1:2958236/150‑1 (MQ=255)
      gTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAg                                                                                                                                              <  2:2318190/150‑1 (MQ=255)
       tAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTTAACCCGCTGCTACGGTCGGGTGTAgg                                                                                                                                             >  2:3016353/1‑150 (MQ=255)
       tAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCCCTCCTACGGTCGGGTGTAgg                                                                                                                                             >  1:4036545/1‑150 (MQ=255)
       tAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTCCTACGGTCGGCGGTAgg                                                                                                                                             >  1:4642581/1‑150 (MQ=255)
        aGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAggg                                                                                                                                            <  1:1050483/150‑1 (MQ=255)
         gTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCTGGAgggg                                                                                                                                           >  2:1966665/1‑150 (MQ=255)
          ttCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGc                                                                                                                                          <  1:4647829/150‑1 (MQ=255)
          ttCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGc                                                                                                                                          <  2:5076002/150‑1 (MQ=255)
          ttCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGc                                                                                                                                          <  2:3468909/150‑1 (MQ=38)
           tCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCg                                                                                                                                         <  2:4961962/150‑1 (MQ=38)
            cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGCCGACCGTAGGAGCgg                                                                                                                                        <  2:305781/150‑1 (MQ=38)
            cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGCCGACCGTAGGAGCgg                                                                                                                                        <  1:837310/150‑1 (MQ=38)
            cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGCCGACCGTAGGAGCgg                                                                                                                                        <  2:5017417/150‑1 (MQ=38)
            cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCgg                                                                                                                                        <  1:3352565/150‑1 (MQ=38)
            cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCgg                                                                                                                                        <  2:1809810/150‑1 (MQ=38)
            cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTACTGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCgg                                                                                                                                        <  1:5173015/150‑1 (MQ=37)
             ggTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCggg                                                                                                                                       <  1:1933295/150‑1 (MQ=255)
              gTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGt                                                                                                                                      <  2:356782/150‑1 (MQ=255)
               tGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGtt                                                                                                                                     <  2:2237224/150‑1 (MQ=255)
                   tCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGGAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTGGCGGGTAAAACCCCACCTACGGTAGGCGGGAGGGGGGGGGTTAcc                                                                                                                                 >  1:3430168/1‑150 (MQ=255)
                   tCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGGGGTGGGGGCGGGTTTAcc                                                                                                                                 >  1:4974126/1‑150 (MQ=255)
                   tCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTAcc                                                                                                                                 <  1:1509604/150‑1 (MQ=255)
                   tCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTAcc                                                                                                                                 >  2:5125793/1‑150 (MQ=255)
                   tCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGGTTAcc                                                                                                                                 >  1:5597364/1‑150 (MQ=255)
                       gCCCACGGCGGGCAGGTGCTGGCGCTGGTGTCGCTTTTGGGTGTATGCCTGATTCTTGTCAGTCTGCATTTTAATCACGCATAGCCTGTGCCGGCCTTTTCGCGGGTAAACCCGCTCCTCCAGCCGACCGTAGGAGCGGTTTTACCcgcg                                                                                                                             <  1:4575739/150‑1 (MQ=25)
                         ccACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTTAACCCCCTCCTACGGGCGGCGGGAGGGGCGGGTTTACCcgcgga                                                                                                                           >  2:4693111/1‑148 (MQ=39)
                           aCGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAaga                                                                                                                         <  2:1624312/150‑1 (MQ=25)
                            cggcggGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAagag                                                                                                                        >  2:5145638/1‑150 (MQ=255)
                              gcggGCTGGTGCAGGCGGTGGTGTCGCTGTTGGGTGGTTGCGGGATTGGTGTCAGTCTGCCTTTTAATCACGCTTAGCCTGTGCCGGCCTGTTCGCGGTTAAACCCGCTCCTACGGTCGGCTTTTGGGGGGGGTTTACCCGCGAAGAGGc                                                                                                                      <  2:2944922/150‑1 (MQ=255)
                                    aGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCggg                                                                                                                                       >  2:5504811/1‑127 (MQ=37)
                                    aGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCggg                                                                                                                                       <  1:5504811/127‑1 (MQ=37)
                                       tGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCGACGGTCGGCTGTAGGGGCGGGGTTACCCGCGAAGAGGCCAGTCCAgg                                                                                                             >  2:1065825/1‑150 (MQ=255)
                                         cAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTAACAGCGAAGTGGCCAGTACAGGtt                                                                                                           >  2:2821758/1‑150 (MQ=255)
                                          aGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTAATCGCCGGTAAAACCCCTCAAACCGTCGGGGTGTCGGGCGGGTTTTCCCCCCAAGAGGGCACTACAAGttt                                                                                                          >  1:3855087/1‑150 (MQ=255)
                                                  gtgtCGCCGGTGGGTGGAGGCCGGATTCGTGTCCGTCGGCCTTTGAAGCCCGCCTCGCCTGTGCGGGCCTCGGCGCGGGTAACCCCGCTCGTCCGGTCGGCGGTAGGGGGGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAgc                                                                                                  <  2:4747245/150‑1 (MQ=255)
                                                  gtgtCGCAGGGGGGTGGCTGCCGGATTCGTGTCAGTCGGCCTTTGAATCCCGCAGCGCGGGGGCCGGCCTCTGCGCGGGTACACCCGCTCCTCCGGGGGGCGGTCGGGGGGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAgc                                                                                                  <  2:4747259/150‑1 (MQ=255)
                                                    gtCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTCCGCGGGTAAACCCGCTCCTACGGTCGGATGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAgcgc                                                                                                <  1:2317373/150‑1 (MQ=255)
                                                     tCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGcgcg                                                                                               >  2:3642417/1‑150 (MQ=255)
                                                      cGCTGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGa                                                                                              <  2:2573179/150‑1 (MQ=255)
                                                             gggTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCCCTCCTACCGGCGGCTGGAGGGGCGGGTTTACCCGCGAAGAGGCCCGTACAGGTTTAAAAAAGCGCGAGGCCAgc                                                                                       >  1:4099053/1‑150 (MQ=37)
                                                              ggTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTAGCGGGTTAAACCCCTCCAAACGGGGGCGGTTTGGGGGGGGTTAAACCCCAACAGGCCCCTACAGGTTTAAAAAATCCCGAGGCCaccc                                                                                      >  1:4857042/1‑147 (MQ=255)
                                                               gTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTAACCGCGAAGAGGCCAGTAAAGGTTTACATAAGCGCGAGGCCAgcgc                                                                                     >  2:4399974/1‑150 (MQ=255)
                                                                      ccGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTCCGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAg                                                                              <  2:5372565/150‑1 (MQ=255)
                                                                      ccGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAg                                                                              <  2:4118428/150‑1 (MQ=255)
                                                                       cGGATTCGTGTCAGTCTGCATTTGAATCACGCATAGCCTGAGCCGTACTCTTCGCGGGTAAAACCCCGCCTACGGGAGGGGGTAGGGGCGGGTTTAACCGCCCAGAGGCCACTAAAGGGTTACAAAAGCCCGAGGCCAGCGCTGCCTacc                                                                             >  2:3523711/1‑148 (MQ=255)
                                                                              gtgtCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCATAGCGGGTAAAACCCCTCCAAAGGTCGGCTGTAGGGGCGGGTTTAACCCCCAAGAGGCCAGTAAACGTTTAAAAAAGCGCGAGGCCACCGCTGCGTAGCGGTtcaa                                                                      >  1:4078730/1‑149 (MQ=255)
                                                                                gtctgtctGCCTTTCAACCCCGCCTCGCCTGTGCGGGCCTCTCCGCGGGAAAACCCGCTCCTCCGGTGGGGTGTGGGGGGGGTTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTtcatca                                                                    <  1:1230213/150‑1 (MQ=255)
                                                                                gtCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTtcatca                                                                    <  2:1749797/150‑1 (MQ=255)
                                                                                      cTGCCTTTGAATCACGCATCGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCGGTCGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGc                                                              <  1:4611382/150‑1 (MQ=255)
                                                                                       ggCCTTTCAAACACGCATCGCCTGTGCCGGCCTCTGCGGGGGTAAACCCGCTCCTCCGGTCGGCTTTCGGGGGGGTTTCACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCg                                                             <  1:312247/149‑1 (MQ=255)
                                                                                           tttGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCGGGAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTTCGTAGCTGTTCATCAGGGTGCGGTAg                                                         >  2:1228351/1‑150 (MQ=255)
                                                                                               aaTCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTgg                                                     <  2:4514131/150‑1 (MQ=255)
                                                                                                 tCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCCCTCCGACCGTCGGCGGTGGGGGCGGGTTTTCCCGCGAACTGGCCCGTACAGGTTTAAAACAGCGCGAGGCCAGCGCGGCGTAGCTGTTAATCACGGTGCGGGGGTTgggg                                                   >  2:984771/1‑150 (MQ=255)
                                                                                                   ccGCTTAGCCTGTGCGGGCCTCTTCGGGGGTAAACCCGCTCCTACGTTGGGCTGTTGGGGGGGTTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGAt                                                 <  1:1686029/149‑1 (MQ=38)
                                                                                                   aCGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCAACGGTCGGCTGTTGGGGCGGGTTTACCCGCGACGAGGCCACTACAGGGTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGGAGTTGGGGAt                                                 >  2:3543008/1‑150 (MQ=38)
                                                                                                      cATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCGACGGGCGGCGGGAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACCGGTTTACATCAGCGCGAGGCCAGCGCTGCGGAGCTGTTCAGCAGGTTGCGGGAGTTGGGGATACg                                              >  2:3166918/1‑150 (MQ=37)
                                                                                                              gggCCGGCCTCTTGGGGGTTAACCCCGCTCCTCGGGGGGGCTGTTGGGGGGGGTTTCCCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGAca                                      <  1:2852731/148‑1 (MQ=255)
                                                                                                                gCCGGCCTCTTCGCGGGTAAACCCGCTCCTCCGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGAcagc                                    <  1:3730249/150‑1 (MQ=255)
                                                                                                                gCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTCGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGAcagc                                    <  2:1135605/150‑1 (MQ=255)
                                                                                                                   ggCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTCGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAgg                                 <  2:641960/150‑1 (MQ=255)
                                                                                                                    gCCTCTTCGCGGGTAACCCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGt                                <  2:3225061/150‑1 (MQ=255)
                                                                                                                        aggagcgGGTTTACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTAcc                            >  2:3033068/5‑150 (MQ=255)
                                                                                                                         ttCGCGGGTAAACCCGCTCCTACGGTCGGCGGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTAccc                           >  2:2464249/1‑150 (MQ=255)
                                                                                                                            gcgGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGc                        >  2:4750474/1‑150 (MQ=255)
                                                                                                                            gcgGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGc                        >  2:747882/1‑150 (MQ=37)
                                                                                                                               gctAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGACGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCCGCGCTGCGTAGCTGTGCATCAGGTTGCGGTCGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCt                     >  2:2577380/3‑150 (MQ=255)
                                                                                                                                gTAAACCCGCTCCTCCGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCtt                    <  2:681762/150‑1 (MQ=255)
                                                                                                                                 tacaCCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTc                   <  2:5455776/147‑1 (MQ=255)
                                                                                                                                 tAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGTTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTc                   <  2:833874/150‑1 (MQ=255)
                                                                                                                                    aCCCGCTCCTCCGGTCGGGTGTCGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGAt                <  1:4454237/150‑1 (MQ=255)
                                                                                                                                     cccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATg               >  1:179201/1‑150 (MQ=255)
                                                                                                                                     cccGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATg               >  2:3886846/1‑150 (MQ=38)
                                                                                                                                     cccGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATg               >  1:5300834/1‑150 (MQ=38)
                                                                                                                                      ccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt              >  2:379786/1‑150 (MQ=255)
                                                                                                                                      ccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt              >  1:2499737/1‑150 (MQ=255)
                                                                                                                                      ccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt              >  1:732776/1‑150 (MQ=255)
                                                                                                                                      ccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt              >  1:3259651/1‑150 (MQ=255)
                                                                                                                                      ccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt              >  1:2194598/1‑150 (MQ=255)
                                                                                                                                      ccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt              >  2:33642/1‑150 (MQ=255)
                                                                                                                                      ccGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt              >  2:3020240/1‑150 (MQ=38)
                                                                                                                                      ccGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt              >  1:2816970/1‑150 (MQ=38)
                                                                                                                                       cGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTc             >  2:759260/1‑150 (MQ=255)
                                                                                                                                       cGCCCCTCCGGTCGGCTGTAGGGGGGGTTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTc             <  1:111573/150‑1 (MQ=255)
                                                                                                                                        gCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCg            >  2:4826804/1‑150 (MQ=255)
                                                                                                                                            ctccGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTgc        <  2:2114978/147‑1 (MQ=255)
                                                                                                                                            cTACGGTCGGCTGTTGGGGCGTGTTTACACTCGAAGAGTCCAGTAGAGGTTTAGATCAGCTCGAGGCCAGCGCTGCGTAGCTGTTCATCAGTTTGCTGTAGTTGGGGATACGCTGCGTCATCAGTTTACCCAGCCCTTCTATGTCGTTgc        <  2:2122912/150‑1 (MQ=255)
                                                                                                                                            cTACGGTCGGCTGTTGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTgc        <  2:4311822/150‑1 (MQ=255)
                                                                                                                                            cTACGGTCGGCTGTAGGGGCGGTTTAACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTgc        <  2:1192294/150‑1 (MQ=255)
                                                                                                                                            cTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTgc        >  2:4977482/1‑150 (MQ=255)
                                                                                                                                            cTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTgc        >  2:4976669/1‑150 (MQ=255)
                                                                                                                                               cGGTCGGCTGTAGGGGCGGTTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGcc     <  1:5018118/150‑1 (MQ=255)
                                                                                                                                                 gTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGACAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAg   >  2:1338621/1‑150 (MQ=255)
                                                                                                                                                 gTCGGCTGTAGGGGCGGGTTAACCCGCGAAGAGGCCAGTACAGGTTTCCATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGGAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGACGTTGCGCCAg   >  2:3194642/1‑150 (MQ=255)
                                                                                                                                                  tCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGt  >  1:3559431/1‑150 (MQ=255)
                                                                                                                                                      gCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAg                                                                              >  2:2767471/1‑71 (MQ=255)
                                                                                                                                                      gCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAg                                                                              <  1:2767471/71‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
CGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGT  >  NC_002947/5989303‑5989598

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