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.1% 191.3 / ‑6.1 53intergenic (‑26/‑206)yhiN/pitAputative oxidoreductase YhiN/metal phosphate:H(+) symporter PitA
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (17/34);  total (18/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

CTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAGGTTGGTAAAGTTGCAGCGCCCACCGCCAGACATAAGGATTTTGCGCCCTGGTTTTTTACCATTATCGATCAGCAGAACCCGGCGTCCTGCCTGACCTGCCAGCGCAGAACAGAACATACCCGCCGCACCAGCGCCTATAATAATGGCATCAAACCTTTCCACGTTGCGCTCCTCTTAGAAAAAACGGGCGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGG  >  NC_000913/3637198‑3637680
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    
tcaaccacgacgactgtctggaaatcgccgcgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggcagattgagtctattattgggtcctttggtgaacgatcaacaaagggcccctagcacacctgattcctctctaactaccctatccgaGTTGAGCTCCTCTTAAAAAAACGGACGTTAATTGTaa                                                                                                                                                                                                                                          >  2:351503‑M2/214‑250 (MQ=255)
     cacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcagtttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgatgcctctctaaataccctagccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTc                                                                                                                                                                                                                                    <  2:338330‑M2/43‑1 (MQ=255)
          cgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAg                                                                                                                                                                                                                               <  1:185667‑M2/48‑1 (MQ=255)
                   ggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCa                                                                                                                                                                                                                       <  1:403178‑M2/56‑1 (MQ=255)
                          gccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATc                                                                                                                                                                                                               <  1:143649‑M2/64‑1 (MQ=255)
                           cagtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATc                                                                                                                                                                                                               >  2:463160‑M2/187‑250 (MQ=255)
                                 tgaaaggtgacgtgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcagcaaagggccactagcacacctgattcctctctaaataccctagccgaGTTGCGCTCCTCTGAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCGTCAATAtt                                                                                                                                                                                                         <  2:406620‑M2/70‑1 (MQ=255)
                                 tgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagtgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATAtt                                                                                                                                                                                                         >  2:293437‑M2/181‑250 (MQ=255)
                                  gaaaggtgacatgggtgacgtgcagcattgtgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTAc                                                                                                                                                                                                       <  2:467250‑M2/72‑1 (MQ=255)
                                   aaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACt                                                                                                                                                                                                      <  1:349817‑M2/73‑1 (MQ=255)
                                   aaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaacaccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACt                                                                                                                                                                                                      <  1:252219‑M2/73‑1 (MQ=255)
                                    aaggtgacatgggtgacgtgcagcattttgccgatgccgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcccttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTa                                                                                                                                                                                                     <  1:151930‑M2/74‑1 (MQ=255)
                                     aggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagagcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTa                                                                                                                                                                                                     >  1:124337‑M2/177‑250 (MQ=255)
                                      ggtggcatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggctactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGGCGTGATTGGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTaaa                                                                                                                                                                                                   <  1:427250‑M2/76‑1 (MQ=255)
                                            atgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGa                                                                                                                                                                                               >  1:247190‑M2/170‑249 (MQ=255)
                                              gggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGt                                                                                                                                                                                            <  1:306012‑M2/83‑1 (MQ=255)
                                                  gacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTAttt                                                                                                                                                                                        <  1:490916‑M2/87‑1 (MQ=255)
                                                       gcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTg                                                                                                                                                                                   >  2:549758‑M2/159‑250 (MQ=255)
                                                       gcagcatgttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattagtgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTTCCTg                                                                                                                                                                                   <  2:227319‑M2/92‑1 (MQ=255)
                                                          gcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAAtt                                                                                                                                                                               <  1:70702‑M2/96‑1 (MQ=255)
                                                           cattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATta                                                                                                                                                                              <  1:9128‑M2/97‑1 (MQ=255)
                                                              tttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATtata                                                                                                                                                                            <  1:259175‑M2/99‑1 (MQ=255)
                                                              tttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATtata                                                                                                                                                                            <  1:214476‑M2/99‑1 (MQ=255)
                                                                 gccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCATCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATaa                                                                                                                                                                         <  1:138231‑M2/102‑1 (MQ=255)
                                                                  ccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATtatata                                                                                                                                                                          >  1:424055‑M2/148‑248 (MQ=255)
                                                                   cgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTAttt                                                                                                                                                                                        <  2:88968‑M2/87‑1 (MQ=255)
                                                                   cgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTAttt                                                                                                                                                                                        >  1:88968‑M2/147‑233 (MQ=255)
                                                                       gacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAatt                                                                                                                                                                  >  2:85653‑M2/143‑251 (MQ=255)
                                                                          gttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAattatt                                                                                                                                                               >  1:130430‑M2/140‑251 (MQ=255)
                                                                                ccccagcgcggcgtgcggcacggtcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGa                                                                                                                                                         <  2:532576‑M2/118‑1 (MQ=255)
                                                                                   cagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTg                                                                                                                                                      >  2:395553‑M2/131‑251 (MQ=255)
                                                                                              gcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATAc                                                                                                                                            <  2:504474‑M2/131‑1 (MQ=255)
                                                                                                  cacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAgggg                                                                                                                                       >  1:75345‑M2/116‑251 (MQ=255)
                                                                                                       gaatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCCCCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCaa                                                                                                                                   >  2:374040‑M2/111‑250 (MQ=255)
                                                                                                            tgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCaa                                                                                                                              <  1:538089‑M2/145‑1 (MQ=255)
                                                                                                                  gcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGCCGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATATTCATGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATCCAGGGGGCAAATCAAAAAAAGt                                                                                                                       >  2:23673‑M2/100‑251 (MQ=255)
                                                                                                                        cgaaggaagattgatgctagtattggggcctttggtgaacgttcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATGCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTAGTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTatat                                                                                                                 <  2:343191‑M2/158‑1 (MQ=255)
                                                                                                                                         tattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAATGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCACATCAAAAAAAGTCTATATTTCACTTTGCccgcgcc                                                                                                >  1:46648‑M2/77‑251 (MQ=255)
                                                                                                                                         tattattggcgcctttggtgaaagatgaacaaagggccactagcacacctgattcctctctaaataccctatccgaATTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGTTCGCACCAGCATCAATATTACTAAAAGGAAGAATTTGCCTGCATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCccgcgcc                                                                                                >  2:471797‑M2/77‑251 (MQ=255)
                                                                                                                                                    cctttggtgaactatcaacaaagggccactagcacacctgattcctctctaaatacccgatccgaGTTGCTCTCCTCTTAAAAAAACGGACGTTAATTGTAAATATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGTCGCAAAAGTCa                                                                                      <  2:200688‑M2/185‑1 (MQ=255)
                                                                                                                                                        tggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGAt                                                                                 >  2:280343‑M2/62‑251 (MQ=255)
                                                                                                                                                                    aacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGt                                                                     <  2:14887‑M2/202‑1 (MQ=255)
                                                                                                                                                                     acaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGtt                                                                    <  1:85653‑M2/203‑1 (MQ=255)
                                                                                                                                                                                  agcacacctgattactctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGTACGTGACTTGTAAAGATGCCTCAGTGGTCGAACCAGCATCAATATTACTAAAAGGAAGTAGTTGCCTGAATTATATAAGATAATTATTGTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCaaaa                                                       <  2:486324‑M2/216‑1 (MQ=255)
                                                                                                                                                                                           gattcctctctaaatgccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATGCCTCAGTGGTCGCAACAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTa                                                <  2:304248‑M2/223‑1 (MQ=255)
                                                                                                                                                                                           gattcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTa                                                <  2:187904‑M2/223‑1 (MQ=255)
                                                                                                                                                                                              tcctctctaaataccctatccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGt                                            <  1:183427‑M2/227‑1 (MQ=255)
                                                                                                                                                                                                               atccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCTCTTTGCCCGCGACGCGAAAGTCTCTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCTTAACGTCCTATGCTACAtttgttt                          <  1:467550‑M2/245‑1 (MQ=255)
                                                                                                                                                                                                                 ccgaGTTGCGCTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACAtttgtttg                         <  1:549758‑M2/246‑1 (MQ=255)
                                                                                                                                                                                                                        gcgcTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCt                  <  1:207170/250‑1 (MQ=255)
                                                                                                                                                                                                                          gcTCCTCTTAAAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGCTAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTgg                >  2:541009/1‑250 (MQ=255)
                                                                                                                                                                                                                                   taaaaaaaCGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAACATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGc          <  1:366622/246‑1 (MQ=255)
                                                                                                                                                                                                                                        aaaCGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCggg  <  1:293437/251‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    
CTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAGGTTGGTAAAGTTGCAGCGCCCACCGCCAGACATAAGGATTTTGCGCCCTGGTTTTTTACCATTATCGATCAGCAGAACCCGGCGTCCTGCCTGACCTGCCAGCGCAGAACAGAACATACCCGCCGCACCAGCGCCTATAATAATGGCATCAAACCTTTCCACGTTGCGCTCCTCTTAGAAAAAACGGGCGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGG  >  NC_000913/3637198‑3637680

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