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,4360GA96.5% 199.8 / NA 59intergenic (‑26/‑206)yhiN/pitAputative oxidoreductase YhiN/metal phosphate:H(+) symporter PitA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/38);  minor base C (1/1);  total (20/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAGGTTGGTAAAGTTGCAGCGCCCACCGCCAGACATAAGGATTTTGCGCCCTGGTTTTTTACCATTATCGATCAGCAGAACCCGGCGTCCTGCCTGACCTGCCAGCGCAGAACAGAACATACCCGCCGCACCAGCGCCTATAATAATGGCATCAAACCTTTCCACGTTGCGCTCCTCTTAGAAAAAACGGGCGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGAT  >  NC_000913/3637188‑3637668
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                        
catgtgcacatcaaccacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGt                                                                                                                                                                                                                                       <  1:550880‑M2/28‑1 (MQ=255)
   gtgcaaatcaaccacgaaggattgctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacccctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGa                                                                                                                                                                                                                                     <  1:581116‑M2/30‑1 (MQ=255)
     gcacatcaaccacgacgactgtctggaaagcgccgtgttgaaaggggacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggaatttgcagtgcttgccgaaggaagattgggtctattattgggtcctttggtgaacgatcaccaaagggccactagcccacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGAcgtgcat                                                                                                                                                                                                                                   <  2:377064‑M2/32‑4 (MQ=255)
        catcaaccacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTa                                                                                                                                                                                                                               <  1:392292‑M2/36‑1 (MQ=255)
             accacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGAt                                                                                                                                                                                                                          <  1:155773‑M2/41‑1 (MQ=255)
                 cgacgactgtctggaaatcgccgtgttgaaangtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCt                                                                                                                                                                                                                      <  2:272305‑M2/45‑1 (MQ=255)
                 cgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCt                                                                                                                                                                                                                      <  1:409819‑M2/45‑1 (MQ=255)
                  gacgattgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgtcgttatcgcccagcgcggcgtgcggcactggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggctactagcacaccagatgcctcgctaaataccctatccgaGGTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTc                                                                                                                                                                                                                     <  2:341087‑M2/46‑1 (MQ=255)
                  gacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggacgattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTT‑CAAAAAACGGCCGTGAATTGTAACAATTCCTc                                                                                                                                                                                                                     >  2:214033‑M2/206‑251 (MQ=255)
                  gacgactgtctggaaagcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTc                                                                                                                                                                                                                     <  2:245582‑M2/46‑1 (MQ=255)
                    cgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAg                                                                                                                                                                                                                   <  2:326409‑M2/48‑1 (MQ=255)
                            tggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcatttggccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGc                                                                                                                                                                                                            <  2:407542‑M2/55‑1 (MQ=255)
                            tggaaatcgccgtgttgaaaggtgacatgggggacgtgcagcattttgccgatgacgtaatcgcccagcgcggcgggcggcacgggcatttgcagtgcttgccgacggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctcactaaataccctctccgaGTTGCGCTCCTCTTA‑AAAAAACGATCGAGAATTGTAAAGAttcatcaggtgacga                                                                                                                                                                                                            >  2:551276‑M2/196‑243 (MQ=255)
                             ggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGc                                                                                                                                                                                                            <  1:184840‑M2/55‑1 (MQ=255)
                                    gccgtgttgaaaggtgacatgggtgacgtgcagcatgttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATc                                                                                                                                                                                                   <  2:239992‑M2/64‑1 (MQ=255)
                                    gccgtgttgaaaggtgaaatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATc                                                                                                                                                                                                   >  2:246583‑M2/188‑251 (MQ=255)
                                     ccgtgttgaaaggggacatgggagacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATACCTCAGTAGTAGCACCAGCATc                                                                                                                                                                                                   >  2:118382‑M2/187‑250 (MQ=255)
                                         ggtgaaaggtgacaggggtgacgtgcagcattttgccgaagacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagatggagtctattagggggtcctttggggaaagatcaacaaagggccactagcacacctgattcctctctaaatcccctatccgaGTTTCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGAGTCCTCAGTGGACGCACCAGCATCAatat                                                                                                                                                                                              <  2:30967‑M2/69‑1 (MQ=255)
                                                ggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTaaa                                                                                                                                                                                       <  1:511487‑M2/76‑1 (MQ=255)
                                                ggggacttgggtgacgtgcagcattttgccgatgacgttagcgcccagctcggcgggcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTaaa                                                                                                                                                                                       <  2:199879‑M2/76‑1 (MQ=255)
                                                 gtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTaaaa                                                                                                                                                                                      <  1:160145‑M2/77‑1 (MQ=255)
                                                   gacatgggtgacgtgcagcattttgccgatgacgctatcgcccagcgcggcgtgcggcacgggcctttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCCGTGGTCGCACCAGCATCAATTTTACTAAAAgg                                                                                                                                                                                    >  2:93045‑M2/173‑251 (MQ=255)
                                                   gacatgggggacgggcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGAATCAATATTACTAAAAgg                                                                                                                                                                                    <  2:477713‑M2/79‑1 (MQ=255)
                                                    acatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGa                                                                                                                                                                                   <  1:470088‑M2/80‑1 (MQ=255)
                                                      atgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAg                                                                                                                                                                                 <  1:381721‑M2/82‑1 (MQ=255)
                                                        gggtgacgtgcagcattttgccgatgacgttatcgaccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGt                                                                                                                                                                                >  1:159004‑M2/168‑250 (MQ=255)
                                                          gtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTAtt                                                                                                                                                                             >  1:507917‑M2/166‑251 (MQ=255)
                                                            gacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTg                                                                                                                                                                           <  1:288377‑M2/88‑1 (MQ=255)
                                                                      attttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATtat                                                                                                                                                                 >  1:363604‑M2/154‑251 (MQ=255)
                                                                             cgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGAt                                                                                                                                                          <  1:295324‑M2/105‑1 (MQ=255)
                                                                              gatgacgttatcgcccaccgcggcgtgaggcacgggcatttgcagtgcttgccgaaggcagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGAt                                                                                                                                                          >  1:450890‑M2/146‑250 (MQ=255)
                                                                               atgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggcagattgagtctataattgggtcctttggtgaacgatcaaaaaagggccactcgcacacctgattcctctctcaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATCAGATa                                                                                                                                                         >  1:322876‑M2/145‑250 (MQ=255)
                                                                                 gacgttatcgcccagcgcggcgtgcggcacgngcatttgcagtgcttgcctaaggaagatttatgctattatggagtccttgggtgaacgatctacaaaggggcactatctctcctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCATTGGTCTCTCCAGCATCAATATTACTAAAAGGAAGTATTTGCCTAAATTATATAAGATAatt                                                                                                                                                      <  2:273592‑M2/109‑1 (MQ=255)
                                                                                     ttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTAttt                                                                                                                                                  >  1:128742‑M2/139‑251 (MQ=255)
                                                                                         cgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTAtttttt                                                                                                                                               <  1:122774‑M2/116‑1 (MQ=255)
                                                                                         cgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTAtttttt                                                                                                                                               <  1:509096‑M2/116‑1 (MQ=255)
                                                                                          gcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGa                                                                                                                                             <  1:39333‑M2/118‑1 (MQ=255)
                                                                                            ccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGt                                                                                                                                           <  1:539046‑M2/120‑1 (MQ=255)
                                                                                                     cgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTGTTTGAGTGAAATCCAt                                                                                                                                  <  1:411344‑M2/129‑1 (MQ=255)
                                                                                                      gtgcggcactggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacccctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAAGTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATa                                                                                                                                 <  2:127092‑M2/130‑1 (MQ=255)
                                                                                                                              ttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCt                                                                                                         >  1:79849‑M2/98‑251 (MQ=255)
                                                                                                                                   gaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTTCTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATAtt                                                                                                    <  1:487288‑M2/159‑1 (MQ=255)
                                                                                                                                      ggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCa                                                                                                 <  1:198185‑M2/162‑1 (MQ=255)
                                                                                                                                                    attattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCccgcgccg                                                                                   >  1:567182‑M2/76‑251 (MQ=255)
                                                                                                                                                        ttgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccaaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCGTCAATATTACTAAAAGGa                                                                                                                                                                                   <  2:493825‑M2/80‑1 (MQ=255)
                                                                                                                                                        ttgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccaaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGa                                                                                                                                                                                   >  1:493825‑M2/72‑151 (MQ=255)
                                                                                                                                                                      gaacgatcaacaaagggccactagcacacctgattcctctctaaataccctaaccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATg                                                                 <  1:214033‑M2/194‑1 (MQ=255)
                                                                                                                                                                           atcaacaaagggccactagcacacctgattcnccgctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAATGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGcgccg                                                            <  2:290357‑M2/199‑1 (MQ=255)
                                                                                                                                                                           atcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGcgccg                                                            >  1:254161‑M2/53‑251 (MQ=255)
                                                                                                                                                                              aacaacgggccactagctcatctgattcctctctaaatcccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATTTCAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAACAAAGTCTATATTTCACTTTGCCTGCGCCGCGAAAGTCACTGATAATGCGCCGCGt                                                         >  2:42967‑M2/50‑251 (MQ=255)
                                                                                                                                                                                 aaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTc                                                       >  1:432426‑M2/47‑250 (MQ=255)
                                                                                                                                                                                            agcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCaaa                                            <  1:573850‑M2/215‑1 (MQ=255)
                                                                                                                                                                                              cacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATg                                         <  2:377617‑M2/218‑1 (MQ=255)
                                                                                                                                                                                               ccacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATgg                                        <  1:13212‑M2/219‑1 (MQ=255)
                                                                                                                                                                                                              ctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGc                         >  2:541005‑M2/18‑251 (MQ=255)
                                                                                                                                                                                                                           cagagttaagcTCATCTTA‑AAAAAACGGACGAGAATTGTCCAGATTCGTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGAATTTGCCTGAATTATATAGGATAATTATCTGTTGAGTTAAATAAATCCTGGGTGCAAATAAAATAAAGCTTATCTTTCACTTTCCCCGCGCCGAGACAGTACCTGATAATGCGCCACGTTCATGTCCTCACAATGGAGTAACGACCTATGCTACActtttttg             >  2:436639/10‑250 (MQ=255)
                                                                                                                                                                                                                            cgaGTTGCGCTACTCTTA‑AAAAAACGGACGTGAATTGTACAGATTCCTCAGTGGTCGCCCCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATACATACCGGGTGCAAATCAAAAAAAGTCTATTTTTCACTTTGCCCGCTCCGCGAAAGTCACTGATAATGCGCCCCGTTCATGTCCTCACAANTGCGGAACGCCCTATGCTACCTTTGTTTGCt           >  2:279237‑M2/4‑251 (MQ=255)
                                                                                                                                                                                                                                   cgcTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCTTATGCTACATTTGTTTGCTGGCCTgg    <  1:511199/251‑1 (MQ=255)
                                                                                                                                                                                                                                     cgcctcTTA‑AAAAAACTGCCGTGAATTGTAAAGATTCCTCAGTGTTCACACCAGCATTAATATTACTAAAAGGAAATATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCAGACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCTTAACGGCCTATGCGACAATTGTTTGCTGGACTgtat  <  2:525432/249‑4 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                        
CGAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAGGTTGGTAAAGTTGCAGCGCCCACCGCCAGACATAAGGATTTTGCGCCCTGGTTTTTTACCATTATCGATCAGCAGAACCCGGCGTCCTGCCTGACCTGCCAGCGCAGAACAGAACATACCCGCCGCACCAGCGCCTATAATAATGGCATCAAACCTTTCCACGTTGCGCTCCTCTTAGAAAAAACGGGCGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGAT  >  NC_000913/3637188‑3637668

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