Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 151,519 A→C 15.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,5190AC15.6% 340.4 / 54.9 139I27M (ATT→ATGyadKputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base A (57/60);  new base C (10/12);  total (67/72)
Fisher's exact test for biased strand distribution p-value = 8.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

Only 100 of 143 total aligned reads displayed.

CACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGT  >  NC_000913/151420‑151618
                                                                                                   |                                                                                                   
cACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCa                                                                                                     <  2:1862461/100‑1 (MQ=255)
 aCCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCaa                                                                                                    >  1:1020819/1‑100 (MQ=255)
 aCCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCaa                                                                                                    >  1:905888/1‑100 (MQ=255)
  ccccAACGACGCCAAGATCAGTATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAAt                                                                                                   <  1:2462468/100‑1 (MQ=255)
  ccccAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAAt                                                                                                   >  1:1381113/1‑100 (MQ=255)
    ccAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATgg                                                                                                 >  2:429995/1‑100 (MQ=255)
     cAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGc                                                                                                >  2:1173344/1‑100 (MQ=255)
          acgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGAc                                                                                           <  1:1179740/100‑1 (MQ=255)
          acgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGAc                                                                                           <  1:1319877/100‑1 (MQ=255)
           cgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACa                                                                                          <  2:2290648/100‑1 (MQ=255)
           cgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACa                                                                                          <  2:364361/100‑1 (MQ=255)
           cgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACa                                                                                          >  2:2369205/1‑100 (MQ=255)
           cgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACa                                                                                          >  1:635414/1‑100 (MQ=255)
           cgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTAGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACa                                                                                          <  2:1586748/100‑1 (MQ=255)
            gCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACAg                                                                                         <  1:710562/100‑1 (MQ=255)
            gCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACAg                                                                                         <  1:667871/100‑1 (MQ=255)
            gCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACAg                                                                                         <  1:576501/100‑1 (MQ=255)
             ccAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACAgg                                                                                        <  2:294937/100‑1 (MQ=255)
             ccAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAgg                                                                                        <  2:2040882/100‑1 (MQ=255)
              cAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACAgga                                                                                       <  1:2049761/100‑2 (MQ=255)
              cAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGc                                                                                       >  1:2108744/1‑100 (MQ=255)
              cAAGATCAATATTACCCTGGGTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACAgga                                                                                       <  1:789773/100‑2 (MQ=255)
               aaGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGc                                                                                       <  2:136799/99‑1 (MQ=255)
               aaGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCg                                                                                      <  1:1697192/100‑1 (MQ=255)
               aaGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCg                                                                                      <  1:1222874/100‑1 (MQ=255)
               aaGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCg                                                                                      <  1:2052813/100‑1 (MQ=255)
                 gATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGAt                                                                                    >  2:1249830/1‑100 (MQ=255)
                  aTCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATg                                                                                   <  2:2297381/100‑1 (MQ=255)
                   tCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGc                                                                                  <  2:2486018/100‑1 (MQ=255)
                   tCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATAGCCTGTCCAGCCATTGCTGGACAggatattg                                                                                  <  1:1574850/100‑7 (MQ=37)
                    cAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAATAAGATCAATATCCTGTCCAGCAATGGCTGGa                                                                                            >  2:1248197/1‑89 (MQ=255)
                      atatTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCaaa                                                                               >  1:99091/1‑100 (MQ=255)
                       tatTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCaaaa                                                                              >  2:348324/1‑100 (MQ=255)
                        atTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAAc                                                                             <  1:1173352/100‑1 (MQ=255)
                          tACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAg                                                                           <  2:837742/100‑1 (MQ=255)
                          tACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAAGGGCTGGACAGGCGATGCAAAACAg                                                                           <  1:2184996/100‑1 (MQ=255)
                           aCCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGt                                                                          >  2:1383611/1‑100 (MQ=255)
                           aCCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTGTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGt                                                                          <  2:2311304/100‑1 (MQ=255)
                            cccTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTa                                                                         <  1:2365338/100‑1 (MQ=255)
                                 gTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGac                                                                    <  2:1213931/100‑1 (MQ=255)
                                 gTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGac                                                                    >  1:2133715/1‑100 (MQ=255)
                                   tACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGaca                                                                   >  2:314479/1‑99 (MQ=255)
                                       gATTCCGCTTTTGGAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTGCAGCAATGGCTGGACAGGCGATGCAGAAAAGTAATGacagaaac                                                              >  1:414649/1‑100 (MQ=255)
                                       gATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaac                                                              >  2:1001153/1‑100 (MQ=255)
                                         ttCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaacag                                                            >  2:153072/1‑100 (MQ=255)
                                           ccGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaacagaa                                                          >  2:2364787/1‑100 (MQ=255)
                                             gCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAatt                                                        >  1:731538/1‑100 (MQ=255)
                                             gCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAatt                                                        <  2:850647/100‑1 (MQ=255)
                                              cTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAATAGAatta                                                       >  1:2213743/1‑100 (MQ=255)
                                              cTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAatta                                                       <  2:215616/100‑1 (MQ=255)
                                                  tGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATtct                                                   >  1:1903215/1‑100 (MQ=255)
                                                      ggTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtc                                               <  2:1806988/100‑1 (MQ=255)
                                                      ggTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtc                                               <  2:972712/100‑1 (MQ=255)
                                                      ggTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtc                                               <  2:594088/100‑1 (MQ=255)
                                                      ggTGCTGATTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtc                                               >  1:1347521/1‑100 (MQ=255)
                                                       ttccgatcttttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtca                                              <  1:560871/92‑1 (MQ=255)
                                                         gCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtcatc                                            <  2:731531/100‑1 (MQ=255)
                                                           tGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAg                                          >  1:954016/1‑100 (MQ=255)
                                                            gTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGTACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTAGTCTCTtcatca                                           >  2:1657260/1‑98 (MQ=255)
                                                              tttttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCtt                                       <  2:1809345/100‑1 (MQ=255)
                                                              tttttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCtt                                       >  1:521295/1‑100 (MQ=255)
                                                              tttttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCtt                                       >  2:1021463/1‑100 (MQ=255)
                                                                tttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTAc                                     >  2:1988613/1‑100 (MQ=255)
                                                                  tACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGt                                   <  1:880402/100‑1 (MQ=255)
                                                                  tACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGt                                   <  1:293798/100‑1 (MQ=255)
                                                                    cATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTg                                 <  1:580836/100‑1 (MQ=255)
                                                                     aTTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGa                                >  1:1571785/1‑100 (MQ=255)
                                                                        ggCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTg                             >  1:319892/1‑100 (MQ=255)
                                                                        ggCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTg                             <  2:1379545/100‑1 (MQ=255)
                                                                          caacaaGGTCGATATACTGTCCAGCCATGGCTGGACATGCGGTGCAAAACAATAATGACAGAAACAGAATTATTCTCTTCATCAGCGTACGTTGAGTgag                           >  2:1657771/1‑98 (MQ=255)
                                                                          caacaaGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTggg                           >  1:2329039/1‑100 (MQ=255)
                                                                            acaaGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGAt                         <  1:597280/100‑1 (MQ=255)
                                                                                ggTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCat                     >  2:2171895/1‑100 (MQ=255)
                                                                                ggTCAATATCCTCTCCAGCAATGGCTGGACAGGCGATGCAAACCAGTAATGACAGAAACAGCATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCat                     >  2:2242159/1‑100 (MQ=255)
                                                                                 gTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCata                    <  2:1907818/100‑1 (MQ=255)
                                                                                  tCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCatat                   <  2:1269395/100‑1 (MQ=255)
                                                                                   cAATATCCTGTCCAGCAATGGCTGGACAGGCTATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACATTGAGTAGGATCCat                     >  2:832691/1‑97 (MQ=255)
                                                                                   cAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCatata                  >  2:43470/1‑100 (MQ=255)
                                                                                      tatCCTGTCCAGCAATGGATGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAg               <  2:1592411/100‑1 (MQ=255)
                                                                                       atattagccAACAATAGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTAGAGTGGGATGCAtatata                <  2:876141/91‑1 (MQ=255)
                                                                                       atCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGt              >  1:1234653/1‑100 (MQ=255)
                                                                                       atCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGt              >  1:148073/1‑100 (MQ=255)
                                                                                        tCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGtt             <  1:2492581/100‑1 (MQ=255)
                                                                                         ccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTc            >  1:1112213/1‑100 (MQ=255)
                                                                                         ccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTc            >  2:555692/1‑100 (MQ=255)
                                                                                         ccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTc            >  1:1512818/1‑100 (MQ=255)
                                                                                         ccTGTCCAGCAATGGCTGAACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTc            >  1:938644/1‑100 (MQ=255)
                                                                                          cTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTcc           >  2:971965/1‑100 (MQ=255)
                                                                                          cTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTcc           <  2:52914/100‑1 (MQ=255)
                                                                                           tGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCt          >  1:433921/1‑100 (MQ=255)
                                                                                           tGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCt          <  1:1837072/100‑1 (MQ=255)
                                                                                           tGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCt          <  2:2056989/100‑1 (MQ=255)
                                                                                            gTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTg         <  2:694121/100‑1 (MQ=255)
                                                                                            gTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTg         >  1:782010/1‑100 (MQ=255)
                                                                                              ccAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGaa       >  2:2193920/1‑100 (MQ=255)
                                                                                              ccAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGaa       <  2:1664069/100‑1 (MQ=255)
                                                                                              ccAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGaa       >  1:791526/1‑100 (MQ=255)
                                                                                               cAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAAt      >  1:65010/1‑100 (MQ=255)
                                                                                               cAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAAAAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAAt      >  1:1048749/1‑100 (MQ=255)
                                                                                                aGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATc     >  1:750624/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGT  >  NC_000913/151420‑151618

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