Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,638,227 A→G 22.0% E196G (GAG→GGG)  pitA → phosphate transporter, low‑affinity, tellurite importer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,638,2270AG22.0% 81.8 / 3.0 56E196G (GAG→GGG) pitAphosphate transporter, low‑affinity, tellurite importer
Reads supporting (aligned to +/- strand):  ref base A (20/21);  new base G (5/9);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 5.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTC‑GGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACG‑CGCCC‑GTATCCACCTGACCCCAGCGGAGCG‑TG‑AAAAGAAAGACGGCAAGAAAAA‑GCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGCGTTTTCGCACGGCGCGAAC‑GATGGTCAG‑AAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGT  >  NC_000913/3638013‑3638449
                                                                                                                                                                                                                         |                                                                                                                                                                                                                                     
aTTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTC‑GGCCTGGTGTTTGCTGGCGGTCTGATTTTTTTTCCCCGTTCCCACTTGATCGGCACCAAGAAAAG‑CGCCC‑CTACCCCCCTGACCCCCACGGGGCG‑GG‑Aaaaggaa                                                                                                                                                                                                                        >  1:463877/1‑228 (MQ=255)
  tGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCGATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTC‑GGCCTGGGGTTTGCTGGCGGTCTGATTTTCTTTCTGCGTCGCTCCTGGGGGGGGCACAAGAAACG‑CGCCC‑GTATCCACCCGACCCCCCCGGGGGG‑GG‑AAaaaaaaga                                                                                                                                                                                                                      >  1:619236/1‑228 (MQ=255)
                        taaacaaagggtgtgatgaCGGGGCGGGTAGTGGGGGATCCCCTCAAATCCCCTGAAGTATTAAGTTTTTCCGGTTCTTGGATCTTTTCCCCGTGTGGC‑GGCCGGGGGTTTGCTGGGGTGCGAATTTTCTTGCTGCGTCGCTACGGGGGCGGCACCAAGAAGCG‑CGCCC‑GTACCCACCGGACCCAAGGGGGGCG‑TG‑AAAAGAAAGACGGCAAGAAAAA‑GCCGCCGTTc                                                                                                                                                                                               <  1:397977/215‑1 (MQ=255)
                              aatggtttaatggCCGGGACGGTGGGGGTGGATGCCCTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTC‑GGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACG‑CGCCC‑GTATCCACCTGACCCCAGCGGAGCG‑TG‑AAAAGAAAGACGGCAAGAAAAA‑GCCGCCGTTCTGGAc                                                                                                                                                                                          <  2:338532/218‑1 (MQ=255)
                              aagggtgtgatgaCCGGGACGGCCGGGGTGGAAGCACTCAATTTCCCGAAAGTATTAAGTATTTCCGGTTCCCTGATCGTTTCCCCTATTGTC‑GGCCGGTTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACG‑CGCCC‑GTATCCACCTGACCCCAGCGGAGCG‑TG‑AAAAGAAAGACGGCAAGAAAAA‑GCCGCCGTTCTGGAc                                                                                                                                                                                          <  2:110958/220‑1 (MQ=255)
                                tGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGGAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTC‑GGCCTGGTGTTTGCTGGCGGTCTGATTTTTTTGCTTCGTCGCTACTGGGGCGGCCCCAAGAAACG‑GGCCC‑CTATTCACCTGACCCCACCGGAGCG‑TG‑CAAAGAAAGACGGCAAGAAAAA‑ACCGCCCTTCtcggccc                                                                                                                                                                                        >  2:34107/1‑221 (MQ=255)
                                tGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTC‑GGCCTGGTGTTTGCTGGCGGTCTGGTTTTCTTGCTGCCTCGCTCCTTGCGCGGCACCAAGAAACT‑CGCCC‑CGTTCCCCCCGCCCCCAACGGAGCG‑GC‑GAAAAAAAAACGGGAAGAAAAA‑accacccctccgggccgg                                                                                                                                                                                       >  1:601671/1‑216 (MQ=255)
                                tGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTC‑GGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACG‑CGCCC‑GTATCCACCTGACCCCCGCGGGGGG‑TG‑AAAAGAAAGACGGCAAGACAAA‑GCCTCCCtt                                                                                                                                                                                                >  2:397977/1‑219 (MQ=255)
                                     tgatgaCCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTC‑GGCCTGGTGTTTGCTGGCGGTCTGATTTTTTTTCTGCGTCGCTACTTGGACGGCACCAAAAAAAG‑CGCCC‑TTTTTCACCTGCCCCCCGCGGCGCG‑TG‑AAAAACAAAACGGACAGAAAAA‑ACCGCCGGTATggggccgaattt                                                                                                                                                                                  >  1:502264/1‑218 (MQ=255)
                                            cGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTC‑GGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCCCTACTGGAGCGGCACCAAGAAACG‑CGCCC‑GGATCCACCCGCACCCAGCGGAGCG‑TT‑GAAAGAAAGACGGCAAGAAAAA‑ACCGCCGTTCTGGACCCGCAtttggcggat                                                                                                                                                                           >  2:206869/1‑220 (MQ=255)
                                                 cGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTC‑GGCCTGGTGTTTGCTGGCGGTCTGATTTTTTTGCCGCGTTGCTACTGGGGGGGGACCAAGGAAAG‑CGCCCGGTTTCCCCCTG‑CACCAACGGAGCG‑TG‑CAAAAAAAAACGGCCAGAAAAACCCCGCCG‑ACTGGCCGCGGATTGCGCTGAGCCtgt                                                                                                                                                                      >  1:473422/1‑226 (MQ=255)
                                                  ttCAGTGGTATATGTACTCCATTCCACGAAAGTATTAGTTATTTGCGTTTCTCTGATCGTCTCCCCTGTTGTC‑GGCCTGTTTTTTGCTGGCGGTCTGATTTTCTTGCTGCGGCCCATCTGGGCCGGCACCAAGAAACG‑CGCCG‑GTTTCCACCTGACCCCAGCGGGGCG‑TG‑AAAAAAAAGACGGCAAGAAAAA‑GCCGCCGTTCTGGACGCGTATTGCGCTGATCCttgc                                                                                                                                                                     <  1:256084/227‑3 (MQ=255)
                                                        agggccagCACTCAATACCCGGAAAGTATTAAGGATTTTCGGTTCTCTGATCGTTTCCCCTATTGTC‑GGCGTGGGGTTTGCTGGCGGTCGGATTTTCTTGCTGCGTCGCTACTGGAGCGGACCCAAGAAACG‑CGCCC‑GTACCCACCTGACCCCAGCGGAGCG‑TG‑AAAAGAAAGACGGCAAGAAAAA‑GCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTa                                                                                                                                                                <  2:60072/220‑1 (MQ=255)
                                                                       tatCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTC‑GGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACG‑CGCCC‑GGATCCACCTGACCCCAGCGGAGCG‑GG‑AAAAGAAAGAAGGCAAGAAAAA‑ACCGCCGTTCTGGACGCGACATGCGCTGATCCTTTCGGCTA‑TCGGCGT‑GGGGgtttt                                                                                                                                              >  2:578429/1‑227 (MQ=255)
                                                                          cccGAAAGTTTTAAGGTTTTCGGGTTCTCTGCTCGTTTCCCCTATTGTC‑GGCCTGGTGTTTGCTGGCGGTCTGATTTTCGTGCTGCGTCGCTACGGGAGCGGCACCAAGAAACG‑CGCCC‑GTATCCACCTGACCCCAGCGGAGCG‑TG‑AAAAGAAAGACGGCAAGAAAAA‑GCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGCGTTTTCGCa                                                                                                                                           <  2:247596/228‑1 (MQ=255)
                                                                           cgggatgttttgatttttttcggttttcggACGGTTTCCCCT‑GTGTCGGGCCGGGTGTTTGCTGGCGGTCTTATTTGTTTGATGCGTCGCTGCTGGGGCGGCACCAAAAAACG‑CGCCC‑GTATCCACCTGACCCGAGCGGAGTG‑TA‑AAAAGAAAAACGGCAAGAAAAC‑GCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGCGTTTTCGCAc                                                                                                                                          <  1:367613/212‑1 (MQ=255)
                                                                            cGAAAGTATTAATTATTTTCGGTTCTCTGATCGTTTCCCCTATTGCC‑GGCCTGGTGTTTGCTGGGGGTCTGATTTTCTTGCTGCGTCGCTACGGGAGCGGCACCAAGAAACG‑CGCCC‑GTATCCACCTGACCCCAGCGGAGCG‑TG‑AAAAGAAAGACGGCAAGAAAAA‑GCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGCGTTTTCGCAc                                                                                                                                          <  2:619236/227‑1 (MQ=255)
                                                                               aagtattaagtattTTCGGTTCTCTGATCGTTTCCCCTATTGTC‑GGCCTGGTGTTTGCTGGCGGTCTTGTTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACG‑CGCCC‑GTATCCACCTGACCCCAGCGGAGCG‑TG‑AAAAGAAAGACGGCAAGAACAA‑ACCGCCGTTCTGGACGCGTCTTGAGCTGAACCTTGCCGCTA‑TCGGCGG‑GGCGTTTTTGccccgac                                                                                                                                      >  1:445393/1‑222 (MQ=255)
                                                                               aagtattaagtattTTCGGTTCTCTGATCGTTTCCCCTATTGTC‑GGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACG‑CGCCC‑GTATCCACCTGACCCCAGCGGGGCG‑TG‑AAAAGAAAGACGGCAAGAAAAA‑ACCGCCGTCTTGGCCGCGTTTTGCGCTGATCCTTTCCGCAA‑TCGGCGT‑TGCGGTTTTGCACcgcg                                                                                                                                      >  1:143293/1‑228 (MQ=255)
                                                                                    ttaagtattTTCGGTTCTCTGATCGTTTCCCCTATTGTC‑GTCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTAATGGTGCGGCACCAAGAAACG‑CGCCC‑GTTTTCCCCTGACCCCCGCGGAGCG‑TG‑ATAAGCAAGAAGGCCAGAAAAA‑CACGCCCCTCTGGACCCCGATTGCCCTGGTTCTTTCCGCCACTCGGCGG‑GCCTtttttccccgcgcggac                                                                                                                                  >  1:420012/1‑215 (MQ=255)
                                                                                            tttCGGTTCTCTGATCGTTTCCCCTATTGTC‑GGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTTCGTCGCTACTGGAGCGGCACCAAGAAACG‑CGCCC‑GTATACACCTGACCCCAGCGGAGGG‑TG‑AAAAGAAAGACGGCCAGAAAAA‑ACCGCCGTTCTCGACCCGGATTGCGCTGTTCACTTCCGCCA‑TCGGAGTGGGCGTTTTccccggccgcgaccactgtca                                                                                                                          >  1:409137/1‑208 (MQ=255)
                                                                                            tttCGGTTCTCTGATCGTTTCCCCTATTGTC‑GGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGGGCGGCACCAAGAAACG‑CGCCC‑GTATCCACCTGACCCCAGCGGAGCG‑TG‑AAAAGAAAAACGGCCAGAAAAA‑GCCCCCGTTCTGGACGCGGCTTTGGCCGAACCCTTCCGCCA‑TCGGCCTGGGGGTTTGCGCCCGGGGGGAAC‑GAAggcg                                                                                                                         >  1:132693/1‑226 (MQ=255)
                                                                                            tttCGGTTCTCTGATCGTTTCCCCTATTGTC‑GGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACG‑CGCCC‑GTATCCACCTGACCCCAGCGGAGCG‑TG‑AAAAGAAAGACGGCAAGAAAAA‑GCCGCCGTTCTGGACCCGCATTGCGCTGATCCTTTCCCCTA‑TCGGCGT‑CGCGTTTTCTCACCGGGCGGAC‑ACTGGTcc                                                                                                                        >  2:488332/1‑227 (MQ=255)
                                                                                            tttCGGTTCTCTGATCGTTTCCCCTATTGTC‑GGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACCTGGGCGGCCCCCAGAAACG‑CGCCC‑GTTTCCCCCTTACCCCCCCGGGGCG‑TG‑AAAAGAAAGACGGCCAGAAAAA‑ACCCGCGGTCTGGACACCGATTGCGCGTATCCTTTCCCCCA‑TCGCGGT‑GGGGTTTTTGccgcgcgcccacaagggtca                                                                                                                         >  1:588956/1‑209 (MQ=255)
                                                                                              tcggttctcgttgtgTTTCCCTTAGTGTC‑GGCCTGGTTGTGGCGGGCGGTTCGGTTTGCGTGCTGCTTGCTTACTGGGGCGGCACAAATGAACG‑CGCCT‑GGTCCCACCTGCCCCCAGAGGAGCG‑TA‑AAAAGAAAGACGGCAAGAAAAA‑GCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGCGTTTTCGCACGGCGCGAAC‑GATGGTCAG‑a                                                                                                                     <  1:509163/214‑1 (MQ=19)
                                                                                                   tctcTGATCGTTTCCCCTATTGTC‑GGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACG‑CGCCC‑GTATCCACCTGACCCCCGCGGAGCG‑TG‑GAAAGAAAGACGGCAAGAAAAA‑GCCGCCCGTCTGGACGCGTATTGCGCTGATCCTTTCCCCTA‑TAGGCGT‑GGCGTTTTCTGACgggccggaacaggtgcaaaaaagc                                                                                                                  >  1:457517/1‑206 (MQ=255)
                                                                                                    tctggaccttttCCCCTGTGTGC‑GGCCGGGTGTTTGCTGGCGGTCTGATTTTCGTGCTGCGTCGCTACGGGAGCGGCACAAAGAAACG‑CGCCG‑GTATCCACCTGACCCAAGGGGAGCG‑GG‑AAAAGAAAGCCGGCAAGAAAAA‑GCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGCGTTTTCGCACGGCGCGAAC‑GATGGTCAG‑AAAGGCa                                                                                                               <  1:488332/219‑1 (MQ=37)
                                                                                                      ctGATCGTTTCCCCTATTGTC‑GGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACG‑CGCCC‑GTAACCACCTGACCCCAGCGGAGCG‑TG‑AAAAGAAAGACGGCAAGAAAAA‑GCCGCCGTTCTGGAc                                                                                                                                                                                          <  2:38057/155‑1 (MQ=255)
                                                                                                        gtttgTTTCCCCTTTTGTG‑GGCGTGGTGTTTGCGGGCGGTTGTATTTTCGTGCTGGGCCGCTAGGGGGGCGGCCACAAAAAACG‑CGCCC‑GTATCCACCGGCCCCCAGCGGGGGT‑TG‑AAAAGAAAGACGGCAAGAAAAA‑GCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGCGTTTTCGCACGGCGCGAAC‑GATGGTCAG‑AAAGGCATTgg                                                                                                           <  1:614542/224‑1 (MQ=25)
                                                                                                         tttgTTCCCCCTTTTGTC‑GGGGTTGTGTTTGCTGGCGGTCTGATTTTGTTGCTGCGGCGCTACTGGAGCGGCACAAAGAAACG‑CGCCC‑GTACCCCCCTGCCCCCAGCGGAGTG‑TG‑AAAAGAAAGACGGCAAGAAAAA‑GCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGCGTTTTCGCACGGCGCGAAC‑GATGGTCAG‑AAAGGCATTGGt                                                                                                          <  1:540536/225‑1 (MQ=38)
                                                                                                         tatgttttccctgtttccgggCCGGGTGTTTGCTGGCTGTCTGTTTTTTTTGCTGCGTCCCTACTGGAGCGGCACCAAAAAACG‑CTCCC‑GTACCCACCTGACCCCGGCGGAGCG‑TG‑AAAAGAAAGACGGCAAGAAAAA‑GCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGCGTTTTCGCACGGCGCGAAC‑GATGGTCAG‑AAAGGCATTgg                                                                                                           <  1:131728/209‑1 (MQ=37)
                                                                                                         aTCTTTTCCCCGTTTGGC‑GGCCTGGTTTTGGGTGGCGGCTGGATTTCTTTGCGCGGTCGCTCCTGGGGGGGCACCAGAAAACG‑GGCCG‑ATATCCACCGGACCCCAGGGGGGGT‑GG‑AAAAGAAAAACGGCAAGAAAAG‑GCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGCGTTTTCGCACGGCGCGAAC‑GATGGTCAG‑AAAGGCATTGGt                                                                                                          <  1:206869/228‑1 (MQ=19)
                                                                                                          ttgTTTCCCCTATTGGC‑GGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACG‑CGCCC‑GTATCCACCTGACCCCAGGGGAGCG‑TG‑AAAAGAAAGACGGCAAGAAAAA‑GCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGCGTTTTCGCACGGCGCGAAC‑GATGGTCAG‑AAAGGCATTGGt                                                                                                          <  2:558355/225‑1 (MQ=255)
                                                                                                           atgtgtataagagacgtgcccTGGTGTTTGCTGTCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACG‑CGCCC‑GTATCCACCTGACCCCAGCGGAGCG‑TG‑AAAAGAAAGACGGCAAGAAAAA‑GCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGCGTTTTCGCACGGCGCGAAC‑GATGGTCAG‑AAAGGCATTGGTc                                                                                                         <  2:154064/209‑1 (MQ=255)
                                                                                                               tCCCCTATTGTC‑GGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACG‑CGCCC‑GTATCCACCTGACCCCAGCGGAGCG‑TG‑AAAAGAAAGACGGCAAGAAAAA‑ACCCCCGGTCTGGACGCCGATTGCGCTGATCCTTTCCGCCA‑TCGGCGT‑GGGGGTTTTCCCCGGCGCGACC‑GACGGGCAG‑AAAAGCacttgggtggt                                                                                                     >  1:558355/1‑217 (MQ=255)
                                                                                                                            ggCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACG‑CGCCC‑GTATCCACCTGACCCCAGCGGAGCG‑TG‑AAAAGAAAGACGGCAAGAAAAA‑GCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGGGTTTTTTCACGGCCCGCAC‑GATGGTCCGAAAAGGCTTTGG‑CTGGGTTTGTTGGGattt                                                                                        >  2:204376/1‑227 (MQ=255)
                                                                                                                                        gCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCCAGAAACG‑CGCCC‑GTATCCACCTGACCCCAGCGGAGCG‑GG‑AAAAGAAAGACGGCAAGAAAAA‑GCCGCCGTTCTCGACGCGGATTGCGCTGATCCTTTCCCCTA‑TCGGCGT‑GGGCTTTTTGCCCGGGGCGCACGGAGGGTAAG‑AAA‑GCAAGGGTCGGGGTCTGTTTGTTATGGTTGGCgcccgg                                                                            >  1:569551/1‑223 (MQ=255)
                                                                                                                                              gggttgtttttCCTGTTGCGTCGCTATGGGAGCGGCCACAAGAAGCG‑CGCCC‑GTCTCCCCCTGACCCGAGGGGGGGG‑TG‑AAAAGAAAGACGGCAAGAAAAG‑GCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGCGTTTTCGCACGGCGCGAAC‑GATGGTCAG‑AAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCcagc                                                                       <  1:418567/220‑1 (MQ=37)
                                                                                                                                                     gtttgtgggcgcgcggacaggggAGGGGCACCAAGAAACG‑GCCCG‑GAAACCAGCTGACCCCAGCGGAGGG‑AA‑AAAAGAAGGTAGTGAAGAAAAG‑GCCGCCGGTCGGGACGAGTATTGCGCTGAGCCTTGCCGGGA‑TCGGCGG‑GGCGTTTTCGCATGGCGCGAAC‑GATGGTCAG‑AAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCg                                                               <  1:461648/207‑1 (MQ=255)
                                                                                                                                                           ttctggggtgCTACGGGGGCGGCAACAAGAAAGG‑CGCCC‑TGAACCACCTGACCCCAG‑GGGGCGTTG‑AAAAGAAATACGGCAAGAAAAA‑GGCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGCGTTTTCGCACGGCGCGAAC‑GATGGTCAG‑AAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCgtggtg                                                          <  1:34107/218‑1 (MQ=37)
                                                                                                                                                                         tGGAGCGGTCCCAAGAAACGCCGCCT‑GT‑CCCACTCGACCCCAGCGGAGCG‑TG‑AAAAGAAAGACGGCAAGAAAAA‑GCCGCCTTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGCGTTTTCGCACGGCGCGAAC‑GATGGTCAG‑AAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACt                                           <  2:454136/228‑1 (MQ=37)
                                                                                                                                                                             ggcgaaacAAAAACCG‑CCCCC‑GTAACCCCCGGACCCC‑GGGGGGCG‑GGAAAAAGAAAGACGGCAAAAAAAA‑GCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGCGTTTTCGCACGGCGCGAAC‑GATGGTCAG‑AAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCt                                       <  1:204376/221‑1 (MQ=25)
                                                                                                                                                                              gggCACAAAGAAACG‑CCCCC‑GTAACCCCCTGACCCCGGCGGAGGT‑GG‑AAAAGAAAGGCGGCAAGAAAAG‑GCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTA‑GGGGAGT‑GGCGTTTTCGCACGGCGCGAAC‑GATGGTCAG‑AAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTa                                      <  2:588956/227‑1 (MQ=255)
                                                                                                                                                                               ggCACCAAGAAACG‑CGCCC‑GTATCCACCTGACCCCAGCGGAGCG‑TG‑AAAAGAAAGACGGCAAGAAAAA‑GCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGCGTTTTCGCACGGCGCGAAC‑GATGGTCAG‑AAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCCTCCTGGTGAACATGAATGCCACTGGCTAc                                     >  2:119258/1‑228 (MQ=255)
                                                                                                                                                                                  aCCAAGAAACG‑CGCCC‑GTAACCACCGGACCCCAGCGGAGCG‑TG‑AAAAAAAAGACGGCAAGAAAAA‑GCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGCGTTTTCGCACGGCGCGAAC‑GATGGTCAG‑AAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGaa                                  <  2:445393/228‑1 (MQ=255)
                                                                                                                                                                                     aaGAAACG‑CGCCC‑GTATCCACCTGGCCCCGGGGGGGTG‑TG‑AAAAGAAAGACGGCAAGAAAAA‑GCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGCGTTTTCGCACGGCGCGAAC‑GATGGTCAG‑AAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATc                               <  2:601671/228‑1 (MQ=255)
                                                                                                                                                                                     aaGAAACG‑CGCCC‑GTATCCACCTGCCCCCAGGGGAGCG‑TG‑AAAAAAAAGACGGCAAGAAAAA‑TCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGCGTTTTCGCACGGCGCGAAC‑GATGGTCAG‑AAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATc                               <  2:497587/228‑1 (MQ=255)
                                                                                                                                                                                      aaaaaaggacccc‑GGATCCAACTGGCCGCATGGGAGCG‑GG‑AAAAGAAAGACGGAAAAAAAAA‑GCCGCCGTTCGGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGCGTTTTCGCACGGCGCGAAC‑GATGGTCAG‑AAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATc                               <  2:596543/218‑1 (MQ=37)
                                                                                                                                                                                             ccccgcggtttacctgtggccccCGCGGAGGG‑TA‑AAAAGAAAGACGGCAAAAAAAA‑GCCGCCGTTCCGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGCGTTTTCGCACGGCGCGAAC‑GATGGTCAG‑AAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCacccgt                         <  1:622560/209‑1 (MQ=37)
                                                                                                                                                                                              ccccc‑GTACCCCCCGAACCCAGGGGGGGCG‑TG‑AAAAAAAAGCCGGAAAGAAAAG‑GCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGCGTTTTCGCACGGCGCGAAC‑GATGGTCAG‑AAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCacccgtac                       <  2:420012/226‑1 (MQ=37)
                                                                                                                                                                                                       tCCACCTGACCCGAGCGGAGCG‑TG‑AAAAGAAAGCCGGCAAGAAAAA‑GCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGCGTTTTCGCACGGCGCGAAC‑GATGGTCAG‑AAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATg               <  2:132693/228‑1 (MQ=255)
                                                                                                                                                                                                          aCCTGACCCCAGCGGAGCG‑TG‑AAAAGAAAGACGGCAAGAAAAA‑GCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGCGTTTTCGCACGGCGCGAAC‑GATGGTCAG‑AAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCCGGGTTCCGGGGGAAAATGAAAGCCCCATGGTACCAAATACACCGCACCCGGGATGcac            >  1:557752/1‑226 (MQ=37)
                                                                                                                                                                                                            cTGACCCCAGCGGAGCG‑TG‑AAAAGAAAGCCGGCAAGACAAA‑GCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGCAGTTTCGCACGGCGCGCAC‑GGTGGTCAA‑AACGGCATTGGTCTGGGTATGTTGGGCTTTGTTGGGGTCACGCCAGCGGGCTTCGTGGTGAACATGAACGCCACTGGGCAGGGAATCacacggcccctggagtacaca          >  1:99864/1‑212 (MQ=255)
                                                                                                                                                                                                            cTGACCCCAGCGGAGCG‑TG‑AAAAGAAAGACGGCAAGGAAAA‑GCCGCCGTTCTGGa                                                                                                                                                                                           >  1:328001/1‑55 (MQ=255)
                                                                                                                                                                                                               aCCCCAGCGGAGCG‑TG‑GAAAGAAAGACGGCAAGAAAAA‑GCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGCGTTTTCGCACCGCGCGAAC‑AGTGGTCAG‑AAAGGCCTTGGGCTGGTTATGTTGGTATTGATTTGCGTCCCGCCAGCAGGCTTCGTGGGGAACATGAATGCACCTGGCTTCGAAAACAACCGGACACCGGATGGCCTcaaa       >  1:249235/1‑227 (MQ=255)
                                                                                                                                                                                                                    aGCGGAGCG‑TG‑AAAAGAAAGACGGCAAGAAAAA‑GCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGCGTTTTCGCACGGCGCGAAC‑GATGGTCAG‑AAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCCCGCCAGCAGGCTTCGTGGGGAACATGAATGCCACTGGCTACGCAATCACCCCGACCCGTGATGCCATCAACAACGt  >  2:606356/1‑228 (MQ=255)
                                                                                                                                                                                                                         |                                                                                                                                                                                                                                     
ATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTC‑GGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACG‑CGCCC‑GTATCCACCTGACCCCAGCGGAGCG‑TG‑AAAAGAAAGACGGCAAGAAAAA‑GCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTA‑TCGGCGT‑GGCGTTTTCGCACGGCGCGAAC‑GATGGTCAG‑AAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGT  >  NC_000913/3638013‑3638449

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