Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 151,522 G→T 9.8% A26A (GCC→GCA yadK ← putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913151,5220GT 9.8% 540.7 / 36.8 162A26A (GCC→GCAyadKputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base G (70/76);  new base T (7/9);  total (77/85)
Fisher's exact test for biased strand distribution p-value = 7.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-02

Only 100 of 165 total aligned reads displayed.

AATATTTGCCTAGAGCCAGCAGCCATGACACCGCTAAGGGCTCTAAAATCAATATGTAACTGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACA  >  NC_000913/151250‑151788
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                          
aaaatTTGCCTAGAGCCAGCAGCCATGACACCGCTAAGGGCTCTAAAATCAATATGTAACTGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATgg                                                                                                                                                                                                                                                                           <  1:132005/271‑1 (MQ=255)
        ccTAGAGCCAGCAGCCATGACACCGCTAATGGCTCTAAAATCAATATGTAACTGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGTTACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCTCCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCGATGGCTGGACAgg                                                                                                                                                                                                                                                                  <  1:155615/275‑1 (MQ=255)
             agCCAGCAGCCATGACACCGCTAAGGGCTCTAAAATCAATATGTAACTGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTCCATTGGCACCATGGTCAATATCCTTTCCAGCAATGGCTGGACAGGCGATg                                                                                                                                                                                                                                                             >  1:24903/1‑275 (MQ=255)
                                                  aaTATTTAACTGATTTAGCACATCGCCAGTTCCATGTAATGCACAGTCGGATAAAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCATGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtca                                                                                                                                                                                                                        <  1:12267/275‑1 (MQ=255)
                                                   atatGTAACTGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGAt                                                                                                                                                                                                                                                              >  4:461094/1‑236 (MQ=255)
                                                       gTAACTGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCaa                                                                                                                                                                                                                                                          <  2:461096/236‑1 (MQ=255)
                                                          aCTGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCAACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGa                                                                                                                                                                                                                                                                      >  4:176607/1‑221 (MQ=255)
                                                                    cacaTCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAgg                                                                                                                                                                                                                                                                  <  2:82012/215‑1 (MQ=255)
                                                                      caTCGCCATTCCCTTGAAATTANNATTCGNATACACTTATAGTGAACTCATTTCCACCTTGTTGATAACTTTCAGGAGTAACTTTACCTGAAAAATATCCCCCCCCAACGACGCCAAAATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTTCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGAAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGa                                                                                                                                                                                                    <  3:24903/275‑1 (MQ=255)
                                                                       aTCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATAGTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGa                                                                                                                                                                                                    <  2:725026/274‑1 (MQ=255)
                                                                                   cTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCACAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAACCAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGtt                                                                                                                                                                                       >  2:279549/1‑275 (MQ=255)
                                                                                              aGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGa                                                                                                                                                                                                                                                                      >  4:566122/1‑185 (MQ=255)
                                                                                                          tGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCANTGGCTGGACAGGCGATGCAAAAAANTANAGACAGAAACAGAaanatnt                                                                                                                                                                                                                              >  4:807495/1‑211 (MQ=255)
                                                                                                          tGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGCCAGAAACAGACTTATTCTCTTACTCAGCTTACGTTTGGTGGGATGCCATTATAGTTCCTGAATCTGTAAGTCTgttagg                                                                                                                                                                >  4:590969/1‑274 (MQ=255)
                                                                                                          tGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAatt                                                                                                                                                                                                                                  >  2:235063/1‑209 (MQ=255)
                                                                                                                  aaCTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATg                                                                                                                                                                                                                                                             >  2:184724/1‑174 (MQ=255)
                                                                                                                  aaCTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGt                                                                                                                                                         <  2:464303/274‑1 (MQ=255)
                                                                                                                  aaCTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTAAGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGGTATAgg                                                                                                                                                          >  1:182788/1‑273 (MQ=255)
                                                                                                                    ctctTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTAc                                                                                                                                                      >  2:641066/1‑275 (MQ=255)
                                                                                                                      ctTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCaaa                                                                                                                                                                                                                                                         <  4:184723/174‑1 (MQ=255)
                                                                                                                       ttttCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAAAANTAATGACAGAAACAGNATTATtctc                                                                                                                                                                                                                            >  4:780177/1‑202 (MQ=255)
                                                                                                                           ccACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtca                                                                                                                                                                                                                        <  2:780179/202‑1 (MQ=255)
                                                                                                                            cACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATg                                                                                                                                                                                                                                                                            <  8:650907/149‑1 (MQ=255)
                                                                                                                                 ggTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAACATCAATATTACCCTGGTTACTGATTCCGCGTTTGCAGGCGCTGTTTTTTACATTGGCAACAAGGTCAATATCTTGTCCAGCAATGGCTGGTCAGGCGATGCAAAACAGTAATGacagaaacaga                                                                                                                                                                                                                                     >  2:290547/1‑183 (MQ=255)
                                                                                                                                   ttGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGa                                                                                                                                                                                                                                                                      <  6:404364/148‑1 (MQ=255)
                                                                                                                                     gATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACa                                                                                                                                                                                                                                                                    <  6:108951/148‑1 (MQ=255)
                                                                                                                                                 gAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAAc                                                                                                                                                                                                                                                       >  8:162874/1‑149 (MQ=255)
                                                                                                                                                  aGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACa                                                                                                                                                                                                                                                      <  8:885735/149‑1 (MQ=255)
                                                                                                                                                       cATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAAt                                                                                                                                                                                                                                                 <  8:721874/149‑1 (MQ=255)
                                                                                                                                                         ttACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATg                                                                                                                                                                                                                                                <  5:511873/148‑1 (MQ=255)
                                                                                                                                                              tGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaa                                                                                                                                                                                                                                          >  8:146289/1‑149 (MQ=255)
                                                                                                                                                              tGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaa                                                                                                                                                                                                                                          >  7:757664/1‑149 (MQ=255)
                                                                                                                                                                    atatCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaacagaa                                                                                                                                                                                                                                    <  7:670186/149‑1 (MQ=255)
                                                                                                                                                                         ccACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAattat                                                                                                                                                                                                                                <  6:37274/148‑1 (MQ=255)
                                                                                                                                                                            ccccAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATtctc                                                                                                                                                                                                                            <  7:376260/149‑1 (MQ=255)
                                                                                                                                                                             cccAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTctct                                                                                                                                                                                                                           <  7:146289/149‑1 (MQ=255)
                                                                                                                                                                              ccAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaacagaa                                                                                                                                                                                                                                    >  6:755260/1‑139 (MQ=255)
                                                                                                                                                                                       ccAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAatt                                                                                                                                                                                                                                  >  4:253184/1‑132 (MQ=255)
                                                                                                                                                                                       ccAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGa                                                                                                                                                                                                          >  2:765494/1‑156 (MQ=255)
                                                                                                                                                                                         aaGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAAGAGAATTATTCTCTTCATCAGCTTa                                                                                                                                                                                                                >  6:595213/1‑148 (MQ=255)
                                                                                                                                                                                         aaGAGCAATCTTACCCGGTTTCTTGATTCCGCCGTTGCAGGTGCTGTTTGTTACATTGGTACCAAGGTAAATAACCGGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTGCCTGAATCTGTAAGTCTGTTCGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGAc                                                                                 <  4:306823/275‑1 (MQ=255)
                                                                                                                                                                                           gATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTgg                                                                                                                                                                                                      <  4:765492/156‑1 (MQ=255)
                                                                                                                                                                                            aTCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGtt                                                                                                                                                                                                            <  7:611194/149‑1 (MQ=255)
                                                                                                                                                                                            aTCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGt                                                                                                                                                                                                             <  5:619071/148‑1 (MQ=255)
                                                                                                                                                                                                    tACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGa                                                                                                                                                                                                    <  8:239466/149‑1 (MQ=255)
                                                                                                                                                                                                     aCCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGAt                                                                                                                                                                                                   <  7:162874/149‑1 (MQ=255)
                                                                                                                                                                                                      cccTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAACAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACt                                                                     >  2:135921/1‑274 (MQ=255)
                                                                                                                                                                                                        cTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTCTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCa                                                                                                                                                                                                >  7:571248/1‑149 (MQ=255)
                                                                                                                                                                                                          ggTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATgg                                                                                                                                                                                                                                                                           <  3:70622/72‑1 (MQ=255)
                                                                                                                                                                                                           gTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGt                                                                                                                                                                                                        >  3:126075/1‑138 (MQ=255)
                                                                                                                                                                                                           gTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGt                                                                                                                                                                                                        >  2:186930/1‑138 (MQ=255)
                                                                                                                                                                                                           gTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGt                                                                                                                                                                        >  1:46654/1‑170 (MQ=255)
                                                                                                                                                                                                           gTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGt                                                                                                                                      >  2:316046/1‑204 (MQ=255)
                                                                                                                                                                                                            ttACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACAg                                                                                                                                                                                                                                                                   >  8:456250/1‑78 (MQ=255)
                                                                                                                                                                                                               cTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAatta                                                                                                                                                                                                                                 <  2:164447/109‑1 (MQ=255)
                                                                                                                                                                                                               cTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGa                                                                                                                                                                                                    <  4:186929/138‑1 (MQ=255)
                                                                                                                                                                                                               cTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCAtatata                                                                                                                                                                                          >  5:356396/1‑148 (MQ=255)
                                                                                                                                                                                                               cTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCAtatata                                                                                                                                                                                          <  5:541739/148‑1 (MQ=255)
                                                                                                                                                                                                               cTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAg                                                                                                                                  <  4:316044/204‑1 (MQ=255)
                                                                                                                                                                                                                 gATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGt                                                                                                                                                                                        <  5:907235/148‑1 (MQ=255)
                                                                                                                                                                                                                      cGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTg                                                                                                                                                                                   <  5:612400/148‑1 (MQ=255)
                                                                                                                                                                                                                         ttttGCAGGTGCTGTTTGTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTa                                                                                                       <  1:110608/221‑1 (MQ=255)
                                                                                                                                                                                                                                  tGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaacaga                                                                                                                                                                                                                                     >  2:17786/1‑86 (MQ=255)
                                                                                                                                                                                                                                    cTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgt                                                                                                                                                                    <  8:465509/149‑1 (MQ=255)
                                                                                                                                                                                                                                      gTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAatta                                                                                                                                                                                                                                 <  4:17786/86‑1 (MQ=255)
                                                                                                                                                                                                                                      gTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgtt                                                                                                                                                                   <  5:386417/148‑1 (MQ=255)
                                                                                                                                                                                                                                         ttttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgttagt                                                                                                                                                                <  5:595213/148‑1 (MQ=255)
                                                                                                                                                                                                                                            tACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGACTCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCAGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGCCCATGCGTGCACTCATCGGATAATTAACGACCTCATTGGCTGCAGTGGAGAt                              >  1:63333/1‑275 (MQ=255)
                                                                                                                                                                                                                                                 tGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTa                                                                                                                                                       <  8:757664/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           gTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAAtg                                                                                                                                              <  5:11795/148‑1 (MQ=255)
                                                                                                                                                                                                                                                             cAATATCCTGTCCAGCAATGGCTGGACAGGCGATGc                                                                                                                                                                                                                                                            >  4:40998/1‑36 (MQ=255)
                                                                                                                                                                                                                                                              aaTATCCTTTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGc            <  2:590971/275‑1 (MQ=255)
                                                                                                                                                                                                                                                              aaTATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGc                                                                                                                                          >  7:350036/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               atatCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGtt                                                                                                                                                        <  5:775931/134‑1 (MQ=255)
                                                                                                                                                                                                                                                               atatCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGt                                                                                                                                                         <  1:186328/133‑1 (MQ=255)
                                                                                                                                                                                                                                                                tatCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTg                                                                                                                                        >  7:809805/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 atCCTGTCCAGCAATGGCTGGACAGGCGATGCaaaa                                                                                                                                                                                                                                                        <  2:40999/36‑1 (MQ=255)
                                                                                                                                                                                                                                                                  gccTGTCCAGCCATTGCTGGACAGGCGATGCAAAGCAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAGGTCTGTTAGTTATAGGTTACTGTAAATGTTGCTGGt                                                                                                                                      >  7:601360/2‑149 (MQ=255)
                                                                                                                                                                                                                                                                     tGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTa                                                                                                                                                                           >  6:618550/1‑109 (MQ=255)
                                                                                                                                                                                                                                                                     tGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGt                                                                                                                                                                                        >  3:51053/1‑96 (MQ=255)
                                                                                                                                                                                                                                                                     tGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTa                                                                                                                          >  2:349595/1‑158 (MQ=255)
                                                                                                                                                                                                                                                                      gTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgt                                                                                                                                                                    >  6:758594/1‑115 (MQ=255)
                                                                                                                                                                                                                                                                      gTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgt                                                                                                                                                                    <  5:758594/115‑1 (MQ=255)
                                                                                                                                                                                                                                                                      gTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCCGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTCCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACCANCAGGCGTGCACTCTTCGGATAATCAACGACCGAATTTGCTGCAGGGGAGTATTTTCTCTCTTGTGATGCTTGTGttt    >  3:170952/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                      gTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGt                                                                                                                                                                                                        >  2:586525/1‑79 (MQ=255)
                                                                                                                                                                                                                                                                      gTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGt                                                                                                                                                                                                        >  1:229852/1‑79 (MQ=255)
                                                                                                                                                                                                                                                                      gTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGt                                                                                                                                                                                                        >  3:18797/1‑79 (MQ=255)
                                                                                                                                                                                                                                                                      gTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCa                                                                                                                                   >  5:824142/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                      gTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGCGGAGATTGTTACCTCTTGTGATGCTGGTGTTa    >  4:153808/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                         cAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGa                                                                                                                                                                                                    <  2:500299/80‑1 (MQ=255)
                                                                                                                                                                                                                                                                         cAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGt                                                                                                                                                         <  4:600622/123‑1 (MQ=255)
                                                                                                                                                                                                                                                                         cAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCg                                                                                                                      <  4:349593/158‑1 (MQ=255)
                                                                                                                                                                                                                                                                          aGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGa                                                                                                                                                                                                    <  3:229851/79‑1 (MQ=255)
                                                                                                                                                                                                                                                                          aGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGa                                                                                                                                                                                                    <  1:18797/79‑1 (MQ=255)
                                                                                                                                                                                                                                                                          aGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGaaa                                                                                                                               <  5:455409/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                           gCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTa                                                                                                                                                                           >  6:461948/1‑103 (MQ=255)
                                                                                                                                                                                                                                                                           gCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAAc                                                                                                                             >  7:747206/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                            cAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAAc                                                                                                                             >  6:613413/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                             aaTGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACt                                                                                                                            >  6:775206/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                               tGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTa                                                                                                                          <  6:824142/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                          
AATATTTGCCTAGAGCCAGCAGCCATGACACCGCTAAGGGCTCTAAAATCAATATGTAACTGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACA  >  NC_000913/151250‑151788

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