Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 832,252 T→A 5.3% intergenic (+16/+216) rhlE → / ← ybiA ATP‑dependent RNA helicase/DUF1768 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913832,2520TA 5.3% 526.4 / 13.4 146intergenic (+16/+216)rhlE/ybiAATP‑dependent RNA helicase/DUF1768 family protein
Reads supporting (aligned to +/- strand):  ref base T (58/80);  new base A (4/4);  total (62/84)
Fisher's exact test for biased strand distribution p-value = 7.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.34e-01

Only 100 of 147 total aligned reads displayed.

CCGGACCCGTCAATCAAAGCCGAACCGATCCAGAACGGTCGCCAGCAACGTGGCGGCGGCGGTCGTGGGCAAGGTGGTGGTCGCGGTCAACAGCAACCACGCCGTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAGGTAGCCCAGTCGATTCT  >  NC_000913/832003‑832526
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                  
ccGGCCCCGTCAATCAAAGCCGACCCGACCCAGAACGGTCGCCAGCAACGTGGCGGCGGCGGTCGTGGGCAAGGTGGTGGTCGCGGTCAACAGCAACCACGCCGTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAAc                                                                                                                                                                                                                                                            <  2:170125/274‑1 (MQ=255)
                     gAAACGATCCAGAAGGGTGGCAGGCAACGTGGCGGCGGGGGTGGTGGGCAAGGGGGTGGTCGCGGGCAACAGCAACCACGCCGTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATgcg                                                                                                                                                                                                                                      <  4:314513/275‑1 (MQ=255)
                           aTCCAGAACGGTCGCCAGCAACGTGGCGGCGGCGGTCGTGGGCAAGGTGGTGGTCGCGGTCAACAGCAACCACGCCGTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAACCCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCAACGc                                                                                                                                                                                                                                >  2:308799/1‑275 (MQ=255)
                                   cGGTCGGCAGCAACGTGGCGGCGGCGGTCGTGGGCAAGGTGGTGGTCGCGGTCAACAGCAACCACGCCGTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAACACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAgg                                                                                                                                                                                                                        <  4:286427/275‑1 (MQ=255)
                                   cGGTCGCCCGGCACGGGGGGGGGGCGGGGNGGGGGAAGGGGGTGGCCGGGGTACCNAGCAACCACGCGGGGGGGAAGGGGGCGGAAAATCTGCAAGCGCGAAACCCGCAGAAAAACCGTCTCGCCGCCTCGGCGATNCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAgg                                                                                                                                                                                                                        <  3:60611/275‑1 (MQ=255)
                                                cGTGGCGGCGGCGGTCGTGGGCAAGGTGGTGGTCGCGGTCAACAGCAACCACGCCGTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGaa                                                                                                                                                                                                              <  2:274745/272‑1 (MQ=255)
                                                                                  gcgGTCAACAGCAACCACGCCGTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGc                                                                                                                                                                          <  2:43025/274‑1 (MQ=255)
                                                                                                ccACGCCGTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAAc                                                                                                                                                                                                                                                 >  4:263836/1‑189 (MQ=255)
                                                                                                    gCCGTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGcc                                                                                                                                                                                                                                             <  2:263836/189‑1 (MQ=255)
                                                                                                    gCCGTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTCAACTCTGCAATATATTTAATTTGCATGCTTTTGTAGGCCGGATAAGGCGGTTAcgcc                                                                                                                                                       >  2:41473/1‑275 (MQ=255)
                                                                                                       gTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTAt                                                                                                                                                                                                                                                                                  <  5:715404/149‑1 (MQ=255)
                                                                                                          gggAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGcc                                                                                                                                                                                                                                                                               >  8:493370/1‑149 (MQ=255)
                                                                                                               ggTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCATAGCCCGGCa                                                                                                                                                                                                                                                                          <  8:1395079/149‑1 (MQ=255)
                                                                                                               ggTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCATAGCCCGGCa                                                                                                                                                                                                                                                                          <  5:14282/149‑1 (MQ=255)
                                                                                                                  ggCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCa                                                                                                                                                                                                                                                                       <  5:1160593/149‑1 (MQ=255)
                                                                                                                  ggCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCa                                                                                                                                                                                                                                                                       <  6:1001497/149‑1 (MQ=255)
                                                                                                                   gcgcAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAg                                                                                                                                                                                                                                                                      <  6:1319159/149‑1 (MQ=255)
                                                                                                                    cgcAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCATAGCCCGGCAtaaaa                                                                                                                                                                                                                                                                     <  7:92831/149‑5 (MQ=255)
                                                                                                                    cgcAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCGTTTATGCCGGGCTATGCCCGGCCCCAGGCTGATGAACAAAATCAAAACTGCCTGATTCCCTACGCTTATCAGGCCTACGTGAACTCTGACATATATTTAATTTCCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATTCTGCATGTca                                                                                                                                        >  4:137274/1‑274 (MQ=38)
                                                                                                                     gcAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCa                                                                                                                                                                                                                                                                       >  6:1285105/1‑146 (MQ=255)
                                                                                                                     gcAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCa                                                                                                                                                                                                                                                                       <  5:1285105/146‑1 (MQ=255)
                                                                                                                     gcAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGc                                                                                                                                                                                                                                                                    >  5:1161047/1‑149 (MQ=255)
                                                                                                                       aaaaTCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCtg                                                                                                                                                                                                                                                                  <  5:197860/149‑1 (MQ=255)
                                                                                                                        aaaTCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCtga                                                                                                                                                                                                                                                                 <  6:896455/149‑1 (MQ=255)
                                                                                                                           tCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAg                                                                                                                                 <  4:302053/274‑1 (MQ=255)
                                                                                                                            cTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGaa                                                                                                                                                                                                                                                             <  8:945203/149‑1 (MQ=255)
                                                                                                                             tGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAAc                                                                                                                                                                                                                                                            <  7:134115/149‑1 (MQ=255)
                                                                                                                             tGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAAc                                                                                                                                                                                                                                                            <  5:856846/149‑1 (MQ=255)
                                                                                                                               cAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACaa                                                                                                                                                                                                                                                          >  8:109546/1‑149 (MQ=255)
                                                                                                                                           cTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGc                                                                                                                                                                                                                                              <  7:758285/149‑1 (MQ=255)
                                                                                                                                             gCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCt                                                                                                                                                                                                                                            >  7:1004656/1‑149 (MQ=255)
                                                                                                                                             gCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCt                                                                                                                                                                                                                                            <  5:306149/149‑1 (MQ=255)
                                                                                                                                              cAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCtga                                                                                                                                                                                                                                                                 >  6:608958/1‑127 (MQ=255)
                                                                                                                                              cAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCtga                                                                                                                                                                                                                                                                 <  5:608958/127‑1 (MQ=255)
                                                                                                                                              cAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGcc                                                                                                                                                                         >  4:318617/1‑215 (MQ=255)
                                                                                                                                                gAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATAAACAAACGCAAAACTGCCTGAt                                                                                                                                                                                                                                         >  7:151763/1‑149 (MQ=255)
                                                                                                                                                 aaaaaCCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATg                                                                                                                                                                                                                                        >  6:927645/1‑149 (MQ=255)
                                                                                                                                                  aaaaCCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGAt                                                                                                                                                                     <  2:318617/215‑1 (MQ=255)
                                                                                                                                                     aCCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGc                                                                                                                                                                                                                                              >  6:718081/1‑139 (MQ=255)
                                                                                                                                                       cGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTAc                                                                                                                                                                                                                                  >  5:199387/1‑149 (MQ=255)
                                                                                                                                                         tctcGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGc                                                                                                                                                                                                                                <  8:302625/149‑1 (MQ=255)
                                                                                                                                                             gccgccTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTAt                                                                                                                                                                                                                            >  8:657227/1‑149 (MQ=255)
                                                                                                                                                               cgccTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCa                                                                                                                                                                                                                          <  6:1024975/149‑1 (MQ=255)
                                                                                                                                                                gccTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAg                                                                                                                                                                                                                         >  5:947747/1‑149 (MQ=255)
                                                                                                                                                                gccTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCAGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAg                                                                                                                                                                                                                         <  6:1331516/149‑1 (MQ=255)
                                                                                                                                                                 ccTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGTTGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTTTGGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTTTTAGCCATTTGAGAAAAAGGGAAAAAATTTAtt                                                                                          >  2:146394/1‑274 (MQ=255)
                                                                                                                                                                   tCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGcc                                                                                                                                                                                                                      >  5:1096289/1‑149 (MQ=255)
                                                                                                                                                                    cGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCt                                                                                                                                                                                                                     <  8:993207/149‑1 (MQ=255)
                                                                                                                                                                        gATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGc                                                                                                                                                                                                                                              >  8:85321/1‑120 (MQ=255)
                                                                                                                                                                        gATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGc                                                                                                                                                                                                                                              <  7:85321/120‑1 (MQ=255)
                                                                                                                                                                          tGCCAAACCGGCAGGCGACCAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGTGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCTCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTAcc                                                                                 <  2:103146/275‑1 (MQ=255)
                                                                                                                                                                           gCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCa                                                                                                                                                                                                                          >  4:367746/1‑137 (MQ=255)
                                                                                                                                                                            ccAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAAc                                                                                                                                                                                                             <  7:203677/149‑1 (MQ=255)
                                                                                                                                                                             cAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTAt                                                                                                                                                                                                                                                                                  <  6:972523/79‑1 (MQ=255)
                                                                                                                                                                                       ggCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAatata                                                                                                                                                                                                  <  8:145068/149‑1 (MQ=255)
                                                                                                                                                                                        gCGAACAACAACGTCGCCGCCGTCCGCGTAAGCCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAatatat                                                                                                                                                                                                 >  8:892878/1‑149 (MQ=255)
                                                                                                                                                                                         cGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATAtt                                                                                                                                                                                                <  6:1453273/149‑1 (MQ=255)
                                                                                                                                                                                         cGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATAtt                                                                                                                                                                                                <  6:1262065/149‑1 (MQ=255)
                                                                                                                                                                                          gAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTg                                                                                                                                                                                               >  6:226784/1‑149 (MQ=255)
                                                                                                                                                                                            acaacaacGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGaa                                                                                                                                                                                             <  6:543888/149‑1 (MQ=255)
                                                                                                                                                                                            acaacaacGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGaa                                                                                                                                                                                             <  7:474447/149‑1 (MQ=255)
                                                                                                                                                                                              aacaacGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAAtt                                                                                                                                                                                           >  8:719749/1‑149 (MQ=255)
                                                                                                                                                                                               acaacGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAAttt                                                                                                                                                                                          <  7:1109391/149‑1 (MQ=255)
                                                                                                                                                                                                caacGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTTCGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTg                                                                                                                                                                                         <  7:892878/149‑1 (MQ=255)
                                                                                                                                                                                                        ccgccgTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAACACTTTTCAGCAATTTGAAAACACGGGAAAAAATTTATCTGTTTTTCCGCCCATAAGACACCGCAACGttccctcc                                                     >  4:42761/1‑267 (MQ=255)
                                                                                                                                                                                                         cgccgTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCg                                                                                                                                                                                                                                                                             <  5:952116/56‑1 (MQ=255)
                                                                                                                                                                                                         cgccgTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCtttt                                                                                                                                                                                <  5:849535/149‑1 (MQ=255)
                                                                                                                                                                                                             gTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAg                                                                                                                                                                            <  6:114252/149‑1 (MQ=255)
                                                                                                                                                                                                              tCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGcc                                                                                                                                                                                                                      >  8:574648/1‑106 (MQ=255)
                                                                                                                                                                                                              tCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGcc                                                                                                                                                                                                                      <  7:574648/106‑1 (MQ=255)
                                                                                                                                                                                                              tCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAgg                                                                                                                                                                           >  8:794376/1‑149 (MQ=255)
                                                                                                                                                                                                                cgcgTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCaa                                                                                                                                                                                                      >  2:121593/1‑120 (MQ=255)
                                                                                                                                                                                                                cgcgTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGAATGCTTTTGTAGGcc                                                                                                                                                                         <  6:1081086/149‑1 (MQ=255)
                                                                                                                                                                                                                   gTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGa                                                                                                                                                                      <  5:800608/149‑1 (MQ=255)
                                                                                                                                                                                                                      aaCCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATaa                                                                                                                                                                   >  7:1093480/1‑149 (MQ=255)
                                                                                                                                                                                                                        ccTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGcc                                                                                                                                                                                                                                             >  8:778038/1‑73 (MQ=255)
                                                                                                                                                                                                                        ccTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGcc                                                                                                                                                                                                                                             <  7:778038/73‑1 (MQ=255)
                                                                                                                                                                                                                         cTGCCGCTGCGCAGTAATCTTTTATGCCGGGTTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGc                                                                                                                                                                <  7:994757/149‑1 (MQ=255)
                                                                                                                                                                                                                            ccGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGtt                                                                                                                                                             <  6:345447/149‑1 (MQ=255)
                                                                                                                                                                                                                             cGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGttt                                                                                                                                                            <  5:1076585/149‑1 (MQ=255)
                                                                                                                                                                                                                                   gcAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTAcgccg                                                                                                                                                      >  7:493309/1‑149 (MQ=255)
                                                                                                                                                                                                                                     aGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCa                                                                                                                                                    >  7:299958/1‑149 (MQ=255)
                                                                                                                                                                                                                                     aGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGAAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATGCGGCGTACAcgcctta                                                                                                                                                    <  8:1352291/149‑4 (MQ=25)
                                                                                                                                                                                                                                       tAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATc                                                                                                                                                  >  8:1288981/1‑149 (MQ=255)
                                                                                                                                                                                                                                       tAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATc                                                                                                                                                  >  7:306576/1‑149 (MQ=255)
                                                                                                                                                                                                                                         aTCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCa                                                                                                                                                                                                                          >  8:1387802/1‑75 (MQ=255)
                                                                                                                                                                                                                                         aTCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCa                                                                                                                                                                                                                          <  7:1387802/75‑1 (MQ=255)
                                                                                                                                                                                                                                             tttATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAt                                                                                                                                            <  5:844671/149‑1 (MQ=255)
                                                                                                                                                                                                                                                aTGCCGGGCATAGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATa                                                                                                                                                                    >  6:203563/1‑122 (MQ=255)
                                                                                                                                                                                                                                                aTGCCGGGCATAGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTc                                                                                                                                         >  8:501497/1‑149 (MQ=255)
                                                                                                                                                                                                                                                 tGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAAtt                                                                                                                                                                                           >  1:20518/1‑98 (MQ=255)
                                                                                                                                                                                                                                                 tGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAttt                                                                                                                                          <  5:552677/147‑1 (MQ=255)
                                                                                                                                                                                                                                                  gCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTcac                                                                                                                                       >  7:908577/1‑149 (MQ=255)
                                                                                                                                                                                                                                                    cGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCacaa                                                                                                                                     <  7:1290203/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     gggCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCa                                                                                                                                                                                       <  3:20518/98‑1 (MQ=255)
                                                                                                                                                                                                                                                      ggCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCacaaca                                                                                                                                   >  8:153741/1‑149 (MQ=255)
                                                                                                                                                                                                                                                       gCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCAcaacaa                                                                                                                                  <  8:562975/149‑1 (MQ=255)
                                                                                                                                                                                                                                                       gCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCAcaacaa                                                                                                                                  <  5:321193/149‑1 (MQ=255)
                                                                                                                                                                                                                                                        cTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAg                                                                                                                                 >  7:630192/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         tATGCTAGGCATCAGACTGAAGAACAAACGGAAAACAGCCGTGTGCGTTATGCGTATCAGGCACACGTGTACACTTCAATATAGTTAATTTGCATGCTTTTGTAGTCCGGATAAGGGGTTTACGCTGCATCCGGAATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTAACCACCCAGAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAGGTAGCCCAGTCGATTCt  <  4:146394/275‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                  
CCGGACCCGTCAATCAAAGCCGAACCGATCCAGAACGGTCGCCAGCAACGTGGCGGCGGCGGTCGTGGGCAAGGTGGTGGTCGCGGTCAACAGCAACCACGCCGTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAGGTAGCCCAGTCGATTCT  >  NC_000913/832003‑832526

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: