Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,169,668 A→C 14.5% intergenic (+57/+25) yegS → / ← gatR phosphatidylglycerol kinase, metal‑dependent/pseudogene, repressor for gat operon, interrupted by IS3, split galactitol utilization operon repressor, fragment 2, split galactitol utilization operon repressor, interrupted

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,169,6680AC14.5% 136.6 / 2.3 59intergenic (+57/+25)yegS/gatRphosphatidylglycerol kinase, metal‑dependent/pseudogene, repressor for gat operon, interrupted by IS3, split galactitol utilization operon repressor, fragment 2, split galactitol utilization operon repressor, interrupted
Reads supporting (aligned to +/- strand):  ref base A (24/24);  new base C (2/9);  total (26/33)
Fisher's exact test for biased strand distribution p-value = 9.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCCGAATATTATCGAAGGCGCTTCGTCGTGGTTTGATATTCAGGCACCACACGACATCACCTTTAATCTTGATGGCGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTAAATTTCGATGAATCGAGAACGGCAATCACTTTATTAGCGGCTGTGACCATCGCACCGCTAATAGAATAACCTTCATTAAAGGT  >  NC_000913/2169444‑2169902
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                          
aaCCCGAATTTTATCGAAGGGGTTTCGTCGTGGTTGGATATTCAGGCACCACAGGACATGACCTTTAATCTTGATGGCGAACCGTTGAGTGGGCAAAATTTTCATAGTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATa                                                                                                                                                                                                                              <  1:226117/239‑1 (MQ=255)
   ccGAATATTATCGAAGGCGCTTCGTCGTGGTTTGATATTCAGGCACCACACGACATCACCTTTAATCTTGATGGCGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAg                                                                                                                                                                                                                           <  2:31867/239‑1 (MQ=255)
   ccGAATAATATCGAAGGCGCTTCGTCGTGGTTTGATATTCAGGCACCACACGACATCACCTTTAATCTTGATGGCGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAg                                                                                                                                                                                                                           <  2:672450/239‑1 (MQ=255)
   ccGAAATTTTTCGAAGGCGCTTGTGTGTGGTTTGATATTCAGGCACCACACGACATCACCTTTAATCTTGATGGCGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAg                                                                                                                                                                                                                           <  1:258880/239‑1 (MQ=255)
     ggattttctccaaGGCCTTTTGTCGTTGTTTTATTTTCAGGCACCACACGGCCTCACCTTTAATCTTGATGGCGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGaa                                                                                                                                                                                                                         <  1:274741/228‑1 (MQ=255)
     gAATATTATCGAAGTCGCTTCGTCGTGGTTTGATATTCAGGCACCACACGACATCACCTTTAATCTTGATGGCGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGaa                                                                                                                                                                                                                         <  2:103504/239‑1 (MQ=255)
     gAATATTATCGAAGGCGCTTCGTCGTGGTTTGATATTCAGGCACCACACGACATCACCTTTAATCTTGATGGCGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGaa                                                                                                                                                                                                                         <  2:433751/239‑1 (MQ=255)
     gAATATTACCGAAGGCGCTTCGTCGTGGTTTGATATTCAGGCACCACACGACATCACCTTTAATCTTGATGGCGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGaa                                                                                                                                                                                                                         <  2:647335/239‑1 (MQ=255)
                     gtttCGTCGGTGTTTTATATTCGGGCACCACCCGACCTCACCTTTAATCTTGATGGCAAAACGGTGAGTGGGCAAAATTTTCATATTGAAATTCGTCCGGAAGCGGTGCGTTGTCGATTACCACCAGATTGGCCACTGTTGCGTTAATCAAACTATTCTCTCATGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGtaa                                                                                                                                                                                                         <  1:70701/237‑1 (MQ=255)
                                  ttGATATTCAGGCCCCACACGACATCACCTTTAATCTTGATGGCGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTAccc                                                                                                                                                                                             <  2:594536/238‑1 (MQ=255)
                                     atatTCAGGCACCACACGACATCACCTTTAATCTTGATGGCGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACCCCGTAATTAATTTTACCcgcg                                                                                                                                                                                          >  2:655667/1‑238 (MQ=255)
                                              cccctccccccATCCCCTTTAATCTTGTGGGCGAACCGTTGAGTGGGCAAAAATTTCATATTGAAATCCTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTTTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTg                                                                                                                                                                                   <  1:866015/226‑1 (MQ=255)
                                               accacACGACATCACCTTTAATCTTGATGGCGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGGGTCTTCTGCaa                                                                                                                                                                                >  2:813467/1‑238 (MQ=255)
                                                   cacGACATCACCTTTAATCTTGATGGCGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACCCCTTATCCGGCATACATACCAATACTACACCGTAATTAATTTTACCCACGCCTTCGGCAATGCCa                                                                                                                                                                           >  1:17072/1‑239 (MQ=255)
                                                                ttAATCTTGATGGCGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACCCCTTATCCGGCATACAGAACACAACTACAACGAACTTAATTTTACCCGCGTCTTCTGCAATGCCCGTTTATCCAcc                                                                                                                                                                >  2:898290/1‑237 (MQ=255)
                                                                ttAAACTTGAGGGCGAACCGTTGAGTGGGCAAAATTTTCATATTTAAATACTTCCGGACGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACctc                                                                                                                                                              <  2:227083/239‑1 (MQ=255)
                                                                   attttgttgGCGACCCGTTGAGTGGGCAAAATTTTCATATTGAAATTCTTCCGGCAGCGTTGCGTTGCCGATTACCCCCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGAGGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGa                                                                                                                                                           <  1:813467/236‑1 (MQ=255)
                                                                        gATGGCGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGAcacac                                                                                                                                                      >  2:304048/1‑238 (MQ=255)
                                                                               aaCCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCCCCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTTCCCCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGcatcat                                                                                                                                                <  2:16817/238‑1 (MQ=255)
                                                                               aaCCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGcatcat                                                                                                                                                <  2:17072/238‑1 (MQ=255)
                                                                                         tGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCAACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATTCCGGATGCGGTTTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGAGTCTTCTGCAATGACAGTTTATCACCCCCTGAACAACCATCATCGGtaaaaaaa                                                                                                                                     >  1:368604/1‑237 (MQ=255)
                                                                                                                     ccGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATACCAGTTTATCCACCTTTTACACAGCATCATCGGTAATAATAATGTCAATTTTTTATATTGgtagta                                                                                                            >  1:722430/1‑236 (MQ=255)
                                                                                                                       ggCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGaa                                                                                                                                                                                                                         <  1:42286/125‑1 (MQ=255)
                                                                                                                       ggCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCAGCATACAGaa                                                                                                                                                                                                                         >  2:42286/1‑125 (MQ=255)
                                                                                                                         cctgcttgctttgGCGATTCCCACCGAATTGCCCCTTTTTGCGTTAATCAACCTATTCTCTCTTGGACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAACCACCTTATCCGGCATACAGAACAATCCTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGAtt                                                                                                     <  1:630001/234‑1 (MQ=255)
                                                                                                                            cgttgcgttgTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCAACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGAACCTGATTaaa                                                                                                  >  1:676330/1‑239 (MQ=255)
                                                                                                                                 cgttgTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGCCTTCTGCCATGCCAGTTTATCCACCTCTAACACAGCATCCCCGGTCATAATAATGTCAATTTTTTCGTTTGGAATAACCGAACTAAAAcccc                                                                                             >  1:298744/1‑237 (MQ=255)
                                                                                                                                      tCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCCCGTCTTCTGCAATGCCAGTTTATCCACCCCTGACACAGCACCATCGGTAATAATAATTTCAATTTTTTCGATTGTAAGAACCGAATTAAAACCGCGACgg                                                                                        >  1:431036/1‑239 (MQ=255)
                                                                                                                                           tACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGa                                                                                                                                                           <  2:562638/167‑1 (MQ=255)
                                                                                                                                           tACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGa                                                                                                                                                           >  1:562638/1‑167 (MQ=255)
                                                                                                                                                        tCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGCAACAGCATCATCGGTATTAATAATGTCAATTTTTTCGATTGGTAGACCCTGTTTAAAACCGGGACGGTAAAATTTCGATAAATCg                                                                      >  1:260831/1‑239 (MQ=255)
                                                                                                                                                        tCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTAAATTTCGATGAATCg                                                                      >  2:546342/1‑239 (MQ=255)
                                                                                                                                                        tCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTAAATTTCGATGAATCg                                                                      >  2:750441/1‑239 (MQ=255)
                                                                                                                                                            cttttGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTAAATTTCGATGAATCGAGaa                                                                  <  2:893283/236‑1 (MQ=255)
                                                                                                                                                                 ggtttAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTTTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTAAATTTCGATGAATCGAGAACGGCa                                                             <  1:546342/236‑1 (MQ=255)
                                                                                                                                                                  cGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTAAATTTCGATGAATCGAGAACGGCaa                                                            >  2:268554/1‑239 (MQ=255)
                                                                                                                                                                  cGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTAAATTTCGATGAATCGAGAACGGCaa                                                            >  2:226353/1‑239 (MQ=255)
                                                                                                                                                                  cGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCCCGACGGTAAAATTTCGATGAATCGAAAACGGCaa                                                            >  1:415216/1‑239 (MQ=255)
                                                                                                                                                                                               gggAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTAAATTTCGATGAATCGCGAACGGCAATCACTTTATTAGAGGCTGTAAcacacca                                >  1:349930/1‑232 (MQ=255)
                                                                                                                                                                                               gggAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTAAATTTCGATGAATCGAGAACGGCAATCACTTTATTAGCGGCTGTGACCATCGc                                >  2:655413/1‑238 (MQ=255)
                                                                                                                                                                                               gggAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTAAAATTTCGATGAATCGAGAACGGCAATCACTTTATTAGAGGcgggaacccacgac                               >  1:23707/1‑226 (MQ=255)
                                                                                                                                                                                                 ggaaacccGTAAGGCCGGGTGCGGTGTAAACCCCTTACCCGGCATAAAGAAAAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTAAATTTCGATGAATCGAGAACGGCAATCACTTTATTAGCGGCTGTGACCATCGCAc                              <  1:750441/232‑1 (MQ=255)
                                                                                                                                                                                                 gagaATCCCTAATGCCCGATGCGGTTGAAACCCCTTATCCGGCATACAGAACAATACTACACCGGAATTACATTTACCCGCGTCTTCTGCAATGCCAGTTTAACCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTAAATTTCGATGAATCGAGAACGGCAATCACTTTATTAGCGGCTGTGACCATCGc                                <  2:492986/236‑1 (MQ=255)
                                                                                                                                                                                                 gaaaacccGTAATGCCGGATGGGGTTTAAACCCCTTACCCGGCATACAGAACAATTCTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTAAATTTCGATGAATCGAGAACGGCAATCACTTTATTAGCGGCTGTGACCATCGc                                <  1:226353/230‑1 (MQ=255)
                                                                                                                                                                                                  aaaaTCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTAAATTTCGATGAATCGAGAACGGCAATCACTTTATTAGCGGCTGTGACCATCGCAcc                             <  2:121147/236‑1 (MQ=255)
                                                                                                                                                                                                    aaTCCGTAATGCCGGATGCGGTGTAAACCCCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTAAATTTCGATGAATCGAGAACGGCAATCACTTTATTAGCGGCTGTGACCATCGCACCGc                           <  2:431036/238‑1 (MQ=255)
                                                                                                                                                                                                    aaTCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGGGACGGTAAAATTTCGCTGAATCGAGAAGGGCAATAACTTTTTTACCGGCTGGTACCATCCCACCGca                          >  1:807424/1‑238 (MQ=255)
                                                                                                                                                                                                     aTCCGTAATGCCGGATGCGGTGTAAACCCCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTAAATTTCGATGAATCGAGAACGGCAATCACTTTATTAGCGGCTGTGACCATCGCACCGCt                          <  2:676330/238‑1 (MQ=255)
                                                                                                                                                                                                     aTCCCTTATTCCGGATGGGGTTTAAACCCCTTACCCGGAATAAAGAAAAATACTACCCCGGAATTAATTTTACCCGTGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTAAATTTCGATGAATCGAGAACGGCAATCACTTTATTAGCGGCTGTGACCATCGCACCGCt                          <  1:448622/238‑1 (MQ=255)
                                                                                                                                                                                                              gCCGGCTGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTAAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTAAATTTCGCTGAATCGTGAACGAGAATAACTTTATTAGCGGCGGTAACCATccccccccaaaaagataaa                >  1:394509/1‑221 (MQ=255)
                                                                                                                                                                                                               ccGGATGCGGTGTAAACCCCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTCCCCGCGTCTTCTGCAATGCCAGTTTATCCCCCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTAAATTTCGATGAATCGAGAACGGCAATCACTTTATTAGCGGCTGTGACCATCGCACCGCTAATAGAATAAc               <  2:777286/239‑1 (MQ=255)
                                                                                                                                                                                                                cGGATGCGGTGTAAACCCCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTAAATTTCGATGAATCGAGAACGGCAATCACTTTATTAGCGGCTGTGACCATCGCACCGCTAATAGAATAAc               <  2:802963/238‑1 (MQ=255)
                                                                                                                                                                                                                cGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATCCTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTAAATTTCGATGAATCGAGAACGGCAATCACTTTATTAGCGGCTGTGACCATCGCACCGCTAATAGAATAAc               <  2:324974/238‑1 (MQ=255)
                                                                                                                                                                                                                cGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTAAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTAAATTTCGATGAATCGAGAACGGTAATCACTTTATTAGCGCCTATGCCCATCGAACCGCTACTAtaataa                >  2:675756/1‑237 (MQ=255)
                                                                                                                                                                                                                  gATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATTCCAGTTTATCCACCTCTGACACAGCCTCATCGGCAATAATAATGTCAATTTTTTCGATTGGTAGACCCTGATTAAAACCGCGACTGTTAAATTTCGATAAATCGAAAACGTAAATCACTTTCGTAACGGCTTTGACCACCGCACCGCTaaaacataaaccat            >  1:177800/1‑227 (MQ=255)
                                                                                                                                                                                                                  gATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTAc                                                                                                           <  1:304634/144‑1 (MQ=255)
                                                                                                                                                                                                                  gATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTAc                                                                                                           >  2:304634/1‑144 (MQ=255)
                                                                                                                                                                                                                          tttaaCACCTTATCCGGCAAACAAAACAATTCTACACCGTAATTAATTTTACCCGCGTTTTCTGCAATGCCAGTTTAACCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTAAATTTCGATGAATCGAGAACGGCAATCACTTTATTAGCGGCTGTGACCATCGCACCGCTAATAGAATAACCTTCATTaaa     <  2:781870/235‑1 (MQ=255)
                                                                                                                                                                                                                            aaaCCCCTTATCCGGCATACAGAACAATCCTACACCGTAATTAATTTTCCCCGCGGCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTAAATTTCGATGAATCGAGAACGGCAATCACTTTATTAGCGGCTGTGACCATCGCGCCGCTAATAGAATAACCTTCATTAAAGGt  <  2:298744/239‑1 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                          
AACCCGAATATTATCGAAGGCGCTTCGTCGTGGTTTGATATTCAGGCACCACACGACATCACCTTTAATCTTGATGGCGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTAAATTTCGATGAATCGAGAACGGCAATCACTTTATTAGCGGCTGTGACCATCGCACCGCTAATAGAATAACCTTCATTAAAGGT  >  NC_000913/2169444‑2169902

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