Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,832,077 G→A 5.8% V541M (GTG→ATG)  yicH → putative inner membrane‑anchored periplasmic AsmA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,832,0770GA 5.8% 193.6 / 4.2 72V541M (GTG→ATG) yicHputative inner membrane‑anchored periplasmic AsmA family protein
Reads supporting (aligned to +/- strand):  ref base G (47/20);  new base A (2/2);  total (50/22)
Fisher's exact test for biased strand distribution p-value = 5.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.21e-01

CGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCA  >  NC_000913/3831839‑3832303
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                    
cGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGCTATTCAGGCCAATGTCCCGTTGAAACCTCCGGTTCCCGGGCAATTGCATGCCg                                                                                                                                                                                                                                      >  1:214312/1‑239 (MQ=255)
       aGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGCATTACCGTTGACTGGCGACGCCAATATTCAGCTTACCGCCAGGGGCGATATTTAGGCCAATGGCCGGTGTAAGACCACGGGTAGCGGGCAAAGGCACGCCCGGAAcgc                                                                                                                                                                                                                               >  2:339618/1‑239 (MQ=255)
        gTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAAcgcc                                                                                                                                                                                                                              >  1:279348/1‑239 (MQ=255)
                aCTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTTAAACCTACGGTTACCGGGCAATTGCATGGCGAGAACGCCGCCAagc                                                                                                                                                                                                                       >  2:469714/1‑238 (MQ=255)
                   gAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATTGCCTGAATTACAGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAAGCCAAGGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGAATGACGTGAAGCCCGAAAAgcggca                                                                                                                                                                                                                   >  1:700572/1‑239 (MQ=255)
                   gAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCACCAAAgaagca                                                                                                                                                                                                                   >  1:586596/1‑239 (MQ=255)
                   gAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATCACGTGAACGCCGCAAAgcagca                                                                                                                                                                                                                   >  1:330412/1‑239 (MQ=255)
                              tCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGCTATTCAGGCCAATGTCCCTTTGAAACCTACGGTTAGCGGGCAAATGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCaa                                                                                                                                                                                                         >  1:491003/1‑238 (MQ=255)
                              tCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGGGAACGCCGCAAACCAGCAAGTGATTCaa                                                                                                                                                                                                         >  1:506475/1‑238 (MQ=255)
                              tCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGCGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGGTATTCCGGCCAATGTCCCGTTGAAACCTACGGTTAAGGGGCCATTGCATGCCGGAAACTCGCCAAAGAAGAAAGAAACTCaaa                                                                                                                                                                                                        >  2:272398/1‑239 (MQ=255)
                                     gCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGaa                                                                                                                                                                                                 >  1:878420/1‑239 (MQ=255)
                                             gattttttCAGAAACGCCCAAACGTCAGACACATAACAGGTGGAATGGACGGGGTGTGCGGGTAAATATTTTGAAACAGGGGGTATGGCCTGAATTCCCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCa                                                                                                                                                                                        <  1:182672/234‑1 (MQ=255)
                                                    tCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTCCCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACACCGCAAAACAGCAAGGGACTACAAACATAAAATCCGGCATCGttt                                                                                                                                                                                  >  2:637098/1‑238 (MQ=255)
                                                                aaagtcaaGACAATTATCGCCCTGAAGGGCCGCGGGGGCCCGGTGAATTTTTGGCAACAGTGGGGTTGGCCTGAATTACCGTGGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGt                                                                                                                                                                     <  1:570954/232‑1 (MQ=255)
                                                                              tatCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGCAGTTACATCGCCGGAg                                                                                                                                                        >  1:762815/1‑238 (MQ=255)
                                                                               atCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATCCTGGCATCGTTTCCACCGGGGAAGTTACATCGACGGAGcc                                                                                                                                                     >  1:59374/1‑239 (MQ=255)
                                                                               atCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGAGACTCAAACCATGAATGCCGGCATCGTTTCCAGCGGTGAAGTAACATCAAAGGAGcc                                                                                                                                                     >  2:103282/1‑239 (MQ=255)
                                                                                    ccTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACa                                                                                                                                                                 >  1:790291/1‑223 (MQ=255)
                                                                                    ccTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACCCAAACCATGAATGCTGGCATCGTTTCCAGCGGGGAAGTTAAATCCACGGAACCGGTgg                                                                                                                                                >  2:714123/1‑238 (MQ=255)
                                                                                       gAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGATTCAAACCATGAATGCTGGCAACGTTTCCAGCGGTGAAGTTACATCGACGGAACCGGTGCggg                                                                                                                                             >  1:581675/1‑238 (MQ=255)
                                                                                        aaTGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCACTTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACACTGAATCTTGGCATCGTTACCCGCGGTGAAGTTAAAACGACGCGGCCGGGGCGGtg                                                                                                                                            >  2:963971/1‑238 (MQ=255)
                                                                                        aaTGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTa                                                                                                                                            >  1:300239/1‑239 (MQ=255)
                                                                                        aaTGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGACCTGAATTACCGTTGACTGGCGCCGGCAATATTCAGATTAACGCCAGTGACGCTATTCAGGCCAATGTCCCGCTGAACCCTACGGTTAGCGAGCAATTGCATGCAGTGAACGCCCCACACCAGCAGGTGACTCAAACCATGCATCTGGCCATCGGTTCCACCGGTAAAGTAACTCCAACGAAGCCGATACggaa                                                                                                                                            >  2:847616/1‑237 (MQ=255)
                                                                                              cgcgGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCCAGGGAACCGGTGCGGAAACaaac                                                                                                                                     >  2:1005975/1‑236 (MQ=255)
                                                                                              cgcgGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGAGACCCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCACCGAGACCGGTGCGGGAAAAAta                                                                                                                                      >  1:506870/1‑238 (MQ=255)
                                                                                                             aaTATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAg                                                                                                                                                                              >  1:877445/1‑185 (MQ=255)
                                                                                                             aaTATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCATGAACGCCGCAAAGCAGCAAGGGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTAAATCGAAGGAGCCGGTGCGGGAAAAATCAATCCCGCCggggag                                                                                                                       >  2:1005855/1‑236 (MQ=255)
                                                                                                               tatTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCCTCGTTTCCAGCGGTGAAGTTCCATCGACGGAGCCGGTGCGGGAAACATTAATCTCCCCGGGTACg                                                                                                                      >  1:621294/1‑238 (MQ=255)
                                                                                                                    tGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCAAGCGGTGATGTGACATTGCCGGTGCCGGTGCGGTAAACTTCAATTGCGCCGGGTACGACTCAg                                                                                                                >  1:295265/1‑239 (MQ=255)
                                                                                                                            tGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCGGACGAAccc                                                                                                        >  1:236207/1‑239 (MQ=255)
                                                                                                                            tGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATACTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGGAAACATCGATCTCGCCGGGTACGACTCAGACGTAccc                                                                                                        >  1:909633/1‑239 (MQ=255)
                                                                                                                                  tGGCCTGAATTTCCTTTGACGGGCCACGGCAATATCCAGTTTCCCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCAt                                                                                                   <  1:502055/238‑1 (MQ=255)
                                                                                                                                    gCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACa                                                                                                                                        >  1:694045/1‑199 (MQ=255)
                                                                                                                                      cTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCACCGGTGAAGTTACATCGACGGGGCCGGTGCGGTAAACATCAATCTCGCCCGGGACGACTCAGACGTACCCGGACTTCCAt                                                                                              >  1:703113/1‑239 (MQ=255)
                                                                                                                                         aaTTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACAt                                                                                                                                                                >  1:169991/1‑171 (MQ=255)
                                                                                                                                         aaTCACCGTGGCCGGGGGACCGAAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGACCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCa                                                                                            <  1:60493/238‑1 (MQ=255)
                                                                                                                                              ccGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCaa                                                                                                                                                                                                         <  1:684485/126‑1 (MQ=255)
                                                                                                                                              ccGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCaa                                                                                                                                                                                                         >  2:684485/1‑126 (MQ=255)
                                                                                                                                               cGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATa                                                                                         >  1:240025/1‑235 (MQ=255)
                                                                                                                                                gtgagcgggCGAAGGAAATTTTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTa                                                                                    <  1:272398/232‑1 (MQ=255)
                                                                                                                                                   gACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCTTTCCTCCAATAAATTCCaa                                                                                 >  1:496975/1‑239 (MQ=255)
                                                                                                                                                      tGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCCGGAACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGt                                                                              >  1:96535/1‑239 (MQ=255)
                                                                                                                                                                  tatTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAACCATCAATATCGCCGGGTAGAAAGCAGACGAACACGCCCTTCCAACAATAGGTTTAAACGTAAATGTCAgcga                                                                  >  2:50796/1‑238 (MQ=255)
                                                                                                                                                                   tttactgttcccGTCAGTGGCGATAGTCAGGCCAAGGTCCCGTTGAAACCTGCGGTTAGCGGGCAATTGCATGCCGTGAAAGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGGTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCa                                                                 <  1:457958/229‑1 (MQ=255)
                                                                                                                                                                    ttcgcttcgccgccgGGGGTGATTTTCAGGCCCATGTCCCGGTGAAACCTCCGGTTACGGGGCAATTGGCTGCCGTCACCGCCGCAAAGCAGCAAGTGACCCAACCCATGAATGTTGCCATCGTTTCCAGCGGTGAAGTTACATCGACGGGGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCAt                                                                <  2:506870/230‑1 (MQ=255)
                                                                                                                                                                    tgtaggttcccGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCAt                                                                <  1:963971/230‑1 (MQ=255)
                                                                                                                                                                         cTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGATGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCt                                                           >  2:1022360/1‑239 (MQ=255)
                                                                                                                                                                            aCCGCCAGTGGCGATTTTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCAt                                                                <  2:67334/231‑1 (MQ=255)
                                                                                                                                                                            aCCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCAt                                                                >  1:67334/1‑231 (MQ=255)
                                                                                                                                                                                 cAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCtt                                                   <  1:103282/239‑1 (MQ=255)
                                                                                                                                                                                 cAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCtt                                                   <  1:714123/239‑1 (MQ=255)
                                                                                                                                                                                   gTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGCAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCCCCGGGTACGACCCAGACGTACCCGGCTTTCCTACAATAGATTACACCGAACTTTTCACGGACTTCCCTTTAGGCTTc                                                  >  2:743377/1‑238 (MQ=255)
                                                                                                                                                                                   gTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCCGACGTACCCGGCATTCCATCAATTGATTACAACGTAATTGGCCGCGCATTCCCTTGTGGCTTc                                                  >  1:963647/1‑238 (MQ=255)
                                                                                                                                                                                   gTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCATGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTca                                                 <  1:1005855/239‑1 (MQ=255)
                                                                                                                                                                                       cGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTcacca                                              <  1:1007510/238‑1 (MQ=255)
                                                                                                                                                                                       cGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCATGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCg                                                                                                                                                    <  2:749068/136‑1 (MQ=255)
                                                                                                                                                                                       cGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCATGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCg                                                                                                                                                    >  1:749068/1‑136 (MQ=255)
                                                                                                                                                                                        gATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTTAGGCTTCaccacc                                            >  1:288200/1‑239 (MQ=255)
                                                                                                                                                                                             tccggCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATGTCATGCCGTGAACGCCGCAAAGCAGCAAGTTACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGccc                                       <  2:470588/236‑1 (MQ=255)
                                                                                                                                                                                             gtgggCCATTGTCCGGGTGAAACCGACCGTTAGCGGGCAATTGCATGCCGGTAACGCCGAAAAGCGGCAAGGGACTAAAACCGGGAATGTGTGGAGAGGTTACAGCGGTGAAGGTACATCGACGGGGCCGGTTCGGTAAACATCAATCTCGCCGGTTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGccc                                       <  2:762815/236‑1 (MQ=255)
                                                                                                                                                                                                  cataTGTCCCTTTGAACCCTTCGGTTAGGGGGCAATTGGATGCCGGTGACGCCGGAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTTAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCt                                  <  2:59374/236‑1 (MQ=255)
                                                                                                                                                                                                      tGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTACTTGTCACCGCCTTCCCTTGAGGCTTCACAACCACCCCCTACCCACt                              >  2:26719/1‑239 (MQ=255)
                                                                                                                                                                                                      tGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAAAATCAATCTCGCCGGGTAGGACCCAGACGTACCCGGCATTCCATAACTAGATTACAACGAAATTGTCAGCGCATTCCCCTGAGGCTTCACCACAAACCCCGGccacat                              >  2:860007/1‑234 (MQ=255)
                                                                                                                                                                                                         cccGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCTGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTc                           <  1:50796/239‑1 (MQ=255)
                                                                                                                                                                                                           cGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGATGCGGTAAACATCAATCTCCCCGGGTACGACTCACGCGTACCCGGCATTCCATAAATAGATTAAAACGTAATTGCCAGCGCATCACCTAGAGGCTTCACCACCAACCCCTCCCCACGTTCa                          >  2:472055/1‑238 (MQ=255)
                                                                                                                                                                                                             ttGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGc                        >  1:680013/1‑238 (MQ=255)
                                                                                                                                                                                                                aaaCCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGGGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCtcc                    >  1:516591/1‑237 (MQ=255)
                                                                                                                                                                                                                   ccTTCGGTTAGCGGGAAATTGCATGGCGTGAACCCCGGAAAGGAGCAAGTGACTCAAACCATGAATGCTTGCATCGTTTCCAGCGGTGAATTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCAGCAATAGATTACAACGTAATTGTCAGCGTATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCg                   <  1:305477/237‑1 (MQ=255)
                                                                                                                                                                                                                   ccTACGGTTAGGGGGCAATTGCATGCCGTGACCGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGaa                 <  1:469714/239‑1 (MQ=255)
                                                                                                                                                                                                                   ccTACGGTTAGCGGGAAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGaa                 <  1:85775/239‑1 (MQ=255)
                                                                                                                                                                                                                            aGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCCCTTTCAGCCTGCGAACCCGCTTGc        >  1:864627/1‑239 (MQ=255)
                                                                                                                                                                                                                            aGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCgctgtct        >  1:697791/1‑233 (MQ=255)
                                                                                                                                                                                                                                  cAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTCCAACGTAATTTTCAGCGCATTCCCTTAGGCCTTCACCCCCAGACCCAGCTCACTTTCAACATGCGAACAGACTTGCAAACCa  >  2:90575/1‑239 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                    
CGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCA  >  NC_000913/3831839‑3832303

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