Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,866,753 A→C 6.4% intergenic (+20/+37) 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,866,7530AC 6.4% 192.5 / 2.5 73intergenic (+20/+37)PP_5140/thyAMerR family transcriptional regulator/thymidylate synthase
Reads supporting (aligned to +/- strand):  ref base A (38/28);  new base C (3/2);  total (42/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGAGCTGCGGCGGCGCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCAGTTTCGG  >  NC_002947/5866605‑5866899
                                                                                                                                                    |                                                                                                                                                   
gtggAGCTGCGGCGGCGCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGaa                                                                                                                                                    <  1:284532/150‑1 (MQ=255)
gtggAGCTGCGGCGGCGCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGaa                                                                                                                                                    <  2:357300/150‑1 (MQ=255)
 tggAGCTGCGGCGGCTCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGCATACCGAACATTGCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAAc                                                                                                                                                   <  2:1051728/150‑1 (MQ=255)
 tggAGCTGCGGCGGCGCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAAc                                                                                                                                                   <  1:1114296/150‑1 (MQ=255)
   gAGCTGCGGCGGCACTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAAccc                                                                                                                                                 >  1:62086/1‑150 (MQ=255)
        gcggcggcgCTGCAACTCGCAGGGGGGTGACCGCGCGATCTTACAGCAAGTCGTGACAAACTGGGCGGGATCGGGACCGGACGCCGAACAGTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGaa                                                                                                                                            <  2:572853/150‑1 (MQ=255)
          ggcggcgCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAaga                                                                                                                                          <  2:522907/150‑1 (MQ=255)
              gcgcTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAACCGGGGCGGCATCGACGCCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTc                                                                                                                                     <  1:797414/150‑1 (MQ=255)
                  tGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGGCAGtt                                                                                                                                 >  2:350180/1‑150 (MQ=255)
                    cAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTAGCGGTACCGGAGACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCt                                                                                                                               >  1:1875281/1‑150 (MQ=255)
                      aCGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTaa                                                                                                                             >  1:1585086/1‑150 (MQ=255)
                            aGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATag                                                                                                                       >  1:620117/1‑150 (MQ=255)
                             ggggAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGAATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATaga                                                                                                                      <  1:1504078/150‑1 (MQ=255)
                                gAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGAt                                                                                                                   >  1:831093/1‑150 (MQ=255)
                                  gTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCt                                                                                                                 >  1:1565915/1‑150 (MQ=255)
                                   tGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTa                                                                                                                >  1:740711/1‑150 (MQ=255)
                                                 tGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCg                                                                                                  >  1:933599/1‑150 (MQ=255)
                                                 tGCAGCAACTGAAGACAAACGGGGCGGTATTGGTACCGGATACCGAATATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACATTAGAGATCTATCAGACCGCCATCg                                                                                                  <  2:142018/150‑1 (MQ=255)
                                                  gcagcaACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCAATGTGCCACAAGCTTCGCCAGCGTTCGTGCCACAGTGTTTTTTCTAACAATAGAGATCTATCAGACCGCCATCgg                                                                                                 >  1:1376072/1‑150 (MQ=255)
                                                     gcaACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGAc                                                                                                          >  2:1351404/1‑138 (MQ=255)
                                                     gcaACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGAc                                                                                                          <  1:1351404/138‑1 (MQ=255)
                                                       aaCTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTg                                                                                            <  2:758041/150‑1 (MQ=255)
                                                         cTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTc                                                                                          <  2:620117/150‑1 (MQ=255)
                                                            gagaCAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCTGGGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTa                                                                                                                              >  1:1786595/1‑111 (MQ=255)
                                                            gagaCAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTa                                                                                                                              <  2:1786595/111‑1 (MQ=255)
                                                            gagaCAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATc                                                                                       <  2:740711/150‑1 (MQ=255)
                                                             agaCAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGGCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCg                                                                                      <  2:1790035/150‑1 (MQ=255)
                                                              gaCAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCgg                                                                                     >  1:402355/1‑150 (MQ=255)
                                                               aCAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGCACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGgc                                                                                    <  1:1430323/150‑1 (MQ=255)
                                                                 aaaCGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGgcgc                                                                                  >  1:926203/1‑150 (MQ=255)
                                                                             aTCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATa                                                                      <  2:1758328/150‑1 (MQ=255)
                                                                                 gTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGCCGCTGCCATCGGCGCATGGTGCTCATTGCCt                                                                  <  2:117455/150‑1 (MQ=255)
                                                                                  tACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGGCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCTGCGCATGGTGCTCATAGCCtt                                                                 >  2:1026846/1‑150 (MQ=255)
                                                                                    ccGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTcc                                                               >  2:1009461/1‑150 (MQ=255)
                                                                                       gATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGCGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTccagc                                                            >  1:68143/1‑150 (MQ=255)
                                                                                       gATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTg                                                                            >  1:1501870/1‑134 (MQ=255)
                                                                                       gATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTg                                                                            <  2:1501870/134‑1 (MQ=255)
                                                                                       gATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTccagc                                                            >  1:398507/1‑150 (MQ=255)
                                                                                            cGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGaa                                                       >  2:300768/1‑150 (MQ=255)
                                                                                             gAACATTCCCATGTGGGGCGGAGTCACTGGCATTGAGCCAAATGCTAAGCGGTCGGCACGGTTAGTGTGGGTGTTCTTACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGaaa                                                      <  1:224087/150‑1 (MQ=255)
                                                                                             gAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGaaa                                                      >  2:806396/1‑150 (MQ=255)
                                                                                             gAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGATGAACCCGCGAAGAGGTCAGGTCTAACAATAGAGATATATCAGACCGCCATCGGCGCTGTCATCGACGCATGGTGCTCATAGCCTTCCAGCCAGaaa                                                      >  2:914045/1‑150 (MQ=255)
                                                                                              aaCATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAAt                                                     <  2:1376072/150‑1 (MQ=255)
                                                                                                 aTTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCAc                                                  <  2:1875281/150‑1 (MQ=255)
                                                                                                     ccATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGaa                                                                                                                                            >  2:586823/1‑57 (MQ=255)
                                                                                                     ccATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGaa                                                                                                                                            <  1:586823/57‑1 (MQ=255)
                                                                                                     ccATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTACGAGGGTGCACCCACGGCCAGGGAAGTAATTTCAATAGTGAGCTATCGGGCGCCCATCTGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCgg                                              >  2:1161822/1‑150 (MQ=255)
                                                                                                      cATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATTAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGc                                             >  2:177392/1‑150 (MQ=255)
                                                                                                      cATGTGGGGCGGAGTCACGCGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGc                                             >  2:1608074/1‑150 (MQ=255)
                                                                                                       aTGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACt                                                 <  2:74777/145‑1 (MQ=255)
                                                                                                       aTGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACt                                                 >  1:74777/1‑145 (MQ=255)
                                                                                                          tgGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCaa                                         >  2:252735/1‑150 (MQ=255)
                                                                                                           ggggCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAAt                                        <  2:1565915/150‑1 (MQ=255)
                                                                                                              gCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCAt                                                                                        <  2:310475/99‑1 (MQ=255)
                                                                                                              gCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCAt                                                                                        >  1:310475/1‑99 (MQ=255)
                                                                                                               cGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAAtcttc                                    <  1:1608074/150‑1 (MQ=255)
                                                                                                                ggAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCCGTAAGAGTTGCGTTTGTCGCATAGCCTTCCAGCCAGAAATCACTCGGCTCAAtcttct                                   >  2:1483145/1‑150 (MQ=255)
                                                                                                                    tCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTAt                                                                                                               >  1:1498857/1‑70 (MQ=255)
                                                                                                                    tCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTAt                                                                                                               <  2:1498857/70‑1 (MQ=255)
                                                                                                                      aCGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGcc                             >  2:454441/1‑150 (MQ=255)
                                                                                                                        gggCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCAc                           >  2:1523603/1‑150 (MQ=255)
                                                                                                                         ggCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACt                          <  1:1009461/150‑1 (MQ=255)
                                                                                                                          gCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTc                         >  1:1681216/1‑150 (MQ=255)
                                                                                                                           cATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCa                                                                                                     <  2:1611778/73‑1 (MQ=255)
                                                                                                                           cATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCa                                                                                                     >  1:1611778/1‑73 (MQ=255)
                                                                                                                            aTTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCg                       >  2:1106013/1‑150 (MQ=255)
                                                                                                                              tGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGc                     <  2:274977/150‑1 (MQ=255)
                                                                                                                                  ccAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGAAAAAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCAt                 >  1:496765/1‑150 (MQ=255)
                                                                                                                                  ccAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCAAGCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCAt                 >  1:1573186/1‑150 (MQ=255)
                                                                                                                                       ctcTTCGCGGGTTCACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCg                                                                                                  >  2:1821859/3‑64 (MQ=37)
                                                                                                                                       ctcTTCGCGGGTTCACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCg                                                                                                  <  1:1821859/62‑1 (MQ=37)
                                                                                                                                       ctcTTCGCGGGTTCACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACtcg            >  2:43906/3‑148 (MQ=255)
                                                                                                                                             gcgGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCAGtt      >  1:1236451/1‑150 (MQ=255)
                                                                                                                                                 gTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTTGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCAGTTTCgg  <  1:549498/150‑1 (MQ=255)
                                                                                                                                                 gTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCAGTTTCgg  >  1:585042/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTGGAGCTGCGGCGGCGCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCAGTTTCGG  >  NC_002947/5866605‑5866899

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