Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 151,688 C→G 9.6% V182L (GTA→CTA)  yadL ← putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913151,6880CG 9.6% 1607.3 / 33.6 508V182L (GTA→CTA) yadLputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base C (259/191);  new base G (3/55);  total (262/246)
Fisher's exact test for biased strand distribution p-value = 9.33e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 526 total aligned reads displayed.

AAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGCTAAATTGTGTTGCACTCCCTGTTGCCGGAGTTACAGTGCTCCAGA  >  NC_000913/151539‑151837
                                                                                                                                                     |                                                                                                                                                     
aaaaCAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTAc                                                                                                                                                       >  2:78741/1‑150 (MQ=255)
aaaaCAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTAc                                                                                                                                                       <  1:2663387/150‑1 (MQ=255)
 aaaCAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACg                                                                                                                                                      >  1:1002661/1‑150 (MQ=255)
  aaCAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACgg                                                                                                                                                     <  2:1979778/150‑1 (MQ=255)
  aaCAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACgg                                                                                                                                                     <  2:9275183/150‑1 (MQ=255)
       tactGACAGACACAGACTTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTTAGTCTGTTAGTTATCGGTTACGGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTa                                                                                                                                                <  1:6452977/147‑1 (MQ=255)
       tAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTa                                                                                                                                                <  1:722606/150‑1 (MQ=255)
        aaTGACAGAAACTGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAATCTTACCCGCAGTGACATTGTTTACGGTTTTAt                                                                                                                                               <  2:5352467/150‑1 (MQ=255)
              agaaacagaaTTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTAttttttt                                                                                                                                         <  1:2234588/150‑1 (MQ=255)
              agaaacagaaTTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTAttttttt                                                                                                                                         <  1:9277656/150‑1 (MQ=255)
               gaaacagaaTTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTTAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTc                                                                                                                                        <  1:1687818/150‑1 (MQ=255)
               gaaacagaaTTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCACAGTGACATTGTTTACGGTTTTATTTTTTTc                                                                                                                                        <  2:4751026/150‑1 (MQ=255)
                 aacagaaTTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCa                                                                                                                                      >  1:1976291/1‑150 (MQ=255)
                    agaaTTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGa                                                                                                                                   >  2:7360621/1‑150 (MQ=255)
                     gaaTTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGAc                                                                                                                                  <  2:1675550/150‑1 (MQ=255)
                      aaTTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGAcc                                                                                                                                 >  1:8894232/1‑150 (MQ=255)
                       attattCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCa                                                                                                                                <  2:5782268/150‑1 (MQ=255)
                         tattCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg                                                                                                                              >  2:3806180/1‑150 (MQ=255)
                         tattCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg                                                                                                                              <  2:2052322/150‑1 (MQ=255)
                           ttCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCg                                                                                                                            >  1:3874818/1‑150 (MQ=255)
                            tctcTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGt                                                                                                                           >  2:5553089/1‑150 (MQ=255)
                            tctcTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGt                                                                                                                           <  1:6107965/150‑1 (MQ=255)
                             ctctTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTg                                                                                                                          <  1:4400336/150‑1 (MQ=255)
                               ctTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCa                                                                                                                        >  1:7764664/1‑150 (MQ=255)
                                  catcaGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTc                                                                                                                     >  1:5017350/1‑150 (MQ=255)
                                    tcaGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCAt                                                                                                                   <  1:8437045/150‑1 (MQ=255)
                                    tcaGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCAt                                                                                                                   >  1:8190507/1‑150 (MQ=255)
                                    tcaGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCAt                                                                                                                   <  1:3754315/150‑1 (MQ=255)
                                        cTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGa                                                                                                               <  2:618893/150‑1 (MQ=255)
                                        cTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGa                                                                                                               <  2:4366994/150‑1 (MQ=255)
                                         ttACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGat                                                                                                              >  1:1349313/1‑150 (MQ=255)
                                         ttACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGat                                                                                                              <  2:5516932/150‑1 (MQ=255)
                                          tACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGata                                                                                                             <  2:7115394/150‑1 (MQ=255)
                                           aCGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGc                                                                                                                             <  2:2628027/133‑1 (MQ=255)
                                             gTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGataata                                                                                                          >  1:7077447/1‑150 (MQ=255)
                                             gTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGataata                                                                                                          >  1:2770503/1‑150 (MQ=255)
                                              ttGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGAtaataa                                                                                                         <  2:3359953/150‑1 (MQ=255)
                                              ttGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGAtaataa                                                                                                         <  1:6278670/150‑1 (MQ=255)
                                               tGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATaaa                                                                                                        <  1:6659337/150‑1 (MQ=255)
                                                gAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAAc                                                                                                       >  2:1747296/1‑150 (MQ=255)
                                                       aTGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCa                                                                                                >  2:3206449/1‑150 (MQ=255)
                                                          cttatataGTTCCTGAATCTGTATGTCTGTTTGTTATTGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACTTTGTTTACGGTTTTATTTTTTTCCACGACCAGGGGTGCACTCTTGGGATAATAAACGACCGCATTg                                                                                             <  2:8456965/148‑1 (MQ=255)
                                                            tatataGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGc                                                                                           >  1:6928904/1‑150 (MQ=255)
                                                                taGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCa                                                                                       >  2:5695455/1‑150 (MQ=255)
                                                                taGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCa                                                                                       >  2:6889596/1‑150 (MQ=255)
                                                                 aGTTCCTGAATCTGTAAGTCTGTTAGTTATCGGTTACTTTAAATGTGGCTGGTGCTGAAAACTTACCCGCTGTGACTTTTTTTAGGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAg                                                                                      <  2:7393874/150‑1 (MQ=255)
                                                                  gTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTAtttttt                                                                                                                                          >  2:4155352/1‑97 (MQ=255)
                                                                      cTGAATCTGTAAGTCTGTTAGTTATCGGTTAGTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCCGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGgag                                                                                 <  2:463824/150‑1 (MQ=255)
                                                                          aTCTGTACGTCGGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTCCCCGCAGTGACATTTTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTg                                                                             <  2:6011565/150‑1 (MQ=255)
                                                                            cTGTAAGTCTGTTAGTTATCGGTTACTGTAAAGGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTTTTTACGTTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGTATAATAAACGACCGCATTGGCTGCAGTGGAGATTGtt                                                                           <  2:3586132/150‑1 (MQ=255)
                                                                             tGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGTTGCAGTGGAGATTGTTa                                                                          >  1:1890786/1‑150 (MQ=255)
                                                                              ttatgtatGTTAGTTATAGTTTACTGTAATTGTGGCTGGTGCATAAAACTTACCCGCAGTGACTTTGTTTAGGTTTTTATTTTTGTCCACGACCAGGGGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTAc                                                                         <  2:8779705/143‑1 (MQ=255)
                                                                              gTAATTCTTTTAGTTATAGGTTACTGTAATTGTGGCTGGTGCAGAAAACTTACCCGCCGTTACATTGTTTAGGTTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTAc                                                                         <  2:8063464/150‑1 (MQ=255)
                                                                               tAAGTCTGTTAGTTATAGGTTACTGTAATTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCCCGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTAcc                                                                        <  2:6250978/150‑1 (MQ=255)
                                                                               tAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTagc                                                                        >  2:8607842/1‑148 (MQ=255)
                                                                               tAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTAcc                                                                        >  1:532403/1‑150 (MQ=255)
                                                                                   tCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCtt                                                                    >  2:1159611/1‑150 (MQ=255)
                                                                                   tCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCtt                                                                    >  2:8907806/1‑150 (MQ=255)
                                                                                    cTGTTAGTTCTCGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTtg                                                                   <  2:6705236/150‑1 (MQ=255)
                                                                                     tGTTAGTTATAGGTTTATTAAATTGTGGCGGGTGCCGAAAACTTACCCGCTGTGACATTGTTTAGGGTTTTATTTTTTACCACGCCCTGGGGTGCTCTCTTGGGTTAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTtgt                                                                  <  2:8009875/150‑1 (MQ=255)
                                                                                     tGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTtgt                                                                  >  2:1280978/1‑150 (MQ=255)
                                                                                      gttagttaTAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTtgtg                                                                 >  2:40487/1‑150 (MQ=255)
                                                                                             ataGGTTACTTTAAATGTGGCTGGTGCATAAAACTAACCCGCTGTGACATTTTTTAGGTTTTTTTTTTTTTCCACGACCCGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTgg                                                          <  2:1004979/150‑1 (MQ=255)
                                                                                              ttgtttACTGTAATTGTGGGTGGTGCATAAAACTAACCCGCCGTGACTTTTTTTTGGTTTTTATTTTTTTCCCCCCCCCGGGGTGCACTCATCGTATAAAACACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGgt                                                         <  1:465447/146‑1 (MQ=255)
                                                                                              taGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGgt                                                         <  2:2770503/150‑1 (MQ=255)
                                                                                                 gTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGtt                                                      >  2:5124407/1‑150 (MQ=255)
                                                                                                 gTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGtt                                                      >  1:4430394/1‑150 (MQ=255)
                                                                                                 gTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGtt                                                      >  1:3591365/1‑150 (MQ=255)
                                                                                                    aCTGTAATTGTGGCTGGTGCATAAAACTTACCCGCCGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTaca                                                   <  2:303543/150‑1 (MQ=255)
                                                                                                    aCTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTaca                                                   <  2:7801569/150‑1 (MQ=255)
                                                                                                      tGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTAcaca                                                 >  1:8051421/1‑150 (MQ=255)
                                                                                                      tGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTAcaca                                                 >  1:5555550/1‑150 (MQ=255)
                                                                                                        tAATTGTGGCTGGTGAACAAAACTTACCCGCAGTGAAATTGTTTAGGTTTTTATTTTTTCCCACGACCCGGCGTGCACTCATCGTATAAAACACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGc                                               <  1:1218827/150‑1 (MQ=255)
                                                                                                        tAAATTTGGATGGTCCACAAAACTAACCCCCAGTCACATTTTTTACGTTTTTATATTTTACCACCACCAGGCGTCCACAAATCGTAAAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGc                                               <  2:8128713/150‑1 (MQ=255)
                                                                                                        tAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGc                                               >  1:6593215/1‑150 (MQ=255)
                                                                                                                gCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGCTAAATtgt                                       >  1:9163614/1‑150 (MQ=255)
                                                                                                                gCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGCTAAATtgt                                       >  2:9096156/1‑150 (MQ=255)
                                                                                                                 cTGGTGCATAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGCTAAATtgtg                                      <  2:747965/150‑1 (MQ=255)
                                                                                                                 cTGGTGCAGAAAACTTACCCGCAGTGAAATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGCTAAATtgtg                                      <  2:1488014/150‑1 (MQ=255)
                                                                                                                    gTGCAGAAAACTTACCCGCAGTTACATTGTTTAGGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGCTAAATTGTGTTg                                   <  2:2902104/150‑1 (MQ=255)
                                                                                                                     tGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGCTAAATTGTGTTGc                                  >  2:4694995/1‑150 (MQ=255)
                                                                                                                        aGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGCTAAATTGTGTTGCACt                               >  1:2052716/1‑150 (MQ=255)
                                                                                                                        aGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGGGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGCTAAATTGTGTTGCACt                               >  1:5760801/1‑150 (MQ=255)
                                                                                                                           aaaCTTCCCCGCTTTTAAATTTTTTACGTTTTTATTTTTTCCCACGACCTGGCGTGCACTCTTCGCATAATACACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGCTAAATTGTGTTGCACTccc                            <  1:8810188/150‑1 (MQ=255)
                                                                                                                           aaaCTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGCTAAATTGTGTTGCACTccc                            >  2:4875411/1‑150 (MQ=255)
                                                                                                                               ttACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGCTAAATTGTGTTGCACTCCCTGtt                        >  2:7758007/1‑150 (MQ=255)
                                                                                                                               ttACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGCTAAATTGTGTTGCACTCCCTGtt                        >  2:8532645/1‑150 (MQ=255)
                                                                                                                               ttACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGGGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGCTAAATTGTGTTGCACTCCCTGtt                        >  1:4307578/1‑150 (MQ=255)
                                                                                                                                tACCCGCATTTACTTTGTTTAGGGTTTTTTTTTTTCCCACGACCTGGGTTGCACTCATCGTTTAATACACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGCTAAATTGTGTTGCACTCCCTGTTg                       <  1:3975427/150‑1 (MQ=255)
                                                                                                                                   ccGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGCTAAATTGTGTTGCACTCCCTGTTGCCg                    >  2:6386329/1‑150 (MQ=255)
                                                                                                                                   ccGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGCTAAATTGTGTTGCACTCCCTGTTGCCg                    >  1:3329353/1‑150 (MQ=255)
                                                                                                                                    cGCATTCACTTTGTTTACGTTTTTTTTTTTTTCCACGACCAGGCGTGCACTCATCGTATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGCTAAATTGTGTTGCACTCCCTGTTGCCgg                   <  1:8574262/150‑1 (MQ=255)
                                                                                                                                     gctgTGACATTTTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACACATCGTATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGCTAAATTGTGTTGCACTCCCTGTTGCCGGa                  <  2:6545318/147‑1 (MQ=255)
                                                                                                                                       aGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGCTAAATTGTGTTGCACTCCCTGTTGCCGGAGt                >  2:2458300/1‑150 (MQ=255)
                                                                                                                                         tGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGCTAAATTGTGTTGCACTCCCTGTTGCCGGAGTTa              >  1:7888515/1‑150 (MQ=255)
                                                                                                                                         tGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGCTAAATTGTGTTGCACTCCCTGTTGCCGGAGTTa              <  2:2156651/150‑1 (MQ=255)
                                                                                                                                           aCATTGTTTCCGTTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATACACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGCTAAATTGTGTTGCACTCCCTGTTGCCGGAGTTACa            <  1:6206567/150‑1 (MQ=255)
                                                                                                                                           aCATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGCTAAATTGTGTTGCACTCCCTGTTGCCGGAGTTACa            <  2:3909640/150‑1 (MQ=255)
                                                                                                                                              ttGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGTATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGCTAAATTGTGTTGCACTCCCTGTTGCCGGAGTTACAGTg         <  1:745670/150‑1 (MQ=255)
                                                                                                                                                     cGGTTTTATTTTTTTCCACGACCAGGTGTGCACTCATCGGATAATACACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGCTAAATTGTGTTGCACTCCCTGTTGCCGGAGTTACAGTGCTCCAGa  <  1:3078084/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGCTAAATTGTGTTGCACTCCCTGTTGCCGGAGTTACAGTGCTCCAGA  >  NC_000913/151539‑151837

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