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

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

GCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAGGTTGGTAAAGTTGCAGCGCCCACCGCCAGACATAAGGATTTTGCGCCCTGGTTTTTTACCATTATCGATCAGCAGAACCCGGCGTCCTGCCTGACCTGCCAGCGCAGAACAGAACATACCCGCCGCACCAGCGCCTATAATAATGGCATCAAACCTTTCCACGTTGCGCTCCTCTTAGAAAAAACGGGCGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCT  >  NC_000913/3637193‑3637682
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
gcacatcaaccacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAatg                                                                                                                                                                                                                                                <  2:90554‑M2/33‑2 (MQ=255)
        accacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTaa                                                                                                                                                                                                                                            >  1:231505‑M2/211‑247 (MQ=255)
             gacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTc                                                                                                                                                                                                                                   >  1:330180‑M2/206‑251 (MQ=255)
                gactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGt                                                                                                                                                                                                                                >  2:246636‑M2/203‑251 (MQ=255)
                gactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGt                                                                                                                                                                                                                                >  1:310389‑M2/203‑251 (MQ=255)
                   tgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagaattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGGCGTGATATGTAAAGATTCCTCAGTGGt                                                                                                                                                                                                                             >  2:27717‑M2/200‑251 (MQ=255)
                        ggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCAc                                                                                                                                                                                                                        <  1:246636‑M2/57‑1 (MQ=255)
                              cgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggttaacgatcaacaaagggccactagcacaccttattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACcatcat                                                                                                                                                                                                                  <  1:283721‑M2/63‑1 (MQ=255)
                                 cgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCaa                                                                                                                                                                                                               <  2:129249‑M2/66‑1 (MQ=255)
                                  gtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaattccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAat                                                                                                                                                                                                              <  2:559916‑M2/67‑1 (MQ=255)
                                  gtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAat                                                                                                                                                                                                              >  2:377015‑M2/185‑251 (MQ=255)
                                  gtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAat                                                                                                                                                                                                              <  1:553044‑M2/67‑1 (MQ=255)
                                       gaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTAc                                                                                                                                                                                                         <  1:720945‑M2/72‑1 (MQ=255)
                                         aaggttacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctagccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATACCTCAGTGGGCGCACCAGCATCAATATTACTa                                                                                                                                                                                                       <  2:484681‑M2/74‑1 (MQ=255)
                                         aaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcactcctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGAGTCCTCAGTGGTCGCACCAGCATCAATATTACTa                                                                                                                                                                                                       <  2:187678‑M2/74‑1 (MQ=255)
                                          aggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTa                                                                                                                                                                                                       >  2:111318‑M2/177‑250 (MQ=255)
                                          aggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTa                                                                                                                                                                                                       >  1:428164‑M2/177‑250 (MQ=255)
                                          aggggacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTaa                                                                                                                                                                                                      <  2:411768‑M2/75‑1 (MQ=255)
                                           ggtgacatgggtggcgtgccgcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattatcgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACTGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTaaa                                                                                                                                                                                                     <  1:231316‑M2/76‑1 (MQ=255)
                                            gtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTaaa                                                                                                                                                                                                     >  1:150806‑M2/175‑250 (MQ=255)
                                              gacatgggcgacgtgcagcattttgacgacgacgttctcgcccagcgcggcgtgcgacacgcgcttttgcagtgcttgccgaaggaagattgtgtctattattgggtcctttggtgaccgatcaacaacgggccactagcacacctgattcctctctaaataccctatccgaGCTGCGCTCCTCTTAACAAAACGGACGTGAATTGTAAAGATTCCTCAGTGCTCGCACCAGCATCACTATTACTAAAAgg                                                                                                                                                                                                  >  1:585928‑M2/173‑251 (MQ=255)
                                                       gacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTAttt                                                                                                                                                                                          <  1:662539‑M2/87‑1 (MQ=255)
                                                       gacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTg                                                                                                                                                                                         >  1:374480‑M2/164‑251 (MQ=255)
                                                        acgtgcagcatgttgccgatgccgttatcgcccatctcggcgggcggcacgggcattttcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgacccatcaacaaagggccactagcacacctgattcctctctaaataccctatcctaGTTGCGCTCCTCGTAAAAAAACGGACGTGAATTGTAAAGATTCCTCCGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGc                                                                                                                                                                                        <  2:481631‑M2/89‑1 (MQ=255)
                                                              agcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAAt                                                                                                                                                                                  <  2:561428‑M2/95‑1 (MQ=255)
                                                                 attttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATtat                                                                                                                                                                               >  1:242966‑M2/154‑251 (MQ=255)
                                                                      gccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATaa                                                                                                                                                                           <  1:11336‑M2/102‑1 (MQ=255)
                                                                      gccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGAATTTGCCTGAATtatata                                                                                                                                                                            >  1:495860‑M2/149‑249 (MQ=255)
                                                                         gatgacgttatcgccccgcgcggcgtgcggcacgggcctttgcagtgcttgccgaaggcagcttgagtctattattgggtcctttggtgaacgatcaacaacgggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAACAAACGGACGTGAAGTGTAACGATTCCTCAGTGGTCGCACCAGCATCATTCTTACTAACAGGAagtatttgcctgactttaatacgat                                                                                                                                                                        >  2:134114‑M2/146‑238 (MQ=255)
                                                                         gatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATa                                                                                                                                                                       <  1:141741‑M2/106‑1 (MQ=255)
                                                                          atgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGTACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATaa                                                                                                                                                                      <  2:310389‑M2/107‑1 (MQ=255)
                                                                               gttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAattatt                                                                                                                                                                 >  1:588547‑M2/140‑251 (MQ=255)
                                                                                    cgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTAtttttt                                                                                                                                                             <  2:141077‑M2/116‑1 (MQ=255)
                                                                                       ccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGt                                                                                                                                                         <  1:698012‑M2/120‑1 (MQ=255)
                                                                                             cggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATc                                                                                                                                                   <  1:323308‑M2/126‑1 (MQ=255)
                                                                                                       cacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAgggg                                                                                                                                         <  1:332384‑M2/136‑1 (MQ=255)
                                                                                                        acgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTTCTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAggggg                                                                                                                                        <  1:436811‑M2/137‑1 (MQ=255)
                                                                                                           ggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTTAGTGAAATCCATACAGGGGGCaa                                                                                                                                     <  2:231505‑M2/140‑1 (MQ=255)
                                                                                                              atttgcagtacttgccgaaggaagaatgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacaccggaatcctctctaaatacactaaccgaGTAGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCATCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGACTTATATAAGATAATTGTTTTTTTAGTGAAATCCATACAGGGGGAAAAt                                                                                                                                   >  2:551052‑M2/109‑250 (MQ=255)
                                                                                                                   cagtgcttgccgaaggaagattgagactattattgggtcctttggcgaacgatcaccaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGAGAATTGTAAAGATTCCTCAGTGGTCGCACCAGAATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCaaaa                                                                                                                              >  2:644564‑M2/104‑250 (MQ=255)
                                                                                                                     gtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAAGCaaaaaaa                                                                                                                           <  2:431620‑M2/150‑1 (MQ=255)
                                                                                                                           gccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAAGTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCtat                                                                                                                     <  2:476241‑M2/156‑1 (MQ=255)
                                                                                                                                aggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTc                                                                                                                <  1:647184‑M2/161‑1 (MQ=255)
                                                                                                                                 ggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTc                                                                                                                >  1:361009‑M2/90‑250 (MQ=255)
                                                                                                                                      attgagtctattattgggtcctttggtgaacgatcaacaaagggacactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAGAGTCTATATTTCACTTTg                                                                                                          <  2:614052‑M2/167‑1 (MQ=255)
                                                                                                                                           gtctattattgggtcctttggtgaacgctcaacaaagggccactagcaaacctgcttcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGACTTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTAATAAAAGGAAGTATTTGCCTGCATTATATACGCTAATTATTTTTTGAGTGAAATACATCCAGGGGACAAATCAAAAAAAGTCTATATTTCACTTTGAccgc                                                                                                     >  2:473161‑M2/80‑251 (MQ=255)
                                                                                                                                               attattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCccgcgcc                                                                                                  <  2:495860‑M2/175‑1 (MQ=255)
                                                                                                                                               attattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCccgcgcc                                                                                                  <  2:428164‑M2/175‑1 (MQ=255)
                                                                                                                                                           tttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCCATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGTGAGAAATCCATACAGGGGGCAAATCAAACAAAGTCTATATTTCACTTTGCACGCGCCGCGAAAGTCACTg                                                                                     >  2:666395‑M2/64‑251 (MQ=255)
                                                                                                                                                               gtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTT‑‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCa                                                                                                                                 <  1:440714‑M2/143‑1 (MQ=255)
                                                                                                                                                               gtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTT‑‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCa                                                                                                                                 >  2:440714‑M2/60‑202 (MQ=255)
                                                                                                                                                                        caacaaagggcctctagtccaccggattcctctctaaatcccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAAGTGTAAAGATGCCTCAGTGGTCGCACCAGCATCAATATGACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGGCGAGATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATaatgcgcagc                                                                         <  2:585928‑M2/200‑4 (MQ=255)
                                                                                                                                                                                    actagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGTCGCGTTCATGTCCTCa                                                            <  2:47525‑M2/213‑1 (MQ=255)
                                                                                                                                                                                      tagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCaa                                                           >  1:68850‑M2/37‑250 (MQ=255)
                                                                                                                                                                                           cacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGGCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGc                                                     <  1:340130‑M2/220‑1 (MQ=255)
                                                                                                                                                                                            cactgatctctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTAGCACCAGCATCAATATTACTAAAAGTAAGTAGTTGCCTGAATTATATAAGATAATTATTGTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCg                                                    <  2:555797‑M2/221‑1 (MQ=255)
                                                                                                                                                                                                gattactctctacataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTCACAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCACGCGCCGAGACAGTCACTGATAATGCGCCACGTTCCTGTCATCAAAATGGCGTAAc                                                >  2:285272‑M2/27‑251 (MQ=255)
                                                                                                                                                                                                       ctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGGAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTAt                                         >  1:249116‑M2/20‑251 (MQ=255)
                                                                                                                                                                                                           aaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATg                                        <  1:611665‑M2/233‑1 (MQ=255)
                                                                                                                                                                                                                      ccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGc                          <  1:318442‑M2/247‑1 (MQ=255)
                                                                                                                                                                                                                      ccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTAGTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACAtttgtttg                           <  1:86456‑M2/246‑1 (MQ=255)
                                                                                                                                                                                                                          gTTGCGCTCCTCTTAAAAAAAAGGACGTGAATTGTAAAGATTCCTCAGAGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTATTTGCTGGc                      >  1:360093/1‑251 (MQ=255)
                                                                                                                                                                                                                           ttGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAattat                                                                                                                                                                  >  1:679473/1‑110 (MQ=255)
                                                                                                                                                                                                                           ttGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGcc                     >  1:256082/1‑251 (MQ=255)
                                                                                                                                                                                                                           ttGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTAGATAAGATAattat                                                                                                                                                                  <  2:679473/110‑1 (MQ=255)
                                                                                                                                                                                                                             gcgcTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTTCTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCt                    <  2:378931/250‑1 (MQ=255)
                                                                                                                                                                                                                                  cctcTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGAtt               <  1:111318/250‑1 (MQ=255)
                                                                                                                                                                                                                                             aaaCGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCggg    >  1:659811/1‑251 (MQ=255)
                                                                                                                                                                                                                                               aCGGACGTGAATTGTAAAGATTACACAGTGGNCGCAACAGCATCAATATTACTCAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTGTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAATTCTATATTTCGCATTGCCCGCACCGCGAAAGTCACTGATAATGCGCGGCGTTCATGTCCTCGAAATGGCGTAACGTCCTATGATACATTTGTTTGCTGGCTTGGCTGTGCATACCGGGCt  <  2:339906/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
GCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAGGTTGGTAAAGTTGCAGCGCCCACCGCCAGACATAAGGATTTTGCGCCCTGGTTTTTTACCATTATCGATCAGCAGAACCCGGCGTCCTGCCTGACCTGCCAGCGCAGAACAGAACATACCCGCCGCACCAGCGCCTATAATAATGGCATCAAACCTTTCCACGTTGCGCTCCTCTTAGAAAAAACGGGCGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCT  >  NC_000913/3637193‑3637682

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