Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,392,617 C→G 7.2% noncoding (35/76 nt) glyY → tRNA‑Gly

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,392,6170CG 7.2% 418.3 / 18.1 122noncoding (35/76 nt)glyYtRNA‑Gly
Reads supporting (aligned to +/- strand):  ref base C (52/61);  new base G (4/5);  total (56/66)
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.49e-01

Only 100 of 125 total aligned reads displayed.

AAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCTGTAATCTTAATTTCAAAGAAACTTCGCACGGTGAATAGTA  >  NC_000913/4392345‑4392802
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                         
aaaaTTAGTTATTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTGGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGcc                                                                                                                                                                                          <  2:487743/274‑1 (MQ=255)
                     atatgtCAGTTGGTAGAACACGATCTTGCCAAGGTCGTGGTCGCGAGTTAGAGTATCGTTTCCCGGTAGAAAATTTAAAAAGTGCTGCAAAGCAAAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGa                                                                                                                                                                    <  1:60968/270‑1 (MQ=18)
                                                      gTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGTAATAGCTCAGTTGGTAGAGCACGAACTTACCAAGGTCGGGGTCGCGATTTCGAGTCCCGTTTCCCGCCCCAAACTTTCAAAGTGTGTTAAGGCACAAACCACCCAAGCGTGACTATCCCGGTTGTTTGACACCGACCTTGCCATGGCCGCGGTCGCGAGTTCGGGTCTCTTTTCCCGTCTCAACTTCTtctctccaa                                                                                                                                   >  4:417315/1‑272 (MQ=17)
                                                                 gatatagagaCTCATTTACCGAACAAAAATCTGTAACGTGCGACAAAACAGCGAGAACCCAATCGGGAAAAGCTCAGTTGGTAGAGCACGACCTTGCAAAAGCCGGGGTCGCGATGTCGAGTCTCGTTTCCCGCCCCAAAATTGGAAAATGCTGTAAGGCAAAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCGCAATAAAATATCca                                                                                                                         <  4:118453/268‑1 (MQ=14)
                                                                           ctcGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAAGAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGAcgcg                                                                                                               <  2:309196/274‑1 (MQ=255)
                                                                                                     aGTGCTGCAAAGCACAGACAACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAAt                                                                                                                                   <  2:503777/228‑1 (MQ=255)
                                                                                                           gcAAAGCACAGACCACCCAAGCGGGAATAGATCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTGGGAAGTTCGGGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACa                                                                              <  2:181645/275‑1 (MQ=255)
                                                                                                             aaaGCCCAGACCACCCAATCGGGAATAGTTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCAAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACaaa                                                                            <  2:392507/275‑1 (MQ=255)
                                                                                                                        caccCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGGCTCGTGTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTaaa                                                                  <  2:417316/274‑1 (MQ=255)
                                                                                                                             aaGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCTTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGc                                                                                                                                                                                           <  5:326141/148‑1 (MQ=255)
                                                                                                                                cGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCaa                                                                                                                                                                                        <  7:776863/148‑1 (MQ=38)
                                                                                                                                    aaTAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTATCAATAAAATCTCCACAGCGACGCGTTGAGTTATTGCTGGTTTTTGTTGTCCCTGACAAACTCTTGTAAACAGAGGTATCcac                                                     >  2:536642/1‑275 (MQ=255)
                                                                                                                                             gTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcgcg                                                                                                                                                                           <  7:153062/148‑1 (MQ=38)
                                                                                                                                                 gTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGtt                                                                                                                                                                       <  7:1096005/148‑1 (MQ=39)
                                                                                                                                                              ttGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGt                                                                                                                                                                        >  5:668003/1‑134 (MQ=40)
                                                                                                                                                              ttGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGt                                                                                                                                                                        <  6:668003/134‑1 (MQ=40)
                                                                                                                                                                   aaGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCt                                                                                                                                                     >  7:920938/1‑148 (MQ=39)
                                                                                                                                                                     ggTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCg                                                                                                                                                                     >  7:530965/1‑130 (MQ=255)
                                                                                                                                                                     ggTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCg                                                                                                                                                                     <  8:530965/130‑1 (MQ=40)
                                                                                                                                                                          gggTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGt                                                                                                                                                                                     >  8:246703/1‑109 (MQ=255)
                                                                                                                                                                          gggTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGt                                                                                                                                                                                     <  7:246703/109‑1 (MQ=255)
                                                                                                                                                                              cgcgAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgg                                                                                                                                                                                  <  5:57883/108‑1 (MQ=255)
                                                                                                                                                                              cgcgAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgg                                                                                                                                                                                  >  6:57883/1‑108 (MQ=255)
                                                                                                                                                                               gcgAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCt                                                                                                  <  2:83622/187‑1 (MQ=255)
                                                                                                                                                                               gcgAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCGCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATc                                                                                                                           <  4:485779/162‑1 (MQ=255)
                                                                                                                                                                                cgAGTTCGAGTCTCGTTTCCCGCTCCACAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTtctc                                                                                                                                        <  8:612427/148‑1 (MQ=255)
                                                                                                                                                                                cgAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGGCAAGGTCgtg                                                                                                                                                                                 >  6:804760/1‑105 (MQ=255)
                                                                                                                                                                                cgAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGGCAAGGTCgcg                                                                                                                                                                                 <  5:804760/107‑3 (MQ=255)
                                                                                                                                                                                 gAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAttcttc                                                                                                                                          <  8:773874/145‑1 (MQ=255)
                                                                                                                                                                                 gAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGGTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAttcttc                                                                                                                                          >  7:773874/1‑145 (MQ=255)
                                                                                                                                                                                    ttCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgg                                                                                                                                                                                  <  8:266130/102‑1 (MQ=255)
                                                                                                                                                                                    ttCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgg                                                                                                                                                                                  >  7:266130/1‑102 (MQ=255)
                                                                                                                                                                                       aattcccggctCCCTCTCCAAAATTGGAAAGTGGTGTAAGGACAAGACCACCCAAGCGGGAATAGTGCAGTTGGTAGAGCCCGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACATCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCTTTAATCTTAATTTCAAAGAAACTTCGCACGGTGAATAGTa  <  4:536641/265‑1 (MQ=37)
                                                                                                                                                                                        aGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCTGTAATCTTAATTTCAAAGAAACTTCGCACGGTGAATAGTa  <  2:530280/274‑1 (MQ=255)
                                                                                                                                                                                        aGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCg                                                                                                                                                                     <  7:932382/111‑1 (MQ=40)
                                                                                                                                                                                        aGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCa                                                                                                                                     >  6:817324/1‑143 (MQ=255)
                                                                                                                                                                                        aGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCa                                                                                                                                     <  5:817324/143‑1 (MQ=255)
                                                                                                                                                                                           ctcGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCa                                                                                                                                     >  7:1058258/1‑140 (MQ=255)
                                                                                                                                                                                           ctcGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCa                                                                                                                                     <  8:1058258/140‑1 (MQ=255)
                                                                                                                                                                                             cGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAttctt                                                                                                                                           >  6:1124745/1‑132 (MQ=255)
                                                                                                                                                                                             cGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAttctt                                                                                                                                           <  5:1124745/132‑1 (MQ=255)
                                                                                                                                                                                               tttCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAttct                                                                                                                                            >  6:836674/1‑129 (MQ=255)
                                                                                                                                                                                                    cGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCg                                                                                                                    >  7:775329/1‑148 (MQ=255)
                                                                                                                                                                                                          aaaaTTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGt                                                                                                                                                                  >  7:283907/1‑96 (MQ=255)
                                                                                                                                                                                                          aaaaTTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGttt                                                                                                                                                           >  7:724561/1‑103 (MQ=255)
                                                                                                                                                                                                          aaaaTTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGttt                                                                                                                                                           <  8:724561/103‑1 (MQ=255)
                                                                                                                                                                                                          aaaaTTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAttcttc                                                                                                                                          >  7:904384/1‑120 (MQ=255)
                                                                                                                                                                                                          aaaaTTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAttcttc                                                                                                                                          <  8:904384/120‑1 (MQ=255)
                                                                                                                                                                                                           aaaTTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGGCAAg                                                                                                                                                                                       >  8:71820/1‑74 (MQ=255)
                                                                                                                                                                                                           aaaTTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGGCAAg                                                                                                                                                                                       <  7:611159/74‑1 (MQ=255)
                                                                                                                                                                                                             aTTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCg                                                                                                                                                                                   <  5:526671/76‑1 (MQ=255)
                                                                                                                                                                                                             aTTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCg                                                                                                                                                                                   >  6:526671/1‑76 (MQ=255)
                                                                                                                                                                                                               ttGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCca                                                                                                                         >  5:105235/1‑132 (MQ=255)
                                                                                                                                                                                                                  aaaGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGa                                                                                                                                                                    <  5:112246/86‑1 (MQ=38)
                                                                                                                                                                                                                  aaaGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGa                                                                                                                                                                    >  6:112246/1‑86 (MQ=40)
                                                                                                                                                                                                                   aaGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTAtt                                                                                                     >  6:471152/1‑148 (MQ=255)
                                                                                                                                                                                                                    aGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTg                                                                                                    <  5:723210/148‑1 (MQ=255)
                                                                                                                                                                                                                          gTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAatat                                                                                                                            >  8:24532/1‑118 (MQ=255)
                                                                                                                                                                                                                          gTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTg                                                                                                    >  7:246652/1‑142 (MQ=255)
                                                                                                                                                                                                                          gTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTg                                                                                                    <  8:246652/142‑1 (MQ=255)
                                                                                                                                                                                                                               gCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTttgt                                                                                         <  8:920938/148‑1 (MQ=255)
                                                                                                                                                                                                                                cacaGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGc                                                                                                                     >  4:468337/1‑119 (MQ=255)
                                                                                                                                                                                                                                cacaGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGa                                                                                                              >  5:879363/1‑126 (MQ=255)
                                                                                                                                                                                                                                cacaGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGa                                                                                                              <  6:879363/126‑1 (MQ=255)
                                                                                                                                                                                                                                  caGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAttcttc                                                                                                                                          >  6:930498/1‑96 (MQ=37)
                                                                                                                                                                                                                                  caGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAttcttc                                                                                                                                          <  5:930498/96‑1 (MQ=37)
                                                                                                                                                                                                                                  caGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCt                                                                                                  >  6:576084/1‑136 (MQ=255)
                                                                                                                                                                                                                                  caGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCt                                                                                                  <  5:576084/136‑1 (MQ=255)
                                                                                                                                                                                                                                   aGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCcaca                                                    >  1:27546/1‑181 (MQ=255)
                                                                                                                                                                                                                                    gACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGAcg                                                                                                                 <  2:468338/119‑1 (MQ=255)
                                                                                                                                                                                                                                     accaccCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAttct                                                                                                                                            >  7:107254/1‑91 (MQ=34)
                                                                                                                                                                                                                                     accaccCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAttct                                                                                                                                            <  8:107254/91‑1 (MQ=34)
                                                                                                                                                                                                                                       caccCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGCCt                                                <  3:27546/181‑1 (MQ=255)
                                                                                                                                                                                                                                           cAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAttct                                                                                                                                            <  6:663908/85‑1 (MQ=35)
                                                                                                                                                                                                                                           cAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAatat                                                                                                                            >  5:1112804/1‑101 (MQ=255)
                                                                                                                                                                                                                                           cAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAatat                                                                                                                            <  6:1112804/101‑1 (MQ=255)
                                                                                                                                                                                                                                                 ggAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATcc                                                                                                                          <  5:573133/97‑1 (MQ=255)
                                                                                                                                                                                                                                                   aaTAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGttttt                                                                                           <  7:531666/126‑1 (MQ=255)
                                                                                                                                                                                                                                                   aaTAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGttttt                                                                                           >  8:531666/1‑126 (MQ=255)
                                                                                                                                                                                                                                                      aGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCtt                                                                       <  7:825458/143‑1 (MQ=255)
                                                                                                                                                                                                                                                      aGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTATCTCTGACAAACTCtt                                                                       >  8:825458/1‑143 (MQ=255)
                                                                                                                                                                                                                                                          cAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGttttt                                                                                           <  5:875033/119‑1 (MQ=255)
                                                                                                                                                                                                                                                          cAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGttttt                                                                                           >  6:200231/1‑119 (MQ=255)
                                                                                                                                                                                                                                                            gTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCGACAGCGACGCGATGCg                                                                                                          >  5:391888/1‑102 (MQ=255)
                                                                                                                                                                                                                                                            gTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCGACAGCGACGCGATGCg                                                                                                          <  6:391888/102‑1 (MQ=255)
                                                                                                                                                                                                                                                             ttGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCcaca                                                                                                                       <  5:328271/88‑1 (MQ=255)
                                                                                                                                                                                                                                                             ttGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCcaca                                                                                                                       >  6:328271/1‑88 (MQ=255)
                                                                                                                                                                                                                                                             ttGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACag                                                               <  5:867310/144‑1 (MQ=255)
                                                                                                                                                                                                                                                             ttGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACag                                                               >  6:867310/1‑144 (MQ=255)
                                                                                                                                                                                                                                                                 tAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGtctc                                                                                   <  5:202491/120‑1 (MQ=255)
                                                                                                                                                                                                                                                                 tAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTaa                                                                   >  7:1031301/1‑136 (MQ=255)
                                                                                                                                                                                                                                                                  agagCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCTGTAATCTTAATTTCAAAGAAACTTCGCACGGTGAAt      >  1:88849/1‑196 (MQ=255)
                                                                                                                                                                                                                                                                   gagCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGtttt                                                                                            >  8:1055154/1‑109 (MQ=255)
                                                                                                                                                                                                                                                                   gagCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGtttt                                                                                            <  7:1055154/109‑1 (MQ=255)
                                                                                                                                                                                                                                                                      cACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCg                                                                                                          >  8:366214/1‑92 (MQ=255)
                                                                                                                                                                                                                                                                      cACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAAc                                                                 <  8:460202/133‑1 (MQ=255)
                                                                                                                                                                                                                                                                      cACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAAc                                                                 >  7:460202/1‑133 (MQ=255)
                                                                                                                                                                                                                                                                      cACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCTGTAATCTTAATTTCAAAGAAACTTCGCACGGTGAATAGTa  <  3:88849/196‑1 (MQ=255)
                                                                                                                                                                                                                                                                       ccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACag                                                               <  8:1055822/133‑1 (MQ=255)
                                                                                                                                                                                                                                                                             ttGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGCCt                                                <  8:805958/143‑1 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                         
AAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCTGTAATCTTAATTTCAAAGAAACTTCGCACGGTGAATAGTA  >  NC_000913/4392345‑4392802

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: