Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,423,383 A→G 8.3% intergenic (+189/+40) yhdZ → / ← rrfF putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,423,3830AG 8.3% 233.2 / 8.9 81intergenic (+189/+40)yhdZ/rrfFputative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base A (33/38);  new base G (3/4);  total (38/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.24e-01

AAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGC  >  NC_000913/3423118‑3423655
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                  
aaGCTGCACCTTATGAAGTTTTTTCGCATCTTAAATTAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGccc                                                                                                                                                                                                                                                                            <  4:265144/274‑1 (MQ=255)
           gatgaATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCTCCGGGCCGTTCGTCCCCAACGTTCAAATCCGCCCCCGGCGGTTTTTGCCTCCTCGGGAGACAGTTCACCGAAAACCCCCAAACAAAACACAAACACCAGTCCTCCgcac                                                                                                                                                                                                                                                                >  4:218951/1‑272 (MQ=255)
             tgaATTTTTTGCGCAGCCTAAATCAGAGCGGACGAGGGCATTTTTAACGCAGGTAATCCATTAATTGAAGGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTg                                                                                                                                                                                                                                                               <  1:39776/274‑1 (MQ=255)
                       gcgcATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtt                                                                                                                                                                                                                                                    <  2:77173/275‑1 (MQ=255)
                                       agcttcagAGGGTATTTTTATCGCAGGTAATCCATTAATTTAATGTTAGTTCGAAAAGCAAAAAGGCACTCCTGTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCTTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTgg                                                                                                                                                                                                                                     <  4:525220/267‑1 (MQ=255)
                                       agCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGGTTTGTCCTACTCGGGAGAGTGTTCACCGCCACACCCACGATAAAACCAAAGGCCCAGTCTTCCCACCGCGCCTTTTTTTTTTTTTGATTTCTgg                                                                                                                                                                                                                                     >  4:108286/1‑274 (MQ=21)
                                       agCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTgg                                                                                                                                                                                                                                     <  2:471148/274‑1 (MQ=255)
                                       agCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTgg                                                                                                                                                                                                                                     <  2:126362/274‑1 (MQ=255)
                                       agCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTATTCCGCCTGAGCCTTTTGTTTTTTTTGTTGTcgtg                                                                                                                                                                                                                                     >  1:151157/1‑271 (MQ=255)
                                       agCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGACGGGCGTCCTGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCCGGAGAGCGTTCACCGACAAACAACAGAAAAAACAAAAGGCCCAGTCTTTCGACCAAACCTTTTGTTTTATTTGATACCTGGc                                                                                                                                                                                                                                    >  1:7816/1‑275 (MQ=18)
                                                  gCATATTTTTCGCAGTAAATCCAGTAATTGAATGTAAATACGAACAGCAAAAAGGCCGTCCGTGGGAGTGGCTGTTCGTGTGAGTTGACAGCTGGAATTTAATGGGGGCCTTCCTGCCCGCCACGCTCCGGGCCGTTGCTTCGCGACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGGGAGTGTTCACCGACAAACAACGGATAAAACAAAAGGCGCAGTCTTCCGACTGAGCCTTGTGTTTTATTTGATGTCTGGCAGTTCCCTAcc                                                                                                                                                                                                                         <  4:163431/275‑2 (MQ=14)
                                                         tATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACt                                                                                                                                                                                                                                                                <  2:107979/229‑1 (MQ=255)
                                                                gggaaTGCATTAAGTTAGTGGTAGGTCGAAAAGGAAAAAGGCTATCCTTTAGTATGGCTGTTCGCTTGATTTGATGGCTGGCAGTTTATGGGTGGCGTCCTGACAGCCACCCGCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACAGGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGGCTGCCGACTGGGCGTTTTGTTTGATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGga                                                                                                                                                                                                           <  4:97154/270‑1 (MQ=14)
                                                                         ttAATTGAATGATCGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTTGTTCGCAACGTTCAAATACGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCAACGACACACAATAGATACAACAAAACGCCCATTCTTTCGACTGAGCCTATTGGTTCATTTTATGTCTGGCACTGCACTACTCTCTCCTGAGGACAcccccacc                                                                                                                                                                                                 >  4:408045/1‑270 (MQ=17)
                                                                                 aTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCCGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCg                                                                                                                                                                                         <  2:188806/274‑1 (MQ=255)
                                                                                 aTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTCCCGCCACCCCTCGGGCCGTTACTTCCCAACGCTACAATCCCCCCCCGGCGGATTTGTCCTACTCGGGCGTGTGTTCGCCGACACACCACCGCTACAACCAAAGGCACAGTTTTCCGACTGAGCCTTTTGTTTTATTTTAACTACGGCAGTTACCTACTCCCGCATGGAGGGACCCCACACtccagacc                                                                                                                                                                                         >  4:63378/1‑267 (MQ=255)
                                                                                  tGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGAAAAAAAAAAAAGACCAGTCTTCCGACTGAGCCTTTGTTTTTTTTTGATGTCTGTCAAGTCCCTACTCCCGCGTTGGGAGAGCCCACACTAcccaccgc                                                                                                                                                                                       >  4:479560/1‑269 (MQ=17)
                                                                                          cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTTCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAAAAAACGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTACCTACTATCGCATGGGGAAACCCCACACTGCCCTCCGCGCTACg                                                                                                                                                                                 >  1:159385/1‑273 (MQ=255)
                                                                                                                ttaggttGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGTGGGGCGTCCTGCCCGCCGCCCTCCGGGCCGTTGCTTCGCAACGTTCAAAGCCGCTCCCGGCGGATTTGTACTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGc                                                                                                                                                           <  1:7741/267‑1 (MQ=255)
                                                                                                                       ggCCTTTCGCATGATTGTATGTCTGGCAGGTGATGGGGGGCGTCCTGCCAGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCGACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGt                                                                                                                                                    <  2:408046/273‑1 (MQ=255)
                                                                                                                        gCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTccc                                                                                                                                                                                                                             >  2:200573/1‑201 (MQ=255)
                                                                                                                          cTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGt                                                                                                                                               <  2:54884/275‑1 (MQ=255)
                                                                                                                            ttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTAct                                                                                                                                                                                                                         <  4:200572/201‑1 (MQ=255)
                                                                                                                              cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACCAACAACAGATAAAACAAAAGGCCCAGTCTTCCGCCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCCTGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTggg                                                                                                                                            >  2:29207/1‑274 (MQ=255)
                                                                                                                                              tGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGGTTTGTCCTACTCGGGGGAGTGTTCACCAACAAACAACCGAAAAAACAAAATGCCCAGTCTTCCCCCTGAGCCTTTTGTTTTTTTTGATGTCTGGCAATTCCCTACTCTCCCATTGGGGAACCCCCACCCACCCTCCGCGCTACGGCCTTTCACCTTTGCGTTCGGCATGGTGTCACGTTGGacccccccccctct                                                                                                                              >  4:418599/1‑264 (MQ=11)
                                                                                                                                              tGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGTCGTTTCACTTCTGAGTTCTGCATGGGGTCCGGTGGGACCACCGCGCTCCTgcc                                                                                                                           >  1:10440/1‑275 (MQ=255)
                                                                                                                                              tGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCCCTTCTGAGTTCGGCATGGGGTCAGGTGGGACCCCCGCGCTACTgcc                                                                                                                           >  1:65649/1‑275 (MQ=37)
                                                                                                                                              tGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTgcc                                                                                                                           >  2:423034/1‑275 (MQ=37)
                                                                                                                                                gCAGTTGATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGGTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTgccgc                                                                                                                         <  1:100479/275‑1 (MQ=37)
                                                                                                                                                  agtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCttt                                                                                                                                                                                                                                                        >  7:47975‑M1/113‑148 (MQ=255)
                                                                                                                                                          gcgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAAtt                                                                                                               <  2:479561/275‑1 (MQ=38)
                                                                                                                                                           cgggctCCCCGCCCGCCCCCCGCCGGGCNNNTGCTTGGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACCCGGGCGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAAtt                                                                                                               <  3:151157/274‑1 (MQ=255)
                                                                                                                                                                cccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacttcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGcc                                                                                                                                                                                                                                                           >  8:1038098‑M1/99‑131 (MQ=255)
                                                                                                                                                                cccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaCTAAAACAAAAGGCCCAGTCTTCCGACTGAGcc                                                                                                                                                                                                                                                           <  7:1038098‑M1/33‑1 (MQ=255)
                                                                                                                                                                     gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACATTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCCCGCATGGGGAGACCCCACCCTACCATCGGCCCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCCGCGGCCTCCGGCCGCCCGGCAACTTCTTTTTTatt                                                                                                     >  4:330466/1‑272 (MQ=11)
                                                                                                                                                                        cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTc                                                                                                                                                                                                                                                                        >  7:661421/1‑110 (MQ=40)
                                                                                                                                                                        cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTg                                                                                                  >  1:71758/1‑275 (MQ=255)
                                                                                                                                                                        cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCAGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTc                                                                                                                                                                                                                                                                        <  8:661421/110‑1 (MQ=255)
                                                                                                                                                                            catgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGa                                                                                                                                                                                                                                                              >  5:491071‑M1/87‑116 (MQ=255)
                                                                                                                                                                            catgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGa                                                                                                                                                                                                                                                              <  6:491071‑M1/30‑1 (MQ=255)
                                                                                                                                                                               cTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAg                                                                                                                                                                                                                                                             >  7:71976/1‑114 (MQ=40)
                                                                                                                                                                               cTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAg                                                                                                                                                                                                                                                             <  8:71976/114‑1 (MQ=40)
                                                                                                                                                                                   gggCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGa                                                                                                                                                                                                                                                                  <  8:713140/105‑1 (MQ=255)
                                                                                                                                                                                   gggCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGa                                                                                                                                                                                                                                                                  >  7:713140/1‑105 (MQ=40)
                                                                                                                                                                                   gggCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGTATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAAc                                                                                      <  2:108287/275‑1 (MQ=255)
                                                                                                                                                                                      ccGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTa                                                                                   <  4:423033/275‑1 (MQ=255)
                                                                                                                                                                                        tttggtccgaaaCGTTACAATCNNGCCCCNGCGTGTTGGTCCACCCCGGGAGGGTTTCAACCGACAAACAACAGATAAAACCAAAGGCCCTGTCTTCCGACTGGGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATTGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTaaa                                                                                 <  3:10440/265‑1 (MQ=17)
                                                                                                                                                                                         gaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTg                                                                                                                                                                                                                                                               >  7:881736‑M1/74‑102 (MQ=255)
                                                                                                                                                                                         gaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTg                                                                                                                                                                                                                                                               <  8:881736‑M1/29‑1 (MQ=255)
                                                                                                                                                                                           gttgtggaaCGTTCAAATCCGATCCCGGCGGCTTTGTCCTTGTCGGGGGAGTGTTCACCGACAAACAACAGATAAAACGAAAGACGCGGTCTTCTGTGTGAGCCGTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGAATGGGGAGAACCAACACTACCATCGGCGTTACGGCGGTTCACTTCTGAGTTCGGCGTGGGGTCAGGTGGGACCACCGTGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAg                                                                              <  4:1242/268‑1 (MQ=255)
                                                                                                                                                                                              tCGCAACGTTACAATCCGCTCCCCGCGGATTTGTCCTACTCGGGCGAGGGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTTGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCGGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAgtgg                                                                           <  2:279283/275‑1 (MQ=255)
                                                                                                                                                                                                  cgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCt                                                                                                                                                                                                                                                          <  6:358880‑M1/34‑1 (MQ=255)
                                                                                                                                                                                                  cgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCt                                                                                                                                                                                                                                                          >  5:358880‑M1/65‑98 (MQ=255)
                                                                                                                                                                                                    cGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTg                                                                      <  2:310494/274‑1 (MQ=255)
                                                                                                                                                                                                      ttCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTg                                                                      >  2:142275/1‑272 (MQ=255)
                                                                                                                                                                                                       tCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGccc                                                                                                                                                                                                                                                                            <  7:629843/75‑1 (MQ=255)
                                                                                                                                                                                                       tCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGccc                                                                                                                                                                                                                                                                            >  8:629843/1‑75 (MQ=255)
                                                                                                                                                                                                       tCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTTGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAACAGTGGTGCTGATAc                                                                  >  2:482259/1‑275 (MQ=255)
                                                                                                                                                                                                        cAAATCCGCTCCCGGCGGATTTGTCCTACGCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATAcc                                                                 <  2:242948/275‑1 (MQ=255)
                                                                                                                                                                                                          aaTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt                                                                                                                                                                                                                                                          >  5:691611/1‑90 (MQ=255)
                                                                                                                                                                                                          aaTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt                                                                                                                                                                                                                                                          >  5:791787/1‑90 (MQ=255)
                                                                                                                                                                                                          aaTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt                                                                                                                                                                                                                                                          <  6:691611/90‑1 (MQ=255)
                                                                                                                                                                                                          aaTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt                                                                                                                                                                                                                                                          <  6:791787/90‑1 (MQ=255)
                                                                                                                                                                                                                  cccGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTg                                                                                                                                                                                                                                                      >  8:149945/1‑86 (MQ=255)
                                                                                                                                                                                                                  cccGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTg                                                                                                                                                                                                                                                      <  7:149945/86‑1 (MQ=255)
                                                                                                                                                                                                                                cctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGAGcc                                                                                                                                                                                                                                                           >  7:926911‑M1/35‑67 (MQ=255)
                                                                                                                                                                                                                                cctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGAGcc                                                                                                                                                                                                                                                           <  8:926911‑M1/33‑1 (MQ=255)
                                                                                                                                                                                                                                 cTACCCGGGTGAGTGTTCACCGCCAAACANNAGATAAAACAAAAGGCCCTGTCTTCCGACGNGGCCTTTTGTTTTATTTGCTGTCTGGCAGTTCCCTACTCTCGCATGGGGGAGACCCCACACTACCATCGGCGCTNCGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCAcc                                         <  3:71758/275‑1 (MQ=21)
                                                                                                                                                                                                                                    cTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCt                                       >  2:152894/1‑273 (MQ=255)
                                                                                                                                                                                                                                      cGGGAGGGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTAc                                    <  2:538484/274‑1 (MQ=255)
                                                                                                                                                                                                                                        ggAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGa                                                                                                                                                                                                                                                                  >  7:474855/1‑52 (MQ=255)
                                                                                                                                                                                                                                        ggAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGa                                                                                                                                                                                                                                                                  <  8:474855/52‑1 (MQ=255)
                                                                                                                                                                                                                                        ggAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCa                                  <  2:345908/274‑1 (MQ=255)
                                                                                                                                                                                                                                         gagaGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAg                                                                                                                                                                                                                                  >  4:570919/1‑83 (MQ=255)
                                                                                                                                                                                                                                             gtgtTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTcc                                                                                                                                                                                                                              <  2:570920/83‑1 (MQ=255)
                                                                                                                                                                                                                                              tgtTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTgc                           <  2:314398/275‑1 (MQ=255)
                                                                                                                                                                                                                                                       acaaacaaCAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGACTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCa                  >  4:451729/1‑275 (MQ=255)
                                                                                                                                                                                                                                                             aacaGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAg            <  4:142274/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                       aCAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGTCGCTACGGCGTTTCACTTCTGAGTTAGGCATGGGGTCAGGTCAGACCAGCGCGCCACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTCATCTTAAAAGTCGTGCTGATACACAGTGTCGCACTGGGGAACTACCCCTTACCAAGGGTGGCCCCTCCCAACCACGCCatattcac  >  4:33032/1‑271 (MQ=18)
                                                                                                                                                                                                                                                                       aCAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAAc                >  4:125905/1‑261 (MQ=255)
                                                                                                                                                                                                                                                                       aCAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCCCACCCTTACCAAGGGTCCGCCCTACCAACTGAGCCATATCAGc  >  4:507444/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                  
AAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGC  >  NC_000913/3423118‑3423655

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