Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,104,180 G→A 7.4% intergenic (+12/‑190) mltC → / → Y75_RS15435 membrane‑bound lytic murein transglycosylase MltC/nucleoside permease NupG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077793,104,1800GA 7.4% 205.8 / 5.1 76intergenic (+12/‑190)mltC/Y75_RS15435membrane‑bound lytic murein transglycosylase MltC/nucleoside permease NupG
Reads supporting (aligned to +/- strand):  ref base G (54/16);  new base A (4/2);  total (58/18)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.93e-01

AACCCAACATCGCGGCGTTATGCCGTCATCACCGCCTATAACGGCGGCGCAGGCAGCGTGCTGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCATCACGATGTGAGGAAATTAACATGAATCTTAAGCTGCAGCTGAAAATCCTCTCTTTTCTGCAGTTC  >  NC_007779/3103947‑3104414
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                          
aaCCCAACATCGCGGCGTTATGCCGTCATCACCGCCTATAACGGCGGCGCAAGCAGCGTGCTGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAAACGcc                                                                                                                                                                                                                                          >  1:509310/1‑236 (MQ=255)
                    tGCCGTCATCACCGCCTATAACGGCGGCGCAGGCAGCGTGCTGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTTCATTAAACGCCCCTGACGATGCCCAGGGGCaaaa                                                                                                                                                                                                                   >  2:824176/1‑239 (MQ=255)
                         tcatcaCCGCCTATAACGGCGGCGCAGGCAGCGTGCTGCGAGTCTTTTCTAATGATAAAATTCAGGCCGCCAAAATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCCGCGGAATCTCGCCGTTATCTTTATAAAAGTAATAACGCGCAAAAATTCTACCGCCCCCGATAATTCCATTAACCGCCCCGGACGATGCTCAGCGGCCAAAATGtt                                                                                                                                                                                                              >  2:487929/1‑239 (MQ=255)
                           atcaCCGCCTATAACGGCGGCGCAGGCAGCGTGCTGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTa                                                                                                                                                                                                             >  1:505892/1‑238 (MQ=255)
                               ccGCCTATAACGGCGGCGCAGGCAGCGTGCTGCGAGTCTTTCCGAATGATAAGATTAGGGCGGCAAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGCCGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCTTGACGATGCTCAGGGGCAAAAATGTTATCcac                                                                                                                                                                                                        <  1:1213328/239‑1 (MQ=255)
                                      taACGGCGGCGCAGGCAGCGTGCTGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCGTCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCATGACGATGCTCAGGGGCAAAAATGTTATCCACATCACaa                                                                                                                                                                                                 >  1:744647/1‑239 (MQ=255)
                                                         gTGCTGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGAACCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTggg                                                                                                                                                                              >  1:1459311/1‑239 (MQ=255)
                                                           gctgcGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACCCCATGACGCCGGGCGATGTTTATCAGACGCTGACGCCCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCACCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGa                                                                                                                                                                             <  1:487477/238‑1 (MQ=255)
                                                                               aTAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCcac                                                                                                                                                        >  1:1163798/1‑239 (MQ=255)
                                                                                      tCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCCCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAAc                                                                                                                                                 >  2:190625/1‑239 (MQ=255)
                                                                                         ggCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCAGTTATCTTTATAAAGTGAATACCACCCAAAACTCCTACAGCCGCCGATAATTACACTAATCGCCCCTGACGATGCCCAGGGGCAAAAATGATATCAACATCCCACTCTCTGTTTGCAAATTGGGAATGGTTGCAATTATTTGCCATAGGTATCaaa                                                                                                                                              >  2:931043/1‑239 (MQ=255)
                                                                                              ccAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAAcc                                                                                                                                         >  1:554674/1‑239 (MQ=255)
                                                                                              ccAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAAcc                                                                                                                                         >  1:48933/1‑239 (MQ=255)
                                                                                              ccAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCGGACGACCCGCCATCCGTCGGCGGAATCTCGCCGTTATCTTTAAAAAGTGAATACCGCGCAAAAATCCTACAGCAGCCGATAATTCCATTAACCACCACTGACGAGGCTCAGGGGCAAAAATGTTCTCCACATCCCAATTTCGTTTTGCAACTTGTGAATCTATGCAAGTATTTGCCATAGGTAACAACAAAcc                                                                                                                                         >  2:1166595/1‑239 (MQ=255)
                                                                                                      attaACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCACCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCcgcg                                                                                                                                 >  1:1639401/1‑239 (MQ=255)
                                                                                                        taACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGaa                                                                                                                               >  1:736624/1‑239 (MQ=255)
                                                                                                        taACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGaa                                                                                                                               >  1:603910/1‑239 (MQ=255)
                                                                                                            acCATGAGCCCGGGCGATGTTTATCAGACGCTGACCACCCGCCATCCCTCTCCGGACTCTCGCCGTTATCTCTATAAAGTCAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCTCCTGACGAAGCTCAGTCGCTCAAGTGTTATCCACATCACAATTTCGTTTTGGCAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACACTAAACCTGTCcgcgcacatg                                                                                                                           >  2:123382/1‑233 (MQ=255)
                                                                                                            acCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTAACCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAAACAGTCCGCGAAGTTg                                                                                                                           >  1:1821730/1‑239 (MQ=255)
                                                                                                              cATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAatt                                                                                                                                                                 >  1:1072417/1‑199 (MQ=255)
                                                                                                              cATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGAt                                                                                                                         >  1:1143225/1‑239 (MQ=255)
                                                                                                              cATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCACCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCcgc                                                                                                                                  >  2:226668/1‑230 (MQ=255)
                                                                                                              cATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCACCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCcgc                                                                                                                                  <  1:226668/230‑1 (MQ=255)
                                                                                                              aaggccgcgGGGCGTTGTTTATCAGCCGCTGCCACCCCCCCACCCCTCGGCGAAATCCCGCCGTTATTTTTATAAGGTGAATCCGGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGAt                                                                                                                         <  1:868311/234‑1 (MQ=255)
                                                                                                               aTGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATa                                                                                                                        >  1:759576/1‑239 (MQ=255)
                                                                                                               aTGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATa                                                                                                                        >  1:671216/1‑239 (MQ=255)
                                                                                                               aTGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTATGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATa                                                                                                                        >  2:1773370/1‑239 (MQ=255)
                                                                                                                tGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAg                                                                                                                       >  2:897812/1‑239 (MQ=255)
                                                                                                                 gACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGa                                                                                                                      >  1:14741/1‑239 (MQ=255)
                                                                                                                      cGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATcc                                                                                                                  >  2:1151420/1‑238 (MQ=255)
                                                                                                                      cGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATcc                                                                                                                  >  2:878700/1‑238 (MQ=255)
                                                                                                                        ggCGAAGTTTATCAGACGCTGACGACACGGCATCACTCCGCGGACTCTAGCCGGTATTTTTATAAAGTGAATACCGAGCACAATTCCTACGGCCGCAGATAATTCCATTAAGCGCCCATGACGATGCTCAGGGGCAAAAATGTTATCCACATCACTATTTCGTTGTGCCAATAGAGAATGTTTGCAATTATTAGCCACAGGTAAAAAAAAAACAGTCCGCTAAGGAGATAGAATccaa                                                                                                                >  2:599871/1‑236 (MQ=255)
                                                                                                                                         gCTGACGACCCGCCATCCCTCTGCGGAATCACGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCaaa                                                                                              >  2:1922267/1‑239 (MQ=255)
                                                                                                                                           tGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAAtg                                                                                            >  2:1583737/1‑239 (MQ=255)
                                                                                                                                               gcccccccATCCCTCGGCGGATTCTCGCCGTTATCTTTATAAAGTTATTCCCCGCCAAAAATCCTCCCCCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCt                                                                                        <  1:1041285/237‑1 (MQ=255)
                                                                                                                                                     ccATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCa                                                                                  >  1:715037/1‑239 (MQ=255)
                                                                                                                                                     ccACCCCCCGGCGAATTCCCGCCTTTATTTTTATAAAGTAAATCCCGCGCAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCa                                                                                  <  1:895603/239‑1 (MQ=255)
                                                                                                                                                                   aaTCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCa                                                                    >  2:569670/1‑239 (MQ=255)
                                                                                                                                                                     tctcGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCcatca                                                                                                             >  1:407861/1‑196 (MQ=255)
                                                                                                                                                                          ccGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCATCACGAt                                                             >  1:1092265/1‑239 (MQ=255)
                                                                                                                                                                                  tttATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCATCACGATGTGAGGaa                                                     >  1:1178923/1‑239 (MQ=255)
                                                                                                                                                                                    tataAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGAAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCATCACGATGTGAGGAAAt                                                   >  2:43755/1‑239 (MQ=255)
                                                                                                                                                                                     atcaaGTGAATACCGCGAAAAAATACTACCCCCGCAGATAATTACATTAACAACCCCTGACGATGGTCAGGGCCAAAAATGTTATCAACATCACAATTTCGTTTTGCAAATTGGGCATGTTTGCAATTATTTGCCACAGCGAAAACAAAAAAAGTGCGCGAAGTTGATGGAATCCCATCATCTCGCACGGTCAGATGTACTTTTTCAAACACTCATCCGGATCACGAATTTAGGAAttt                                                  >  2:1362125/4‑239 (MQ=255)
                                                                                                                                                                                     ataAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGAGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCACTTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCATCACGATGTGAGGAAAtt                                                  >  2:130331/1‑239 (MQ=255)
                                                                                                                                                                                           tGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCATCACGATGTGAGGAAATTAACATg                                            >  2:137057/1‑239 (MQ=255)
                                                                                                                                                                                            gactACCGCGCAAAAATCCTACAGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCACATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCTTCACGATGTGAGCAAATTAATATGa                                           >  2:453811/4‑239 (MQ=255)
                                                                                                                                                                                            gAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAaca                                                                              >  1:645451/1‑204 (MQ=255)
                                                                                                                                                                                            gAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCATCACGATGTGAGGAAATTAACATGa                                           >  1:35011/1‑239 (MQ=255)
                                                                                                                                                                                            aactcccccaaaaaaaCCCTCCCGCCGCCGATAATCCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCATCACGATGTGAGGAAATTAACAt                                             <  2:1459311/226‑1 (MQ=255)
                                                                                                                                                                                             aaTACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCATCACGATGTGAGGAAATTAACATGaa                                          >  1:1379900/1‑239 (MQ=255)
                                                                                                                                                                                                   gcgcAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAatt                                                                                                                                                                 >  1:1327755/1‑114 (MQ=255)
                                                                                                                                                                                                    cgcAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGa                                                                                                                                <  1:520485/146‑1 (MQ=255)
                                                                                                                                                                                                    cgcAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGa                                                                                                                                >  2:520485/1‑146 (MQ=255)
                                                                                                                                                                                                    cgcAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCATCACGATGTGAGGAAATTAACATGAATCTTaa                                    >  2:847482/1‑238 (MQ=255)
                                                                                                                                                                                                    cgcAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCATCACGATGTGAGGAAATTAACATGAATCTTaa                                    >  1:101427/1‑238 (MQ=255)
                                                                                                                                                                                                    cgcAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCATCACGATGTGAGGAAATTAACATGAATCTTaa                                    >  1:1868735/1‑238 (MQ=255)
                                                                                                                                                                                                    cgcAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCATCACGATGTGAGGAAATTAACATGAATCTTaa                                    >  2:1038322/1‑238 (MQ=255)
                                                                                                                                                                                                    cgcAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCATCACGATGTGAGGAAATTAACATGAATCTTaa                                    >  1:1578720/1‑238 (MQ=255)
                                                                                                                                                                                                    cgcAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCATCACGATGTGAGGAAATTAACATGAATCTTaa                                    >  1:1270252/1‑238 (MQ=255)
                                                                                                                                                                                                    cgcAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCATCACGATGTGAGGAAATTAACATGAATCTTaa                                    >  1:635485/1‑238 (MQ=255)
                                                                                                                                                                                                    cgcAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTAATCCGCATCACGATGTGAGGAAATTAACATGAATCTTaa                                    >  2:150183/1‑238 (MQ=255)
                                                                                                                                                                                                    cgcAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAAACAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCATCACGATGTGAGGAAATTAACATGAATCTTaa                                    >  1:610656/1‑238 (MQ=255)
                                                                                                                                                                                                    cgcAAAAATCCTACCGCCGCCGATAATTCCATAAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAAt                                                                                             >  1:285699/1‑181 (MQ=255)
                                                                                                                                                                                                                cccgccgccgATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCATCACGATGTGAGGAAATTAACATGAATCTTAAGCTGCAGCTGaaa                       <  1:897812/238‑1 (MQ=255)
                                                                                                                                                                                                                 ccgccgccgATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGc                                                                                                       >  2:1458297/1‑158 (MQ=255)
                                                                                                                                                                                                                 ccgccgccgATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGc                                                                                                       <  1:1458297/158‑1 (MQ=255)
                                                                                                                                                                                                                    ccgccaATATTCCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCATCACGATGTGAGGAAATTAACATGAATCTTAAGCTGCAGCTGAAAATcc                   <  1:802355/239‑1 (MQ=255)
                                                                                                                                                                                                                    acgccgATAATTCCATTAACCGCCCCTGACGATGCTCAGGAGCAAAAAAGCTATACACATCACAATTTCGTTTTCCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACAAGTCCGAGAAGTTGATAGAATACCATCATCTCGCACGATCAAATGTGCTTTTTCAAACACTCATCCGCATCACGCTGTGAGGAAATGACCATGACTCTTCAGTTGCAGCTGAAAAt                     >  2:416210/2‑237 (MQ=255)
                                                                                                                                                                                                                     cgcctaaatttCAATTACCCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCATCACGATGTGAGGAAATTAACATGAATCTTAAGCTGCAGCTGAAAATCct                  <  2:603910/230‑1 (MQ=255)
                                                                                                                                                                                                                          aacaccccATCAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTATTGCAAATTGTGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACTCTCATCCGCATCACGATGTGAGGAAATTAACATGAATCTTAAGCTGCAGCTGAAAATCCTCTCtt             <  2:434826/233‑1 (MQ=255)
                                                                                                                                                                                                                          aTAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCATCACGATGTGAGGAAATTAACATGAATCTTAAGCTGCAGCTGAAAATCCTCTCtt             <  1:1600809/239‑1 (MQ=255)
                                                                                                                                                                                                                            aaTCCCAGTAGCCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCATCACGATGTGAGGAAATTAACATGAATCTTAAGCTGCAGCTGAAAATCCTCGCtgtt           <  2:554674/239‑1 (MQ=255)
                                                                                                                                                                                                                              ttCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCATCACGATGTGAGGAAATTAACATGAATCTTAAGCTGCAGCTGAAAATCCTCTCTTTTCt         <  1:847482/239‑1 (MQ=255)
                                                                                                                                                                                                                               tCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCATCACGATGTGAGGAAATTAACATGAATCTTAAGCTGCAGCTGAAAATCCTCTCTTTTCTg        <  1:1151420/239‑1 (MQ=255)
                                                                                                                                                                                                                                     aaCCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCATCACGATGTGAGGAAATTAACATGAATCTTAAGCTGCAGCTGAAAATCCTCTCTTTTCTGCAGTtc  <  2:759576/239‑1 (MQ=255)
                                                                                                                                                                                                                                      aCCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCATCACGATGTGAGGAAATTAACATGAATCTTAAGCTGCAGCTGAAAATCCTCTCTTTTCTGCAGt    >  1:144514/1‑236 (MQ=255)
                                                                                                                                                                                                                                         gCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGCATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCATCACGATGTGAGGAAATTAACATGAATCTTAAGCTGCAGCTGAAAATCCTCTCTTTTCTGCAGt    >  1:77044/1‑233 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                          
AACCCAACATCGCGGCGTTATGCCGTCATCACCGCCTATAACGGCGGCGCAGGCAGCGTGCTGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCATCATCTCGCACGGTCAAATGTGCTTTTTCAAACACTCATCCGCATCACGATGTGAGGAAATTAACATGAATCTTAAGCTGCAGCTGAAAATCCTCTCTTTTCTGCAGTTC  >  NC_007779/3103947‑3104414

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