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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913151,5220GT17.3% 528.4 / 87.0 203A26A (GCC→GCAyadKputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base G (74/92);  new base T (20/15);  total (95/108)
Fisher's exact test for biased strand distribution p-value = 1.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 207 total aligned reads displayed.

CCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAG  >  NC_000913/151424‑151621
                                                                                                  |                                                                                                   
ccAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATgg                                                                                                    >  2:2047394/1‑100 (MQ=255)
 cAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGc                                                                                                   >  2:465853/1‑100 (MQ=255)
 cAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTGTTTACATTGGCAACAATGTCAATATCCTGTCCAGCAATGGc                                                                                                   >  1:2044585/1‑100 (MQ=255)
  aaCGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCt                                                                                                  <  1:3666675/100‑1 (MQ=255)
   acgacgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCtn                                                                                                 <  1:4033896/100‑2 (MQ=255)
     gacgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGa                                                                                               <  2:40965/100‑1 (MQ=255)
      acgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATCGCTGGAc                                                                                              <  2:2308293/100‑1 (MQ=255)
       cgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACa                                                                                             <  2:3959180/100‑1 (MQ=255)
       cgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACa                                                                                             <  2:1298777/100‑1 (MQ=255)
       cgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACa                                                                                             <  2:2397909/100‑1 (MQ=255)
       cgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACa                                                                                             <  2:303644/100‑1 (MQ=255)
       cgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACa                                                                                             <  2:1743227/100‑1 (MQ=255)
       cgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACa                                                                                             >  2:2086105/1‑100 (MQ=255)
       cgCCAAGATCAATATTACACTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACa                                                                                             >  1:224549/1‑100 (MQ=255)
        gCCAAGATCAATATTCCCCTGGTTACGTATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGCCAATATCCTGTCCAGCCATTGCt                                                                                                  <  1:1943713/94‑1 (MQ=255)
        gCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACAg                                                                                            <  2:2203637/100‑1 (MQ=255)
        gCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACAg                                                                                            <  1:2754590/100‑1 (MQ=255)
          cAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGc                                                                                          <  1:18202/100‑1 (MQ=255)
          cAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGc                                                                                          >  1:1163051/1‑100 (MQ=255)
          cAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCGGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACAgga                                                                                          <  1:3470388/100‑2 (MQ=255)
            aGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACAggata                                                                                        <  2:1622318/100‑4 (MQ=39)
             gATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAGGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGAt                                                                                       >  2:2542539/1‑100 (MQ=255)
             gATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGAt                                                                                       >  2:395979/1‑100 (MQ=255)
              aTCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATg                                                                                      <  2:1181316/100‑1 (MQ=255)
              aTCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATg                                                                                      <  2:4027177/100‑1 (MQ=255)
              aTCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATg                                                                                      <  2:3715090/100‑1 (MQ=255)
                    atTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAAc                                                                                <  1:3021980/100‑1 (MQ=255)
                      ttctcTGGTTACTGACTCCGCTTTTGCAGGTGTTGTTTTTTACTTTGGTAACAATTTCAATATCCTGCCCAGCAATGGCTGGACAGGCGATGCAAAACAg                                                                              <  1:6921/96‑1 (MQ=255)
                      tACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAg                                                                              <  2:813987/100‑1 (MQ=255)
                      tACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAg                                                                              >  2:3291014/1‑100 (MQ=255)
                      tACCCCGGTTACTGATTCCTCTTTTTCAGGTGCTGTTTTTTACATTGTCAACAAGGCCAATATCCCGTCCAGCAATGGCTCGACAGGCGATGCCAAAc                                                                                <  1:228880/98‑1 (MQ=255)
                        cccTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTa                                                                            >  1:1347205/1‑100 (MQ=255)
                        cccTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTa                                                                            >  2:2425695/1‑100 (MQ=255)
                        cccTCGTTACTGATTCCGCTTTTGCAGGTGATGTTTTTTACATTGGCAACAAAGTCAATATACTGTCCACCAATGGCTGGACAGGCGATTCaaa                                                                                  >  1:847866/1‑94 (MQ=255)
                          cTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAAt                                                                          <  1:4308648/100‑1 (MQ=255)
                             gTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGac                                                                       <  2:385608/100‑1 (MQ=255)
                               tACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacag                                                                     <  2:2770997/100‑1 (MQ=255)
                               tACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacag                                                                     <  2:2121989/100‑1 (MQ=255)
                               tACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacag                                                                     >  1:155935/1‑100 (MQ=255)
                               tACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacag                                                                     <  1:3303476/100‑1 (MQ=255)
                                aCTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacaga                                                                    >  1:142345/1‑100 (MQ=255)
                                aCTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacaga                                                                    <  2:2715048/100‑1 (MQ=255)
                                  tGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaa                                                                  >  2:3699272/1‑100 (MQ=255)
                                   gATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACACGGTCAATATCCTGTCCAGCAATAGCTCGACTGGCGATGCAAACCAGTCATGacagaaa                                                                  >  1:3556796/1‑99 (MQ=255)
                                   gATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaac                                                                 >  2:2133783/1‑100 (MQ=255)
                                    aTTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaaca                                                                <  2:2634116/100‑1 (MQ=255)
                                      tCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaacaga                                                              >  2:1233825/1‑100 (MQ=255)
                                          cTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAatta                                                          <  2:3387048/100‑1 (MQ=255)
                                           ttttGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAattat                                                         <  1:3984305/100‑1 (MQ=255)
                                              tGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATtct                                                      >  1:1272533/1‑100 (MQ=255)
                                                  ggTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtc                                                  >  2:1950103/1‑100 (MQ=255)
                                                    tGCTGTTTTTTACATTGGCAACAAGGTAAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAACACAGTAATGACAGAAACAGAATTATTCTCTtcat                                                <  1:4366778/100‑1 (MQ=255)
                                                     gCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtcatc                                               >  2:889540/1‑100 (MQ=255)
                                                     gCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtcatc                                               >  1:1392258/1‑100 (MQ=255)
                                                      cTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtcatca                                              <  1:3358456/100‑1 (MQ=255)
                                                      cTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtcatca                                              >  1:1094629/1‑100 (MQ=255)
                                                       tGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAg                                             >  2:710899/1‑100 (MQ=255)
                                                        gTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGc                                            <  2:3584351/100‑1 (MQ=255)
                                                          tttttACATTGGCAACAAGGTCCATATCCTGTCCAGCCATGGCTGGATAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCtt                                          <  1:2975937/100‑1 (MQ=255)
                                                               aCATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGtt                                     <  1:1029697/100‑1 (MQ=255)
                                                               aCATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGtt                                     <  2:2977939/100‑1 (MQ=255)
                                                                 aTTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGa                                   <  1:4064855/100‑1 (MQ=255)
                                                                 aTTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCCGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGa                                   <  2:3084409/100‑1 (MQ=255)
                                                                   tGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGt                                 <  1:4158949/100‑1 (MQ=255)
                                                                    ggCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTg                                <  2:1347210/100‑1 (MQ=255)
                                                                    ggCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTg                                <  2:1388423/100‑1 (MQ=255)
                                                                      caacaaGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTggg                              <  2:1679075/100‑1 (MQ=255)
                                                                      aaacaaGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTggg                              >  1:3034382/2‑100 (MQ=255)
                                                                           aGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCa                         <  1:2217797/100‑1 (MQ=255)
                                                                           aGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCa                         >  1:3860961/1‑100 (MQ=255)
                                                                             ggcAATATCCTGTGCAGCAATGGCTGGACAGACGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCata                       >  2:2963068/3‑100 (MQ=255)
                                                                             gTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCata                       <  2:859795/100‑1 (MQ=255)
                                                                             gTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCata                       <  2:1064962/100‑1 (MQ=255)
                                                                              tCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCatat                      >  1:1057607/1‑100 (MQ=255)
                                                                                cttccgatctCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCatatat                    <  1:1763786/91‑1 (MQ=255)
                                                                                catcgccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCatatat                    >  1:3959211/6‑100 (MQ=38)
                                                                                catcgccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCatatat                    >  2:432318/6‑100 (MQ=38)
                                                                                 atatCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCAtatata                   <  2:256371/100‑1 (MQ=255)
                                                                                    tCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGtt                >  1:2655972/1‑100 (MQ=255)
                                                                                    tCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGtt                >  2:1336422/1‑100 (MQ=255)
                                                                                    gccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGtt                >  2:3203091/2‑100 (MQ=255)
                                                                                    gccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGtt                >  2:31418/2‑100 (MQ=255)
                                                                                    gccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGtt                >  1:2069088/2‑100 (MQ=255)
                                                                                    gccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGtt                >  2:2794890/2‑100 (MQ=255)
                                                                                     ccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTc               >  2:4257217/1‑100 (MQ=255)
                                                                                     ccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTc               >  1:524444/1‑100 (MQ=255)
                                                                                     ccTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTc               >  1:2540944/1‑100 (MQ=255)
                                                                                     ccTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTc               <  1:1789737/100‑1 (MQ=255)
                                                                                      cTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTcc              <  2:1112293/100‑1 (MQ=255)
                                                                                      cTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTcc              >  2:1029180/1‑100 (MQ=255)
                                                                                      cTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTcc              >  1:3201428/1‑100 (MQ=255)
                                                                                         tCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGa           <  2:3336951/100‑1 (MQ=255)
                                                                                         tCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGa           <  2:4081082/100‑1 (MQ=255)
                                                                                          ccAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGaa          <  1:2722172/100‑1 (MQ=255)
                                                                                          ccAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGaa          >  1:4081023/1‑100 (MQ=255)
                                                                                          ccAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGATcggaag          >  2:1763792/1‑95 (MQ=255)
                                                                                               aaTGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGt     >  2:1847508/1‑100 (MQ=255)
                                                                                                aTGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTa    <  2:3271287/100‑1 (MQ=255)
                                                                                                aTGGCGGGACAGGCGATGGAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTGACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTa    <  2:3990200/100‑1 (MQ=255)
                                                                                                  ggCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAg  >  2:421682/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
CCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAG  >  NC_000913/151424‑151621

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: