Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,169,668 A→C 15.6% 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,6680AC15.6% 125.6 / 2.6 56intergenic (+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 (21/24);  new base C (2/9);  total (23/33)
Fisher's exact test for biased strand distribution p-value = 1.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

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