Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,637,436 G→A intergenic (‑26/‑206) yhiN ← / → pitA putative oxidoreductase YhiN/metal phosphate:H(+) symporter PitA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,637,4360GA98.0% 157.0 / NA 50intergenic (‑26/‑206)yhiN/pitAputative oxidoreductase YhiN/metal phosphate:H(+) symporter PitA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (21/28);  minor base C (0/1);  total (21/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAGGTTGGTAAAGTTGCAGCGCCCACCGCCAGACATAAGGATTTTGCGCCCTGGTTTTTTACCATTATCGATCAGCAGAACCCGGCGTCCTGCCTGACCTGCCAGCGCAGAACAGAACATACCCGCCGCACCAGCGCCTATAATAATGGCATCAAACCTTTCCACGTTGCGCTCCTCTTAGAAAAAACGGGCGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCG  >  NC_000913/3637193‑3637678
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  
gcacatcaaccacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAAt                                                                                                                                                                                                                                             <  1:223328‑M2/32‑1 (MQ=255)
    atcaacaacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTa                                                                                                                                                                                                                                         >  1:25488‑M2/215‑250 (MQ=255)
      caaccacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAg                                                                                                                                                                                                                                      <  1:467323‑M2/39‑1 (MQ=255)
           acgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTcc                                                                                                                                                                                                                                 >  2:494169‑M2/208‑251 (MQ=255)
             gacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTc                                                                                                                                                                                                                               <  2:239165‑M2/46‑1 (MQ=255)
               cgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTc                                                                                                                                                                                                                               >  1:15975‑M2/204‑249 (MQ=255)
                    gtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTc                                                                                                                                                                                                                        <  1:408353‑M2/53‑1 (MQ=255)
                              cgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCAt                                                                                                                                                                                                              >  1:419025‑M2/189‑251 (MQ=255)
                               gccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATc                                                                                                                                                                                                             <  1:38392‑M2/64‑1 (MQ=255)
                                gcgttttgaaaggtgacatgggttcggtgcagccttttgccgatgacggtatcgcccaggccggcgtgcggctcgggcatttgccgtgcttgccgacggaagattgagtctattaatgggtcctttggtgaacgagcaacaaagggccactagctcacctgattcctctcgaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGCCGTGAATTGTAAAGATTCCTCAGTGGTCGCACCTGCATCa                                                                                                                                                                                                            <  2:5403‑M2/65‑1 (MQ=255)
                                ccgtgtcgaaaggtgacatgggtgccgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AACAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCCCCAGCATc                                                                                                                                                                                                             >  2:405793‑M2/187‑250 (MQ=255)
                                    gttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAatat                                                                                                                                                                                                        <  1:78403‑M2/69‑1 (MQ=255)
                                         aaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACt                                                                                                                                                                                                    >  1:249209‑M2/178‑250 (MQ=255)
                                         aaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTa                                                                                                                                                                                                   <  1:327395‑M2/74‑1 (MQ=255)
                                             tgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTaaa                                                                                                                                                                                                 >  2:133419‑M2/174‑249 (MQ=255)
                                                catgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCGCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATTAATATTACTAAAAGGaa                                                                                                                                                                                            <  1:384374‑M2/81‑1 (MQ=255)
                                                          gtgcagcagtttgccgatgacgttatcgctcagcggggcgggcggcacgggcatttgcagtgcgtgccgaatgaagattgagtctattattgggtcctttggtgaacgatcaacaatgggccactagcacccctgattccgctctaaattccctatccgaGTTGCGTTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATGCATCAGTGGTCGCACCTGCCTCAATATTACTAAAAGTAAGTATTTGCCt                                                                                                                                                                                  <  2:53585‑M2/91‑1 (MQ=255)
                                                                 attttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagatttagtctatgagtgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctcgctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATtat                                                                                                                                                                           <  2:17818‑M2/98‑1 (MQ=255)
                                                                            gacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAat                                                                                                                                                                 <  1:83949‑M2/108‑1 (MQ=255)
                                                                            gacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAat                                                                                                                                                                 <  1:420035‑M2/108‑1 (MQ=255)
                                                                                            gcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCATCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTGTTGAGTGAAAt                                                                                                                                                <  1:173114‑M2/125‑1 (MQ=255)
                                                                                                cgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCAt                                                                                                                                            <  1:45242‑M2/129‑1 (MQ=255)
                                                                                                     ggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAgg                                                                                                                                       <  1:107594‑M2/134‑1 (MQ=255)
                                                                                                       cacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAgggg                                                                                                                                     <  1:343497‑M2/136‑1 (MQ=255)
                                                                                                          gggcatttgcagtgcttgccgaaggaagattgcgtctattattgggtcctgttgtgaacgatcaacaaagggccactagcacacctgattcctctctaaattccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGAACCAGCAACAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCa                                                                                                                                  <  2:59611‑M2/139‑1 (MQ=255)
                                                                                                                    agtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCaaaaaa                                                                                                                        <  2:375456‑M2/149‑1 (MQ=255)
                                                                                                                           gccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaaggtccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGTACGTTAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCtat                                                                                                                 <  2:278576‑M2/156‑1 (MQ=255)
                                                                                                                              gaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGCTTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGTCAAATCAAAAAAAGTCTATAtt                                                                                                              >  2:371760‑M2/93‑251 (MQ=255)
                                                                                                                              gaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATAtt                                                                                                              <  2:484991‑M2/159‑1 (MQ=255)
                                                                                                                                 ggaagattgcgtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTc                                                                                                            >  2:233829‑M2/90‑250 (MQ=255)
                                                                                                                                        tgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGc                                                                                                     <  1:286857‑M2/168‑1 (MQ=255)
                                                                                                                                               attattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCccgcgcc                                                                                              <  1:133419‑M2/175‑1 (MQ=255)
                                                                                                                                                    tgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGaa                                                                                         >  1:10172‑M2/71‑250 (MQ=255)
                                                                                                                                                                    cgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATgcgc                                                                        >  2:55395‑M2/55‑251 (MQ=255)
                                                                                                                                                                         aacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGt                                                                   <  1:278487‑M2/202‑1 (MQ=255)
                                                                                                                                                                                  ccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGACAGTCACTGATAATGCGCCGCGTTCATGTCCt                                                          >  2:42972‑M2/41‑251 (MQ=255)
                                                                                                                                                                                    actagcacagctgattcctctctaaataccctatcctaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGCGGTACCACCAGCATCAATATTACTAAAAGGAAGTATTGGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCTATACAGGGGGCAAATCAAACAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTc                                                         >  1:124682‑M2/39‑250 (MQ=255)
                                                                                                                                                                                        gcacacctgattcctctctaaataccctagccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAAt                                                    <  2:162291‑M2/217‑1 (MQ=255)
                                                                                                                                                                                           cacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGc                                                 >  1:4717‑M2/32‑251 (MQ=255)
                                                                                                                                                                                            acctgattcctctctaaataccctatcagaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAACGATTCCTCAGTGGTCGCACACGCATCAATATTACTAAAAGGCAGTCTTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATACATACAGGGGGCAACTCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGAGTTCATGTCCTCAAAATGGCg                                                >  2:273457‑M2/31‑251 (MQ=255)
                                                                                                                                                                                                gattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAAc                                            >  1:234552‑M2/27‑251 (MQ=255)
                                                                                                                                                                                                           aaataccctatccgaGTTGCGCTCCTCTTA‑CAAAAACGGACGTGAATTGTAAAGATTCCTCCGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAGGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCt                                  >  2:314073‑M2/16‑250 (MQ=255)
                                                                                                                                                                                                             ttaccctatacgaGTTACGCTCCTCTTA‑AACAAACGGACGTGAATTGTAAAGATTCCTCAGCGGTCTATCCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATACGATACTTATTTTTTGAGTGCAATCCATACAGCGGACAAATCAAAAAAAGCCTATATTTCACTTTGCCCGCGCCGCGCAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGATAc                                >  2:128307‑M2/14‑250 (MQ=255)
                                                                                                                                                                                                                       cgaGTAGCGCTCCTCTTA‑AAAAAACGTACGTGAATTGTAAAGATTGCTCAGTGGTCGCTCCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTGTTTTAGTGAAAGCCGTACAGCGGGCAAATCAAAAAAAGTCTATACTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCGATGCTACATTTTTTTGCt                     <  2:67440‑M2/248‑1 (MQ=255)
                                                                                                                                                                                                                        gaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCCGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAATAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGAAACGTCCTATGCTACATTTGTTTGCTg                    >  2:237075‑M2/3‑251 (MQ=255)
                                                                                                                                                                                                                         aGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCCGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTgg                   >  2:435578‑M2/2‑251 (MQ=255)
                                                                                                                                                                                                                             gcgcTCATCTTA‑AAAAAACGGACGTTAATTGTAAAGATTCCTCAGTGGTAGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCTAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTg               <  2:505401/251‑1 (MQ=255)
                                                                                                                                                                                                                                  cctcTTA‑AAAAAACGGACGTGAATTGTAAAGCGTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGGCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGGCCGCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGAttt          <  2:255068/251‑1 (MQ=255)
                                                                                                                                                                                                                                         aaaaaaaCGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATAc    <  1:480683/250‑1 (MQ=255)
                                                                                                                                                                                                                                           aaaaaCGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCg  >  1:415268/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  
GCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAGGTTGGTAAAGTTGCAGCGCCCACCGCCAGACATAAGGATTTTGCGCCCTGGTTTTTTACCATTATCGATCAGCAGAACCCGGCGTCCTGCCTGACCTGCCAGCGCAGAACAGAACATACCCGCCGCACCAGCGCCTATAATAATGGCATCAAACCTTTCCACGTTGCGCTCCTCTTAGAAAAAACGGGCGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCG  >  NC_000913/3637193‑3637678

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