Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,637,426 Δ1 bp intergenic (‑16/‑216) 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,4260G.98.4% 289.7 / ‑4.3 64intergenic (‑16/‑216)yhiN/pitAputative oxidoreductase YhiN/metal phosphate:H(+) symporter PitA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (21/42);  minor base C (1/0);  total (22/42)
Fisher's exact test for biased strand distribution p-value = 3.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGTGCGAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAGGTTGGTAAAGTTGCAGCGCCCACCGCCAGACATAAGGATTTTGCGCCCTGGTTTTTTACCATTATCGATCAGCAGAACCCGGCGTCCTGCCTGACCTGCCAGCGCAGAACAGAACATACCCGCCGCACCAGCGCCTATAATAATGGCATCAAACCTTTCCACGTTGCGCTCCTCTTAGAAAAAACGGGCGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGAT  >  NC_000913/3637182‑3637668
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
ccgctgcatgtgcacatcaaccacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctcgaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAAc                                                                                                                                                                                                                                             <  2:458705‑M2/22‑1 (MQ=255)
acgctgcatgtgcacatcaaccacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAAc                                                                                                                                                                                                                                             <  1:361352‑M2/22‑1 (MQ=255)
  gctgcaagtgctcatcaaccacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattgtgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtcgattattgggtcctttggtgaacgatcaacaaatggccactagcactcctgattcctctctaaataccctatccgaGTTGCGCTCCTCATA‑AAAAAACgg                                                                                                                                                                                                                                           <  2:564181‑M2/24‑1 (MQ=255)
     gcatgtgcacatcaaccacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtggacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAaacggacg                                                                                                                                                                                                                                        >  1:377617‑M2/225‑248 (MQ=255)
     gcatgtgcacatcaaccacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAaacggacg                                                                                                                                                                                                                                        >  1:61689‑M2/225‑248 (MQ=255)
      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)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
AACGTGCGAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAGGTTGGTAAAGTTGCAGCGCCCACCGCCAGACATAAGGATTTTGCGCCCTGGTTTTTTACCATTATCGATCAGCAGAACCCGGCGTCCTGCCTGACCTGCCAGCGCAGAACAGAACATACCCGCCGCACCAGCGCCTATAATAATGGCATCAAACCTTTCCACGTTGCGCTCCTCTTAGAAAAAACGGGCGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGAT  >  NC_000913/3637182‑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: ‑