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
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                        
catgtgcacatcaaccacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGt                                                                                                                                                                                                                                       <  1:550880‑M2/28‑1 (MQ=255)
   gtgcaaatcaaccacgaaggattgctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacccctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGa                                                                                                                                                                                                                                     <  1:581116‑M2/30‑1 (MQ=255)
     gcacatcaaccacgacgactgtctggaaagcgccgtgttgaaaggggacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggaatttgcagtgcttgccgaaggaagattgggtctattattgggtcctttggtgaacgatcaccaaagggccactagcccacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGAcgtgcat                                                                                                                                                                                                                                   <  2:377064‑M2/32‑4 (MQ=255)
        catcaaccacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTa                                                                                                                                                                                                                               <  1:392292‑M2/36‑1 (MQ=255)
             accacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGAt                                                                                                                                                                                                                          <  1:155773‑M2/41‑1 (MQ=255)
                 cgacgactgtctggaaatcgccgtgttgaaangtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCt                                                                                                                                                                                                                      <  2:272305‑M2/45‑1 (MQ=255)
                 cgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCt                                                                                                                                                                                                                      <  1:409819‑M2/45‑1 (MQ=255)
                  gacgattgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgtcgttatcgcccagcgcggcgtgcggcactggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggctactagcacaccagatgcctcgctaaataccctatccgaGGTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTc                                                                                                                                                                                                                     <  2:341087‑M2/46‑1 (MQ=255)
                  gacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggacgattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTCAAAAAACGGCCGTGAATTGTAACAATTCCTc                                                                                                                                                                                                                     >  2:214033‑M2/206‑251 (MQ=255)
                  gacgactgtctggaaagcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTc                                                                                                                                                                                                                     <  2:245582‑M2/46‑1 (MQ=255)
                    cgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAg                                                                                                                                                                                                                   <  2:326409‑M2/48‑1 (MQ=255)
                            tggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcatttggccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGc                                                                                                                                                                                                            <  2:407542‑M2/55‑1 (MQ=255)
                            tggaaatcgccgtgttgaaaggtgacatgggggacgtgcagcattttgccgatgacgtaatcgcccagcgcggcgggcggcacgggcatttgcagtgcttgccgacggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctcactaaataccctctccgaGTTGCGCTCCTCTTAAAAAAACGATCGAGAATTGTAAAGAttcatcaggtgacga                                                                                                                                                                                                            >  2:551276‑M2/196‑243 (MQ=255)
                             ggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGc                                                                                                                                                                                                            <  1:184840‑M2/55‑1 (MQ=255)
                                    gccgtgttgaaaggtgacatgggtgacgtgcagcatgttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATc                                                                                                                                                                                                   <  2:239992‑M2/64‑1 (MQ=255)
                                    gccgtgttgaaaggtgaaatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATc                                                                                                                                                                                                   >  2:246583‑M2/188‑251 (MQ=255)
                                     ccgtgttgaaaggggacatgggagacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATACCTCAGTAGTAGCACCAGCATc                                                                                                                                                                                                   >  2:118382‑M2/187‑250 (MQ=255)
                                         ggtgaaaggtgacaggggtgacgtgcagcattttgccgaagacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagatggagtctattagggggtcctttggggaaagatcaacaaagggccactagcacacctgattcctctctaaatcccctatccgaGTTTCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGAGTCCTCAGTGGACGCACCAGCATCAatat                                                                                                                                                                                              <  2:30967‑M2/69‑1 (MQ=255)
                                                ggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTaaa                                                                                                                                                                                       <  1:511487‑M2/76‑1 (MQ=255)
                                                ggggacttgggtgacgtgcagcattttgccgatgacgttagcgcccagctcggcgggcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTaaa                                                                                                                                                                                       <  2:199879‑M2/76‑1 (MQ=255)
                                                 gtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTaaaa                                                                                                                                                                                      <  1:160145‑M2/77‑1 (MQ=255)
                                                   gacatgggtgacgtgcagcattttgccgatgacgctatcgcccagcgcggcgtgcggcacgggcctttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCCGTGGTCGCACCAGCATCAATTTTACTAAAAgg                                                                                                                                                                                    >  2:93045‑M2/173‑251 (MQ=255)
                                                   gacatgggggacgggcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGAATCAATATTACTAAAAgg                                                                                                                                                                                    <  2:477713‑M2/79‑1 (MQ=255)
                                                    acatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGa                                                                                                                                                                                   <  1:470088‑M2/80‑1 (MQ=255)
                                                      atgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAg                                                                                                                                                                                 <  1:381721‑M2/82‑1 (MQ=255)
                                                        gggtgacgtgcagcattttgccgatgacgttatcgaccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGt                                                                                                                                                                                >  1:159004‑M2/168‑250 (MQ=255)
                                                          gtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTAtt                                                                                                                                                                             >  1:507917‑M2/166‑251 (MQ=255)
                                                            gacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTg                                                                                                                                                                           <  1:288377‑M2/88‑1 (MQ=255)
                                                                      attttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATtat                                                                                                                                                                 >  1:363604‑M2/154‑251 (MQ=255)
                                                                             cgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGAt                                                                                                                                                          <  1:295324‑M2/105‑1 (MQ=255)
                                                                              gatgacgttatcgcccaccgcggcgtgaggcacgggcatttgcagtgcttgccgaaggcagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGAt                                                                                                                                                          >  1:450890‑M2/146‑250 (MQ=255)
                                                                               atgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggcagattgagtctataattgggtcctttggtgaacgatcaaaaaagggccactcgcacacctgattcctctctcaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATCAGATa                                                                                                                                                         >  1:322876‑M2/145‑250 (MQ=255)
                                                                                 gacgttatcgcccagcgcggcgtgcggcacgngcatttgcagtgcttgcctaaggaagatttatgctattatggagtccttgggtgaacgatctacaaaggggcactatctctcctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCATTGGTCTCTCCAGCATCAATATTACTAAAAGGAAGTATTTGCCTAAATTATATAAGATAatt                                                                                                                                                      <  2:273592‑M2/109‑1 (MQ=255)
                                                                                     ttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTAttt                                                                                                                                                  >  1:128742‑M2/139‑251 (MQ=255)
                                                                                         cgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTAtttttt                                                                                                                                               <  1:122774‑M2/116‑1 (MQ=255)
                                                                                         cgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTAtttttt                                                                                                                                               <  1:509096‑M2/116‑1 (MQ=255)
                                                                                          gcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGa                                                                                                                                             <  1:39333‑M2/118‑1 (MQ=255)
                                                                                            ccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGt                                                                                                                                           <  1:539046‑M2/120‑1 (MQ=255)
                                                                                                     cgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTGTTTGAGTGAAATCCAt                                                                                                                                  <  1:411344‑M2/129‑1 (MQ=255)
                                                                                                      gtgcggcactggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacccctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAAGTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATa                                                                                                                                 <  2:127092‑M2/130‑1 (MQ=255)
                                                                                                                              ttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCt                                                                                                         >  1:79849‑M2/98‑251 (MQ=255)
                                                                                                                                   gaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTTCTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATAtt                                                                                                    <  1:487288‑M2/159‑1 (MQ=255)
                                                                                                                                      ggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCa                                                                                                 <  1:198185‑M2/162‑1 (MQ=255)
                                                                                                                                                    attattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCccgcgccg                                                                                   >  1:567182‑M2/76‑251 (MQ=255)
                                                                                                                                                        ttgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccaaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCGTCAATATTACTAAAAGGa                                                                                                                                                                                   <  2:493825‑M2/80‑1 (MQ=255)
                                                                                                                                                        ttgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccaaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGa                                                                                                                                                                                   >  1:493825‑M2/72‑151 (MQ=255)
                                                                                                                                                                      gaacgatcaacaaagggccactagcacacctgattcctctctaaataccctaaccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATg                                                                 <  1:214033‑M2/194‑1 (MQ=255)
                                                                                                                                                                           atcaacaaagggccactagcacacctgattcnccgctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAATGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGcgccg                                                            <  2:290357‑M2/199‑1 (MQ=255)
                                                                                                                                                                           atcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGcgccg                                                            >  1:254161‑M2/53‑251 (MQ=255)
                                                                                                                                                                              aacaacgggccactagctcatctgattcctctctaaatcccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATTTCAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAACAAAGTCTATATTTCACTTTGCCTGCGCCGCGAAAGTCACTGATAATGCGCCGCGt                                                         >  2:42967‑M2/50‑251 (MQ=255)
                                                                                                                                                                                 aaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTc                                                       >  1:432426‑M2/47‑250 (MQ=255)
                                                                                                                                                                                            agcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCaaa                                            <  1:573850‑M2/215‑1 (MQ=255)
                                                                                                                                                                                              cacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATg                                         <  2:377617‑M2/218‑1 (MQ=255)
                                                                                                                                                                                               ccacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATgg                                        <  1:13212‑M2/219‑1 (MQ=255)
                                                                                                                                                                                                              ctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGc                         >  2:541005‑M2/18‑251 (MQ=255)
                                                                                                                                                                                                                           cagagttaagcTCATCTTAAAAAAACGGACGAGAATTGTCCAGATTCGTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGAATTTGCCTGAATTATATAGGATAATTATCTGTTGAGTTAAATAAATCCTGGGTGCAAATAAAATAAAGCTTATCTTTCACTTTCCCCGCGCCGAGACAGTACCTGATAATGCGCCACGTTCATGTCCTCACAATGGAGTAACGACCTATGCTACActtttttg             >  2:436639/10‑250 (MQ=255)
                                                                                                                                                                                                                            cgaGTTGCGCTACTCTTAAAAAAACGGACGTGAATTGTACAGATTCCTCAGTGGTCGCCCCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATACATACCGGGTGCAAATCAAAAAAAGTCTATTTTTCACTTTGCCCGCTCCGCGAAAGTCACTGATAATGCGCCCCGTTCATGTCCTCACAANTGCGGAACGCCCTATGCTACCTTTGTTTGCt           >  2:279237‑M2/4‑251 (MQ=255)
                                                                                                                                                                                                                                   cgcTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCTTATGCTACATTTGTTTGCTGGCCTgg    <  1:511199/251‑1 (MQ=255)
                                                                                                                                                                                                                                     cgcctcTTAAAAAAACTGCCGTGAATTGTAAAGATTCCTCAGTGTTCACACCAGCATTAATATTACTAAAAGGAAATATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCAGACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCTTAACGGCCTATGCGACAATTGTTTGCTGGACTgtat  <  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: