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
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  
gcacatcaaccacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAAt                                                                                                                                                                                                                                             <  1:223328‑M2/32‑1 (MQ=255)
    atcaacaacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTa                                                                                                                                                                                                                                         >  1:25488‑M2/215‑250 (MQ=255)
      caaccacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAg                                                                                                                                                                                                                                      <  1:467323‑M2/39‑1 (MQ=255)
           acgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTcc                                                                                                                                                                                                                                 >  2:494169‑M2/208‑251 (MQ=255)
             gacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTc                                                                                                                                                                                                                               <  2:239165‑M2/46‑1 (MQ=255)
               cgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTc                                                                                                                                                                                                                               >  1:15975‑M2/204‑249 (MQ=255)
                    gtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTc                                                                                                                                                                                                                        <  1:408353‑M2/53‑1 (MQ=255)
                              cgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCAt                                                                                                                                                                                                              >  1:419025‑M2/189‑251 (MQ=255)
                               gccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATc                                                                                                                                                                                                             <  1:38392‑M2/64‑1 (MQ=255)
                                gcgttttgaaaggtgacatgggttcggtgcagccttttgccgatgacggtatcgcccaggccggcgtgcggctcgggcatttgccgtgcttgccgacggaagattgagtctattaatgggtcctttggtgaacgagcaacaaagggccactagctcacctgattcctctcgaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGCCGTGAATTGTAAAGATTCCTCAGTGGTCGCACCTGCATCa                                                                                                                                                                                                            <  2:5403‑M2/65‑1 (MQ=255)
                                ccgtgtcgaaaggtgacatgggtgccgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAACAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCCCCAGCATc                                                                                                                                                                                                             >  2:405793‑M2/187‑250 (MQ=255)
                                    gttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAatat                                                                                                                                                                                                        <  1:78403‑M2/69‑1 (MQ=255)
                                         aaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACt                                                                                                                                                                                                    >  1:249209‑M2/178‑250 (MQ=255)
                                         aaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTa                                                                                                                                                                                                   <  1:327395‑M2/74‑1 (MQ=255)
                                             tgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTaaa                                                                                                                                                                                                 >  2:133419‑M2/174‑249 (MQ=255)
                                                catgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCGCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATTAATATTACTAAAAGGaa                                                                                                                                                                                            <  1:384374‑M2/81‑1 (MQ=255)
                                                          gtgcagcagtttgccgatgacgttatcgctcagcggggcgggcggcacgggcatttgcagtgcgtgccgaatgaagattgagtctattattgggtcctttggtgaacgatcaacaatgggccactagcacccctgattccgctctaaattccctatccgaGTTGCGTTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATGCATCAGTGGTCGCACCTGCCTCAATATTACTAAAAGTAAGTATTTGCCt                                                                                                                                                                                  <  2:53585‑M2/91‑1 (MQ=255)
                                                                 attttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagatttagtctatgagtgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctcgctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATtat                                                                                                                                                                           <  2:17818‑M2/98‑1 (MQ=255)
                                                                            gacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAat                                                                                                                                                                 <  1:83949‑M2/108‑1 (MQ=255)
                                                                            gacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAat                                                                                                                                                                 <  1:420035‑M2/108‑1 (MQ=255)
                                                                                            gcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCATCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTGTTGAGTGAAAt                                                                                                                                                <  1:173114‑M2/125‑1 (MQ=255)
                                                                                                cgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCAt                                                                                                                                            <  1:45242‑M2/129‑1 (MQ=255)
                                                                                                     ggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAgg                                                                                                                                       <  1:107594‑M2/134‑1 (MQ=255)
                                                                                                       cacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAgggg                                                                                                                                     <  1:343497‑M2/136‑1 (MQ=255)
                                                                                                          gggcatttgcagtgcttgccgaaggaagattgcgtctattattgggtcctgttgtgaacgatcaacaaagggccactagcacacctgattcctctctaaattccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGAACCAGCAACAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCa                                                                                                                                  <  2:59611‑M2/139‑1 (MQ=255)
                                                                                                                    agtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCaaaaaa                                                                                                                        <  2:375456‑M2/149‑1 (MQ=255)
                                                                                                                           gccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaaggtccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGTACGTTAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCtat                                                                                                                 <  2:278576‑M2/156‑1 (MQ=255)
                                                                                                                              gaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGCTTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGTCAAATCAAAAAAAGTCTATAtt                                                                                                              >  2:371760‑M2/93‑251 (MQ=255)
                                                                                                                              gaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATAtt                                                                                                              <  2:484991‑M2/159‑1 (MQ=255)
                                                                                                                                 ggaagattgcgtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTc                                                                                                            >  2:233829‑M2/90‑250 (MQ=255)
                                                                                                                                        tgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGc                                                                                                     <  1:286857‑M2/168‑1 (MQ=255)
                                                                                                                                               attattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCccgcgcc                                                                                              <  1:133419‑M2/175‑1 (MQ=255)
                                                                                                                                                    tgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGaa                                                                                         >  1:10172‑M2/71‑250 (MQ=255)
                                                                                                                                                                    cgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATgcgc                                                                        >  2:55395‑M2/55‑251 (MQ=255)
                                                                                                                                                                         aacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGt                                                                   <  1:278487‑M2/202‑1 (MQ=255)
                                                                                                                                                                                  ccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGACAGTCACTGATAATGCGCCGCGTTCATGTCCt                                                          >  2:42972‑M2/41‑251 (MQ=255)
                                                                                                                                                                                    actagcacagctgattcctctctaaataccctatcctaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGCGGTACCACCAGCATCAATATTACTAAAAGGAAGTATTGGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCTATACAGGGGGCAAATCAAACAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTc                                                         >  1:124682‑M2/39‑250 (MQ=255)
                                                                                                                                                                                        gcacacctgattcctctctaaataccctagccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAAt                                                    <  2:162291‑M2/217‑1 (MQ=255)
                                                                                                                                                                                           cacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGc                                                 >  1:4717‑M2/32‑251 (MQ=255)
                                                                                                                                                                                            acctgattcctctctaaataccctatcagaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAACGATTCCTCAGTGGTCGCACACGCATCAATATTACTAAAAGGCAGTCTTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATACATACAGGGGGCAACTCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGAGTTCATGTCCTCAAAATGGCg                                                >  2:273457‑M2/31‑251 (MQ=255)
                                                                                                                                                                                                gattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAAc                                            >  1:234552‑M2/27‑251 (MQ=255)
                                                                                                                                                                                                           aaataccctatccgaGTTGCGCTCCTCTTACAAAAACGGACGTGAATTGTAAAGATTCCTCCGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAGGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCt                                  >  2:314073‑M2/16‑250 (MQ=255)
                                                                                                                                                                                                             ttaccctatacgaGTTACGCTCCTCTTAAACAAACGGACGTGAATTGTAAAGATTCCTCAGCGGTCTATCCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATACGATACTTATTTTTTGAGTGCAATCCATACAGCGGACAAATCAAAAAAAGCCTATATTTCACTTTGCCCGCGCCGCGCAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGATAc                                >  2:128307‑M2/14‑250 (MQ=255)
                                                                                                                                                                                                                       cgaGTAGCGCTCCTCTTAAAAAAACGTACGTGAATTGTAAAGATTGCTCAGTGGTCGCTCCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTGTTTTAGTGAAAGCCGTACAGCGGGCAAATCAAAAAAAGTCTATACTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCGATGCTACATTTTTTTGCt                     <  2:67440‑M2/248‑1 (MQ=255)
                                                                                                                                                                                                                        gaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCCGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAATAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGAAACGTCCTATGCTACATTTGTTTGCTg                    >  2:237075‑M2/3‑251 (MQ=255)
                                                                                                                                                                                                                         aGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCCGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTgg                   >  2:435578‑M2/2‑251 (MQ=255)
                                                                                                                                                                                                                             gcgcTCATCTTAAAAAAACGGACGTTAATTGTAAAGATTCCTCAGTGGTAGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCTAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTg               <  2:505401/251‑1 (MQ=255)
                                                                                                                                                                                                                                  cctcTTAAAAAAACGGACGTGAATTGTAAAGCGTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGGCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGGCCGCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGAttt          <  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: