Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,971,154 G→A 5.6% intergenic (+21/‑117) ygdQ → / → ygdR UPF0053 family inner membrane protein/DUF903 family verified lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,971,1540GA 5.6% 756.4 / 21.8 231intergenic (+21/‑117)ygdQ/ygdRUPF0053 family inner membrane protein/DUF903 family verified lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (112/106);  new base A (4/9);  total (116/115)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.21e-02

Only 100 of 239 total aligned reads displayed.

TTTATTATGATGGCAGCCGTGGTGATTGCCGTAGGCGTGATGATGTTCGCCGCGCGCTCGATTGGTGATTTTGTCGAACGCCATCCTTCGGTAAAAATGCTGGCGCTCTCTTTCCTGATTCTGGTGGGCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAATTTCTGCAGTCGGACTGGCGTTTGCTGTTTCCGGGTGTTCCAGTGATTACGTCATGGCGACCAAAGATGGCCGTATGATTTTGACCGATGGAAAACCTGAAATTGATGATGATACCGGGCTGGTGAGTTATCAC  >  NC_000913/2970882‑2971426
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                
tttATTATGATGGCAGCCGTGGTGATTGCCGTAGGCGTGATGATGTTCGCCGCGCGCTCGATTGGTGATTTTGTCGAACGCCATCCTTCGGTAAAAATGCTGGCGCTCTCTTTCCTGATTCTGGTGGGCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCggg                                                                                                                                                                                                                                                                                <  2:520190/275‑1 (MQ=255)
              aGCCGTGGTGATTGCCGTAGGCGTGATGATGTTCGCCGCGCGCTCGATTGGTGATTTTGTCGAACGCCATCCTTCGGTAAAAATGCTGGCGCTCTCTTTCCTGATTCTGGTGGGCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCggg                                                                                                                                                                                                                                                                                <  2:268636/261‑1 (MQ=255)
              aGCCGTGGTGATTGCCGTAGGCGTGATGATGTTCGCCGCGCGCTCGATTGGTGATTTTGTCGAACGCCATCCTTCGGTAAAAATGCTGGCGCTCTCTTTCCTGATTCTGGTGGGCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCggg                                                                                                                                                                                                                                                                                <  2:422091/261‑1 (MQ=255)
                            ccGTAGGCGTGATGGTGTTCGCCGCGGGCTCGATTGGTGATTTTGTCGAACGCCATCCTTCGGTAAAAATGCTGGCGCTCTCTTTCCTGATTCTGGTGGGCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAGCAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTca                                                                                                                                                                                                                                                     <  2:149629/274‑1 (MQ=255)
                            ccGGAGGCGTGTTGGTGTTTGCGGGGCGTGCGATTGGGGATGTTGTCGAACGCCATCCTTCGGGAAAAATGCTGGCGCTGTCTTTCCTGATTCTGGTTGGCTTTACCATGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACGTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTAGCAAAAAAAAGAATCAGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTcg                                                                                                                                                                                                                                                     <  4:194004/274‑2 (MQ=255)
                                aGGCGTGATGATGTTCGCCGCGCGCTCGATTGGTGAGTTTGTCGAACGCCATCCTTCGGTAAAAATGCTGGCGCTCTCTTTCCTGATTCTGGTGGGCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGtt                                                                                                                                                                                                                                                <  2:74282/275‑1 (MQ=255)
                                      gatgatGTTCGCCGCGCGCTCGATTGGTGATTTTGTCGAACGCCATCCTTCGGTAAAAATGCTGGCGCTCTCTTTCCTGATTCTGGTGGGCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTca                                                                                                                                                                                                                                                     <  4:549039/264‑1 (MQ=255)
                                                   gcgcgcTCGATTGGTGATTTTGTCGAACGCCATCCTTCGGTAAAAATGCTGGCGCTCTCTTTCCTGATTCTGGTGGGCTTTACCCTGATTCTGGAAAGTTTCGACAGCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTc                                                                                                                                                                                                                              <  1:89208/274‑1 (MQ=255)
                                                    cgcgcTCGATTGGTGATTTTGTCGAACGCCATCCTTCGGTAAAAATGCTGGCGCTCTCTTTCCTGATTCTGGTGGGCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCa                                                                                                                                                                                                                             <  4:237768/274‑1 (MQ=255)
                                                                                                     ggCGCTCTCTTTCCTGATTCTGGTGGGCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCt                                                                                                                                                                                                                                 >  2:205621/1‑221 (MQ=255)
                                                                                                                          ggtggGCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCtt                                                                                                                                                                                                                                <  2:217002/201‑1 (MQ=255)
                                                                                                                             gggCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCgg                                                                                                                                                                                                                                                                                 <  7:846874/149‑1 (MQ=255)
                                                                                                                                 tttACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAgg                                                                                                                                                                                                                                                                             <  8:659324/149‑1 (MQ=255)
                                                                                                                                 tttACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAgg                                                                                                                                                                                                                                                                             <  6:393740/149‑1 (MQ=255)
                                                                                                                                      ccTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTa                                                                                                                                                                                                                                                                        >  8:273212/1‑149 (MQ=255)
                                                                                                                                      ccTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTa                                                                                                                                                                                                                                                                        >  7:412749/1‑149 (MQ=255)
                                                                                                                                      ccTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTa                                                                                                                                                                                                                                                                        <  5:667434/149‑1 (MQ=255)
                                                                                                                                       cTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTac                                                                                                                                                                                                                                                                       >  7:602209/1‑149 (MQ=255)
                                                                                                                                       cTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTac                                                                                                                                                                                                                                                                       <  5:1020009/149‑1 (MQ=255)
                                                                                                                                        tGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTact                                                                                                                                                                                                                                                                      >  6:839616/1‑149 (MQ=255)
                                                                                                                                        tGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTact                                                                                                                                                                                                                                                                      >  5:338310/1‑149 (MQ=255)
                                                                                                                                        tGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTact                                                                                                                                                                                                                                                                      >  6:304093/1‑149 (MQ=255)
                                                                                                                                            tCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTacggagt                                                                                                                                                                                                                                                                  <  5:397055/149‑6 (MQ=255)
                                                                                                                                              tGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTacgga                                                                                                                                                                                                                                                                    <  5:588581/145‑4 (MQ=255)
                                                                                                                                              tGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTacgga                                                                                                                                                                                                                                                                    >  6:588581/1‑142 (MQ=255)
                                                                                                                                                gAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCggg                                                                                                                                                                                                                                                                                >  2:239615/1‑131 (MQ=255)
                                                                                                                                                      ttCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTAccc                                                                                                                                                                                                                                                        >  5:769148/1‑149 (MQ=255)
                                                                                                                                                       tCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCt                                                                                                                                                                                                                                                       <  7:972891/149‑1 (MQ=255)
                                                                                                                                                               cACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCttt                                                                                                                                                                                                                               >  4:208258/1‑165 (MQ=255)
                                                                                                                                                                 cGTACCTGAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTaa                                                                                                                                                                                                                                             <  6:788229/149‑1 (MQ=255)
                                                                                                                                                                 cGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTaa                                                                                                                                                                                                                                             <  7:826307/149‑1 (MQ=255)
                                                                                                                                                                     ccGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGAtt                                                                                                                                                                                                                                         >  6:472422/1‑149 (MQ=255)
                                                                                                                                                                      cGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGAttt                                                                                                                                                                                                                                        >  5:886453/1‑149 (MQ=255)
                                                                                                                                                                      cGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGAttt                                                                                                                                                                                                                                        >  7:869363/1‑149 (MQ=255)
                                                                                                                                                                      cGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGAttt                                                                                                                                                                                                                                        >  7:775118/1‑149 (MQ=255)
                                                                                                                                                                       gAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGAtttt                                                                                                                                                                                                                                       >  7:518142/1‑149 (MQ=255)
                                                                                                                                                                        aaaGGTTACATCTACTTCGTGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTc                                                                                                                                                                                                                                      >  6:709760/1‑149 (MQ=255)
                                                                                                                                                                        aaaGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTc                                                                                                                                                                                                                                      >  5:873673/1‑149 (MQ=255)
                                                                                                                                                                        aaaGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTc                                                                                                                                                                                                                                      <  6:122976/149‑1 (MQ=255)
                                                                                                                                                                                 aTCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCa                                                                                                                                                                                                                             >  7:592075/1‑149 (MQ=255)
                                                                                                                                                                                     aCTTCGCGATGTTCTTATCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGat                                                                                                                                                                                                                         >  6:368718/1‑149 (MQ=255)
                                                                                                                                                                                      cTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTg                                                                                                                                                                                                                                   >  4:576887/1‑138 (MQ=255)
                                                                                                                                                                                      cTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTATGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGCCTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTACTAATTTCTGCAGTCGGTCTGGCGTTTGCTGTTTCCGGgttgtccctta                                                                                           >  4:332143/1‑265 (MQ=255)
                                                                                                                                                                                       ttCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGatat                                                                                                                                                                                                                       >  7:638431/1‑149 (MQ=255)
                                                                                                                                                                                        tCGCGATGTGCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGatata                                                                                                                                                                                                                      <  5:472422/149‑1 (MQ=255)
                                                                                                                                                                                         cgcgATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCCCAGTTTAAGATTTTCCTGCTTTCAAGatatat                                                                                                                                                                                                                     <  6:1012640/149‑1 (MQ=255)
                                                                                                                                                                                          gcgATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCttt                                                                                                                                                                                                                               <  2:576887/138‑1 (MQ=255)
                                                                                                                                                                                             aTGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACg                                                                                                                                                                                                                 >  8:386665/1‑149 (MQ=255)
                                                                                                                                                                                                  cttctCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGt                                                                                                                                                                                                            <  8:620869/149‑1 (MQ=255)
                                                                                                                                                                                                  cttctCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGt                                                                                                                                                                                                            <  7:267258/149‑1 (MQ=255)
                                                                                                                                                                                                   ttctCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGtt                                                                                                                                                                                                           <  5:936070/149‑1 (MQ=255)
                                                                                                                                                                                                      tctATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTtat                                                                                                                                                                                                        <  7:633246/149‑1 (MQ=255)
                                                                                                                                                                                                        tATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCt                                                                                                                                                                                                                                    >  4:289521/1‑119 (MQ=255)
                                                                                                                                                                                                        tATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATaa                                                                                                                                                                                                      >  6:185195/1‑149 (MQ=255)
                                                                                                                                                                                                        tATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATaa                                                                                                                                                                                                      >  8:412127/1‑149 (MQ=255)
                                                                                                                                                                                                            gCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCtt                                                                                                                                                                                                                                <  2:289521/119‑1 (MQ=255)
                                                                                                                                                                                                             cGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTAcc                                                                                                                                                                                                                                                         >  2:662835/1‑93 (MQ=255)
                                                                                                                                                                                                                  gAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGt                                                                                                                                                                                                                                                 >  1:39046/1‑96 (MQ=255)
                                                                                                                                                                                                                   aaaGCCGCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAATTTCTGCAGTCGGACTGGCGTTTGCTGTTTCCGGGTGTTCCAGTGATTACGTCATGGCGACCAAAGATGGCCGTa                                                             <  1:151736/275‑1 (MQ=255)
                                                                                                                                                                                                                     aGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATc                                                                                                                                                                                         <  7:708760/149‑1 (MQ=255)
                                                                                                                                                                                                                      gCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCggg                                                                                                                                                                                                                                                                                <  4:312352/61‑1 (MQ=255)
                                                                                                                                                                                                                      gCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTaa                                                                                                                                                                                                                                             <  3:39046/96‑1 (MQ=255)
                                                                                                                                                                                                                       ccTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAt                                                                                                                                                                                       >  7:1038142/1‑149 (MQ=255)
                                                                                                                                                                                                                        cTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAta                                                                                                                                                                                      <  6:992216/149‑1 (MQ=255)
                                                                                                                                                                                                                        cTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAta                                                                                                                                                                                      <  7:741793/149‑1 (MQ=255)
                                                                                                                                                                                                                         tCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtag                                                                                                                                                                                     <  8:869363/149‑1 (MQ=255)
                                                                                                                                                                                                                         tCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtag                                                                                                                                                                                     <  8:802309/149‑1 (MQ=255)
                                                                                                                                                                                                                           aaCTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAAc                                                                                                                                                                                                                  >  7:349559/1‑118 (MQ=255)
                                                                                                                                                                                                                           aaCTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAAc                                                                                                                                                                                                                  <  8:349559/118‑1 (MQ=255)
                                                                                                                                                                                                                           aaCTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtaggt                                                                                                                                                                                   >  6:2767/1‑149 (MQ=255)
                                                                                                                                                                                                                           aaCTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtaggt                                                                                                                                                                                   >  8:476772/1‑149 (MQ=255)
                                                                                                                                                                                                                              ttGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCt                                                                                                                                                                                                                                                       <  4:682812/78‑1 (MQ=255)
                                                                                                                                                                                                                              ttGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATc                                                                                                                                                                                         >  7:496029/1‑140 (MQ=255)
                                                                                                                                                                                                                              ttGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATc                                                                                                                                                                                         <  8:496029/140‑1 (MQ=255)
                                                                                                                                                                                                                                    cGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAAt                                                                                                                                                                          <  8:518142/149‑1 (MQ=255)
                                                                                                                                                                                                                                      caacaaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGa                                                                                                                                                                        >  7:290931/1‑149 (MQ=255)
                                                                                                                                                                                                                                       aacaaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAg                                                                                                                                                                       <  7:386665/149‑1 (MQ=255)
                                                                                                                                                                                                                                         caaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACAATTATCATAGGTAGGCCAAATGAGGa                                                                                                                                                                     <  6:910438/149‑1 (MQ=255)
                                                                                                                                                                                                                                            aaaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATaa                                                                                                                                                                  >  7:1076798/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     gCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTAGTCACTCCTTTTACCCTCACAGTTTAAGCTTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGaa                                                                                                                                                         >  7:906752/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     gCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATgca                                                                                                                                                         >  8:250512/1‑147 (MQ=255)
                                                                                                                                                                                                                                                     gCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGCATGaa                                                                                                                                                         >  6:855226/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          gATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAAt                                                                                                                                                    <  8:603999/149‑1 (MQ=255)
                                                                                                                                                                                                                                                             aCTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATggg                                                                                                                                                 <  8:314028/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               tCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATggg                                                                                                                                                 >  5:822062/1‑147 (MQ=255)
                                                                                                                                                                                                                                                                ccGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAATTTCTGCAGTCGGACTGGCGTTTGCTGt                                                                                                            >  4:26616/1‑183 (MQ=255)
                                                                                                                                                                                                                                                                 cGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATaa                                                                                                                                                                                                                   >  3:165531/1‑79 (MQ=255)
                                                                                                                                                                                                                                                                   tACGCTTTCCTGCGGGAGGGCGTACTCACTCCTTTTGCCCTCACAGTTCAAGATTTTCCTGCTTTCACGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCTAAATGAGGTTACCCGGATGAAAAAATGGGCAGTAATAATTTCTGCAGTCGGACTGGCGTTTGCTGTTTCCGGGTGTTCAAGTGTTTCCGTCATGGCTACCAAAGACGGCCGTTACATTTTGCCCGATGGAAAACTGGAAATTGATGATGTTACCGGGTTg              >  4:86935/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                   tACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGa                                                                                                                                                                                                 <  8:1052902/95‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aCGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGtaat                                                                                                                                          <  8:775118/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                     cGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGc                                                                                                                                                                                                                                  >  3:103923/1‑60 (MQ=255)
                                                                                                                                                                                                                                                                     cGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTc                                                                                                                                                                                                               <  1:165531/79‑1 (MQ=255)
                                                                                                                                                                                                                                                                     cGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGa                                                                                                                                                                        <  5:268879/118‑1 (MQ=255)
                                                                                                                                                                                                                                                                        tctcCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAAtt                                                                                                                                      >  8:20019/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         ctcCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTc                                                                                                                                                                                                                              <  1:103923/60‑1 (MQ=255)
                                                                                                                                                                                                                                                                          tcCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAATTTc                                                                                                                                    >  6:461463/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                          tcCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAATTTc                                                                                                                                    <  8:74174/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                            cTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACta                                                                                                                                                                                              <  2:565514/89‑1 (MQ=255)
                                                                                                                                                                                                                                                                              gCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATggg                                                                                                                                                 <  2:395911/132‑1 (MQ=255)
                                                                                                                                                                                                                                                                               cGGGATGGTGTACACACTTCTTTTACCCTCACAGTTTACGATTTTCCTGCTTTCAAGATATAGAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATACTTTATGCAGTCGGACTGGCGTTGGCTGTTCCCGCGTGTTCCAGTGATTACGTCATGGCGACCAAAGAGGGCCCTATGATTTTGACCGATGGAAAACCTGCACTTGATGATGATTCCGGCATGGTCAGTTATCAc  >  4:104212/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                
TTTATTATGATGGCAGCCGTGGTGATTGCCGTAGGCGTGATGATGTTCGCCGCGCGCTCGATTGGTGATTTTGTCGAACGCCATCCTTCGGTAAAAATGCTGGCGCTCTCTTTCCTGATTCTGGTGGGCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAATTTCTGCAGTCGGACTGGCGTTTGCTGTTTCCGGGTGTTCCAGTGATTACGTCATGGCGACCAAAGATGGCCGTATGATTTTGACCGATGGAAAACCTGAAATTGATGATGATACCGGGCTGGTGAGTTATCAC  >  NC_000913/2970882‑2971426

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑