Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 151,519 A→C 5.6% I27M (ATT→ATG)  yadK ← putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913151,5190AC 5.6% 696.1 / 26.4 198I27M (ATT→ATG) yadKputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base A (86/100);  new base C (3/8);  total (89/109)
Fisher's exact test for biased strand distribution p-value = 3.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

Only 100 of 200 total aligned reads displayed.

AGCCAGCAGCCATGACACCGCTAAGGGCTCTAAAATCAATATGTAACTGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGC  >  NC_000913/151263‑151792
                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                 
    agcagcCATGACACCGCTAAGGGCTCTAAAATCAATATGTAACTGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGTTTCCGTTTTTGCAGGAGACGTTTTTTACATTGGGAACAAGGTCAATATTCTGTCCATCAAGGGCTGGAGCGGCGATGCaaa                                                                                                                                                                                                                                                             >  4:26689/1‑275 (MQ=255)
         ccATGACACCGCTAAGGGCTCTAAAATCAATATGTAACTGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGt                                                                                                                                                                                                                                                        >  2:188962/1‑275 (MQ=255)
                 cctcTAACGGCTCTAAAATCAATATGTAACTGATTTAGCACATCGCCAGTTCACTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGac                                                                                                                                                                                                                                                  <  1:73437/270‑1 (MQ=255)
                                     aaTATGTAACTGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACGCTTTTCCACCTGGGTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGGCGGCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCGGTCCAGCAATGGATGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtca                                                                                                                                                                                                                            <  2:4906/275‑1 (MQ=255)
                                                                 gTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCatata                                                                                                                                                                                                <  2:20785/275‑1 (MQ=255)
                                                                                        tACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTg                                                                                                                                                                         <  4:188962/275‑1 (MQ=255)
                                                                                             tGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTCCGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAg                                                                                                                                                                             >  4:220352/1‑266 (MQ=255)
                                                                                             tGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCAAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAACCAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgttagt                                                                                                                                                                    >  2:169090/1‑275 (MQ=255)
                                                                                                 aGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGac                                                                                                                                                                                                                                                  <  2:253090/193‑1 (MQ=255)
                                                                                                 aGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTtata                                                                                                                                                                >  2:160889/1‑275 (MQ=255)
                                                                                                     aaCTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCtt                                                                                                                                                                                                                              >  4:150529/1‑209 (MQ=255)
                                                                                                         ctTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtcatc                                                                                                                                                                                                                          <  2:150529/209‑1 (MQ=255)
                                                                                                               cACCTGGTTGATAACTTTCAGGAGTAACATGACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGCTTACTGATTCCGCTTTTGCAGGTGCTGCTTTTTACATTGGCAACAAGGTCAGTATCCTGTCCAGCAATg                                                                                                                                                                                                                                                                                <  8:1219781/149‑1 (MQ=255)
                                                                                                                aCCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATgg                                                                                                                                                                                                                                                                               >  5:1108452/1‑149 (MQ=255)
                                                                                                                  cTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCt                                                                                                                                                                                                                                                                             >  7:259297/1‑149 (MQ=255)
                                                                                                                   tGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTg                                                                                                                                                                                                                                                                            <  6:1415575/149‑1 (MQ=255)
                                                                                                                     gTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTCCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGa                                                                                                                                                                                                                                                                          <  7:985222/149‑1 (MQ=255)
                                                                                                                      ttGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGAc                                                                                                                                                                                                                                                                         <  8:755772/149‑1 (MQ=255)
                                                                                                                      ttGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGAc                                                                                                                                                                                                                                                                         <  7:1363874/149‑1 (MQ=255)
                                                                                                                      ttGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGAc                                                                                                                                                                                                                                                                         <  8:907259/149‑1 (MQ=255)
                                                                                                                      ttGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGAc                                                                                                                                                                                                                                                                         <  6:926356/149‑1 (MQ=255)
                                                                                                                       tGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACa                                                                                                                                                                                                                                                                        <  5:991839/149‑1 (MQ=255)
                                                                                                                       tGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCACAACGACGCAAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACa                                                                                                                                                                                                                                                                        <  5:670666/149‑1 (MQ=255)
                                                                                                                        gATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAg                                                                                                                                                                                                                                                                       <  8:340665/149‑1 (MQ=255)
                                                                                                                        gATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAg                                                                                                                                                                                                                                                                       >  8:1308526/1‑149 (MQ=255)
                                                                                                                              tttCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATg                                                                                                                                                                                                                                                                 >  8:405871/1‑149 (MQ=255)
                                                                                                                              tttCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATg                                                                                                                                                                                                                                                                 >  8:223046/1‑149 (MQ=255)
                                                                                                                              tttCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATg                                                                                                                                                                                                                                                                 <  5:352575/149‑1 (MQ=255)
                                                                                                                                      gTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAg                                                                                                                                                                                                                                                         >  6:909747/1‑149 (MQ=255)
                                                                                                                                      gTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAg                                                                                                                                                                                                                                                         >  6:132721/1‑149 (MQ=255)
                                                                                                                                         aCATTACCTGAAAAATATCCCTCCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTaa                                                                                                                                                                                                                                                      <  7:1195729/149‑1 (MQ=255)
                                                                                                                                          cATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAAt                                                                                                                                                                                                                                                     <  6:520641/149‑1 (MQ=255)
                                                                                                                                           aTTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGc                                                                                                                                                                                                                                                                <  6:1386801/137‑1 (MQ=255)
                                                                                                                                           aTTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtcatc                                                                                                                                                                                                                          <  4:172191/175‑1 (MQ=255)
                                                                                                                                              aCCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGaca                                                                                                                                                                                                                                                 <  6:1235051/149‑1 (MQ=255)
                                                                                                                                                 tGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGc                                                                                                                                                                                                                                                                     >  7:1452290/1‑126 (MQ=255)
                                                                                                                                                 tGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGc                                                                                                                                                                                                                                                                     <  8:1452290/126‑1 (MQ=255)
                                                                                                                                                  gAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaa                                                                                                                                                                                                                                             <  8:1071548/149‑1 (MQ=255)
                                                                                                                                                  gAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaa                                                                                                                                                                                                                                             >  7:458621/1‑149 (MQ=255)
                                                                                                                                                  gAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaa                                                                                                                                                                                                                                             <  6:756006/149‑1 (MQ=255)
                                                                                                                                                       atatCCCACCCCAACGACGCCAAGATTAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaacagaa                                                                                                                                                                                                                                        <  8:60552/149‑1 (MQ=255)
                                                                                                                                                       atatCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaacagaa                                                                                                                                                                                                                                        <  5:1075394/149‑1 (MQ=255)
                                                                                                                                                         atCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCTTGTCCAGCCATTGCTg                                                                                                                                                                                                                                                                            >  5:931650/1‑111 (MQ=255)
                                                                                                                                                         atCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCTTGTCCAGCCATTGCTg                                                                                                                                                                                                                                                                            <  6:931650/111‑1 (MQ=255)
                                                                                                                                                          tCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAACCAGAatta                                                                                                                                                                                                                                     >  7:1327265/1‑149 (MQ=255)
                                                                                                                                                            ccACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCaa                                                                                                                                                                                                                                                                                  <  5:1183861/102‑1 (MQ=255)
                                                                                                                                                            ccACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCaa                                                                                                                                                                                                                                                                                  >  6:1183861/1‑102 (MQ=255)
                                                                                                                                                            ccACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAattatt                                                                                                                                                                                                                                   >  6:277778/1‑149 (MQ=255)
                                                                                                                                                            ccACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAattatt                                                                                                                                                                                                                                   >  8:1444666/1‑149 (MQ=255)
                                                                                                                                                            ccACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAattatt                                                                                                                                                                                                                                   >  5:1403412/1‑149 (MQ=255)
                                                                                                                                                              aCCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATtct                                                                                                                                                                                                                                 >  6:839121/1‑149 (MQ=255)
                                                                                                                                                                   aaCGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtca                                                                                                                                                                                                                            <  5:132721/149‑1 (MQ=255)
                                                                                                                                                                    acgacgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtcat                                                                                                                                                                                                                           <  5:909747/149‑1 (MQ=255)
                                                                                                                                                                      gacgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtcatca                                                                                                                                                                                                                         >  6:617058/1‑149 (MQ=255)
                                                                                                                                                                            aaGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTAc                                                                                                                                                                                                                   <  7:1074659/149‑1 (MQ=255)
                                                                                                                                                                                   atatTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGt                                                                                                                                                                                                            >  7:487066/1‑149 (MQ=255)
                                                                                                                                                                                        aCCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGAt                                                                                                                                                                                                       <  7:346375/149‑1 (MQ=255)
                                                                                                                                                                                         cccTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTc                                                                        >  2:319686/1‑275 (MQ=255)
                                                                                                                                                                                             ggTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCata                                                                                                                                                                                                  <  6:1462892/149‑1 (MQ=255)
                                                                                                                                                                                             ggTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGc                                                                                                                                              <  4:37631/201‑1 (MQ=255)
                                                                                                                                                                                              gTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCatat                                                                                                                                                                                                 >  7:753390/1‑149 (MQ=255)
                                                                                                                                                                                              gTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTNACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATTTGNGTGGTGCAGAAAACTTACCCGCAGTGACATTTTTTACGGTTTTTTNTNTTTCTACGACCCGGCgatcaccccttg                                                                    >  3:34918/1‑263 (MQ=255)
                                                                                                                                                                                               ttACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCatata                                                                                                                                                                                                <  7:333884/149‑1 (MQ=255)
                                                                                                                                                                                                tACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCatatat                                                                                                                                                                                               <  7:405871/149‑1 (MQ=255)
                                                                                                                                                                                                  cTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAg                                                                                                                                                                                             <  7:938920/149‑1 (MQ=255)
                                                                                                                                                                                                   taaTTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGc                                                                                                                                                                                                                                                                >  8:464403/3‑81 (MQ=255)
                                                                                                                                                                                                   tGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGt                                                                                                                                                                                            >  6:969124/1‑149 (MQ=255)
                                                                                                                                                                                                   tGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGt                                                                                                                                                                                            >  5:373892/1‑149 (MQ=255)
                                                                                                                                                                                                         cGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGa                                                                                                                                                                                      >  7:604561/1‑149 (MQ=255)
                                                                                                                                                                                                                 cAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTaa                                                                                                                                                                              >  6:1096505/1‑149 (MQ=255)
                                                                                                                                                                                                                 cAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTaa                                                                                                                                                                              <  5:797621/149‑1 (MQ=255)
                                                                                                                                                                                                                       cTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTacagattc                                                                                                                                                                        >  7:1307848/1‑142 (MQ=255)
                                                                                                                                                                                                                       cTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgt                                                                                                                                                                        <  5:580418/149‑1 (MQ=255)
                                                                                                                                                                                                                        tGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgtt                                                                                                                                                                       >  5:138170/1‑149 (MQ=255)
                                                                                                                                                                                                                          ttttttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAg                                                                                                                                                                                                                                                         <  7:28985/65‑1 (MQ=255)
                                                                                                                                                                                                                          ttttttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTg                                       <  4:319686/275‑1 (MQ=255)
                                                                                                                                                                                                                                aCATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAg                                                                                                                                                               <  6:393866/149‑1 (MQ=255)
                                                                                                                                                                                                                                aCATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAg                                                                                                                                                               <  5:636551/149‑1 (MQ=255)
                                                                                                                                                                                                                                aCATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAg                                                                                                                                                               <  6:1273737/149‑1 (MQ=255)
                                                                                                                                                                                                                                   ttGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGtt                                                                                                                                                                                           <  8:880579/118‑1 (MQ=255)
                                                                                                                                                                                                                                   ttGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGtt                                                                                                                                                            >  7:885061/1‑149 (MQ=255)
                                                                                                                                                                                                                                        aacaaGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTa                                                                                                                                                           <  6:726997/145‑1 (MQ=255)
                                                                                                                                                                                                                                        aacaaGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGt                                                                                                                                                       <  7:1253577/149‑1 (MQ=255)
                                                                                                                                                                                                                                                  atatCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCt                                                                                                                                                                          >  8:1014465/1‑120 (MQ=255)
                                                                                                                                                                                                                                                  atatCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCt                                                                                                                                             <  7:727926/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     tCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAAtgt                                                                                                                                                 >  8:277073/1‑142 (MQ=255)
                                                                                                                                                                                                                                                      ccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAACGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgt                                                                                                                                                                        >  7:781/1‑118 (MQ=255)
                                                                                                                                                                                                                                                      ccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAACGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgt                                                                                                                                                                        <  8:781/118‑1 (MQ=255)
                                                                                                                                                                                                                                                         gTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgtt                                                                                                                                                                       >  2:370652/1‑116 (MQ=255)
                                                                                                                                                                                                                                                          tCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGa                                                                                                                                     >  8:554513/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          tCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGa                                                                                                                                     <  6:578774/149‑1 (MQ=255)
                                                                                                                                                                                                                                                             aGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgt                                                                                                                                                                        <  5:25032/111‑1 (MQ=255)
                                                                                                                                                                                                                                                             aGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgt                                                                                                                                                                        >  6:25032/1‑111 (MQ=255)
                                                                                                                                                                                                                                                             aGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGaaaa                                                                                                                                  >  5:1253337/1‑149 (MQ=255)
                                                                                                                                                                                                                                                             aGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGaaaa                                                                                                                                  <  6:656169/149‑1 (MQ=255)
                                                                                                                                                                                                                                                             aGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGaaaa                                                                                                                                  >  5:791072/1‑149 (MQ=255)
                                                                                                                                                                                                                                                             aGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGaaaa                                                                                                                                  >  5:1472054/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              gCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAAc                                                                                                                                 >  8:171174/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              gCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAAc                                                                                                                                 >  7:138507/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               cAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGAAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGc  >  2:268769/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                 
AGCCAGCAGCCATGACACCGCTAAGGGCTCTAAAATCAATATGTAACTGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGCTGGTGTTACACAGC  >  NC_000913/151263‑151792

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