Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,868,674 G→T 5.2% intergenic (+19/+38) PP_5140 → / ← thyA MerR family transcriptional regulator/thymidylate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,868,6740GT 5.2% 728.5 / 20.9 210intergenic (+19/+38)PP_5140/thyAMerR family transcriptional regulator/thymidylate synthase
Reads supporting (aligned to +/- strand):  ref base G (106/93);  new base T (9/2);  total (115/95)
Fisher's exact test for biased strand distribution p-value = 1.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

Only 100 of 214 total aligned reads displayed.

TGGTGGAGCTGCGGCGGCGCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCAGTTTCGGCA  >  NC_002947/5868525‑5868823
                                                                                                                                                     |                                                                                                                                                     
tggtggAGCTGCGGCGGCGCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTg                                                                                                                                                       >  2:338662/1‑150 (MQ=255)
 ggtggAGCTGCGGCGGCGCTGCAACGCGCAGGGGAGTGACTGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACCTTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGa                                                                                                                                                      >  2:712183/1‑150 (MQ=255)
   tggAGCTGCGGCGGCGCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAAc                                                                                                                                                    <  1:3198075/150‑1 (MQ=255)
      aGCTGCGGCGGCGCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCcg                                                                                                                                                 <  1:1305863/150‑1 (MQ=255)
         tGCGGCGGCGCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGa                                                                                                                                              >  1:2406093/1‑150 (MQ=255)
         tGCGGCGGCGCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGa                                                                                                                                              <  1:3901651/150‑1 (MQ=255)
          gcggcggcgCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTTCACCCGCGaa                                                                                                                                             <  2:266406/150‑1 (MQ=255)
          gcggcggcgCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGaa                                                                                                                                             <  2:2921842/150‑1 (MQ=255)
          gcggcggcgCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGaa                                                                                                                                             >  1:654087/1‑150 (MQ=255)
          gcggcggcgCTGCAACGCGCAGCGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTAGCGGGTGAACCCGCGaa                                                                                                                                             >  2:1498050/1‑150 (MQ=255)
           cggcggcgCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAag                                                                                                                                            >  2:1996744/1‑150 (MQ=255)
            ggcggcgCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAaga                                                                                                                                           <  1:346851/150‑1 (MQ=255)
            ggcggcgCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAaga                                                                                                                                           <  1:3069684/150‑1 (MQ=255)
             gcggcgCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAagag                                                                                                                                          <  1:360764/150‑1 (MQ=255)
             gcggcgCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAagag                                                                                                                                          <  2:3698050/150‑1 (MQ=255)
                 cgcTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCa                                                                                                                                      <  1:4505524/150‑1 (MQ=255)
                 cgcTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCa                                                                                                                                      <  2:2696214/150‑1 (MQ=255)
                               ggggAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATaga                                                                                                                        <  2:1206330/150‑1 (MQ=255)
                               ggggAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATaga                                                                                                                        >  2:400732/1‑150 (MQ=255)
                               ggggAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATaga                                                                                                                        <  2:463556/150‑1 (MQ=255)
                               ggggAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATaga                                                                                                                        <  1:942448/150‑1 (MQ=255)
                               ggggAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATaga                                                                                                                        >  1:3732282/1‑150 (MQ=255)
                                  gAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGAt                                                                                                                     <  2:1052719/150‑1 (MQ=255)
                                    gTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGGGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCt                                                                                                                   >  1:2328569/1‑150 (MQ=255)
                                       aaTGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGGGGTCAGTTCTAACAATAGAGATCTATc                                                                                                                <  1:4241094/150‑1 (MQ=255)
                                        aTGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCa                                                                                                               <  2:1610198/150‑1 (MQ=255)
                                         tGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAg                                                                                                              >  1:4153176/1‑150 (MQ=255)
                                         tGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAg                                                                                                              >  2:2400409/1‑150 (MQ=255)
                                          gcgcGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGa                                                                                                             >  2:734036/1‑150 (MQ=255)
                                          gcgcGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGa                                                                                                             >  2:4429904/1‑150 (MQ=255)
                                           cgcgATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGAc                                                                                                            >  1:3934488/1‑150 (MQ=255)
                                            gcgATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGAcc                                                                                                           >  2:3918912/1‑150 (MQ=255)
                                            gcgATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGAcc                                                                                                           >  2:430557/1‑150 (MQ=255)
                                              gATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGc                                                                                                         >  2:4332670/1‑150 (MQ=255)
                                              gATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGc                                                                                                         >  2:3970919/1‑150 (MQ=255)
                                                tCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCa                                                                                                       <  2:2121754/150‑1 (MQ=255)
                                                  ttGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATc                                                                                                     >  1:3440054/1‑150 (MQ=255)
                                                    gcagcaACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGa                                                                                                                                              <  2:1325933/107‑1 (MQ=255)
                                                      agcaACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGgcg                                                                                                 >  1:3865149/1‑150 (MQ=255)
                                                      agcaACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGgcg                                                                                                 >  1:1353884/1‑150 (MQ=255)
                                                         aaCTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTg                                                                                              >  2:3909160/1‑150 (MQ=255)
                                                         aaCTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTg                                                                                              >  2:81707/1‑150 (MQ=255)
                                                                 aCAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGTCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGgc                                                                                      <  2:2716149/150‑1 (MQ=255)
                                                                 aCAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGgc                                                                                      <  2:1177201/150‑1 (MQ=255)
                                                                   aaaCGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGgcgc                                                                                    >  1:51672/1‑150 (MQ=255)
                                                                   aaaCGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGACATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGgcgc                                                                                    >  2:1917712/1‑150 (MQ=255)
                                                                    aaCGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCa                                                                                   >  2:4164666/1‑150 (MQ=255)
                                                                     aCGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCAt                                                                                  <  1:430162/150‑1 (MQ=255)
                                                                     aCGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCAt                                                                                  >  2:752879/1‑150 (MQ=255)
                                                                              tATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCAt                                                                         <  1:3020530/150‑1 (MQ=255)
                                                                                tCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAg                                                                       >  1:585615/1‑150 (MQ=255)
                                                                                    tACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCtt                                                                   >  2:1435975/1‑150 (MQ=255)
                                                                                       cGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTcca                                                                <  1:4370907/150‑1 (MQ=255)
                                                                                       cGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTcca                                                                <  1:3522565/150‑1 (MQ=255)
                                                                                        ggATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTccag                                                               >  2:3583193/1‑150 (MQ=255)
                                                                                        ggATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTccag                                                               >  2:1955459/1‑150 (MQ=255)
                                                                                              cGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGaa                                                         >  2:4144176/1‑150 (MQ=255)
                                                                                              cGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGaa                                                         >  2:2505572/1‑150 (MQ=255)
                                                                                              cGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGaa                                                         <  1:1029276/150‑1 (MQ=255)
                                                                                              cGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGaa                                                         >  2:1887388/1‑150 (MQ=255)
                                                                                                   aTTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCAc                                                    >  1:2936762/1‑150 (MQ=255)
                                                                                                   aTTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCAc                                                    <  2:3866133/150‑1 (MQ=255)
                                                                                                     tCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTc                                                  <  1:1896368/150‑1 (MQ=255)
                                                                                                      cccATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCg                                                 <  1:1981273/150‑1 (MQ=255)
                                                                                                              gggCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAAtc                                         <  2:926551/150‑1 (MQ=255)
                                                                                                                gCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGACATCACTCGGCTCAAtctt                                       >  2:3193589/1‑150 (MQ=255)
                                                                                                                 cGGTGTCACGGGCAGTGAGGCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAAtcttc                                      <  1:3055538/150‑1 (MQ=255)
                                                                                                                   gAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTtctc                                    >  1:1921068/1‑150 (MQ=255)
                                                                                                                   gAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTtctc                                    >  1:1686838/1‑150 (MQ=255)
                                                                                                                        aCGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGcc                               >  2:3011638/1‑150 (MQ=255)
                                                                                                                        aCGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGcc                               >  2:2135589/1‑150 (MQ=255)
                                                                                                                        aCGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGcc                               >  2:2995748/1‑150 (MQ=255)
                                                                                                                        aCGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGcc                               >  1:4204495/1‑150 (MQ=255)
                                                                                                                        aCGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGcc                               <  2:903178/150‑1 (MQ=255)
                                                                                                                            gCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTc                           >  2:2774540/1‑150 (MQ=255)
                                                                                                                               ttGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCgg                        <  1:4332670/150‑1 (MQ=255)
                                                                                                                                tGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGc                       <  1:609846/150‑1 (MQ=255)
                                                                                                                                tGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGc                       <  1:3258423/150‑1 (MQ=255)
                                                                                                                                tGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGc                       <  2:3838162/150‑1 (MQ=255)
                                                                                                                                  aGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTc                     >  2:2205196/1‑150 (MQ=255)
                                                                                                                                     cAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATa                  <  1:124419/150‑1 (MQ=255)
                                                                                                                                     cAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATa                  >  2:2581782/1‑150 (MQ=255)
                                                                                                                                      aaaaGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATAc                 >  2:201283/1‑150 (MQ=255)
                                                                                                                                      aaaaGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATAc                 >  2:1822026/1‑150 (MQ=255)
                                                                                                                                      aaaaGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATAc                 <  1:741088/150‑1 (MQ=255)
                                                                                                                                      aaaaGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATAc                 >  2:346707/1‑150 (MQ=255)
                                                                                                                                          tcTTCGCGGGTTCACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGc             >  1:2675856/2‑150 (MQ=255)
                                                                                                                                          tcTTCGCGGGTTCACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGc             >  1:1938470/2‑150 (MQ=255)
                                                                                                                                          tcTTCGCGGGTTCACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGc             >  2:2641818/2‑150 (MQ=255)
                                                                                                                                          tcTTCGCGGGTTCACCAGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCAAAGGAGCATGGTGCTCATAGACTTCCAGCCAGAAATCACTCGGCTCAAACTTCTCCAGCCACTCCGGCTCATACTTGc             >  1:4103631/2‑150 (MQ=255)
                                                                                                                                          gCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGc             <  2:661834/150‑1 (MQ=255)
                                                                                                                                            ttCGCGGGTTCACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCCCATACTTGCCa           >  2:4289271/1‑150 (MQ=255)
                                                                                                                                            ttCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCa           >  1:4413684/1‑150 (MQ=255)
                                                                                                                                               gcgGGTTCACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCAGtt        >  1:2696337/1‑150 (MQ=255)
                                                                                                                                               gcgGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCAGtt        >  1:1092043/1‑150 (MQ=255)
                                                                                                                                                cgGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCAGttt       <  1:338662/150‑1 (MQ=255)
                                                                                                                                                   gTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCAGTTTCgg    <  1:1087304/150‑1 (MQ=255)
                                                                                                                                                   gTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCAGTTTCgg    <  1:1052123/150‑1 (MQ=255)
                                                                                                                                                   gTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCAGTTTCgg    <  1:1007734/150‑1 (MQ=255)
                                                                                                                                                   gTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCAGTTTCgg    >  2:2612467/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGTGGAGCTGCGGCGGCGCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCAGTTTCGGCA  >  NC_002947/5868525‑5868823

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: