Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 151,519 A→C 17.3% 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,5190AC17.3% 493.9 / 90.3 202I27M (ATT→ATGyadKputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base A (77/89);  new base C (19/17);  total (96/106)
Fisher's exact test for biased strand distribution p-value = 5.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 205 total aligned reads displayed.

CACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGT  >  NC_000913/151420‑151618
                                                                                                   |                                                                                                   
cACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATACCGCTTTTGCAGGAGCTGTTTTTTACATTGCCAACAAGGTCAATATCCTGTCCAGCa                                                                                                     >  1:2091267/1‑100 (MQ=255)
  ccccAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAAt                                                                                                   >  1:1008936/1‑100 (MQ=255)
  ccccAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAAt                                                                                                   <  2:1676319/100‑1 (MQ=255)
   cccAACGACGCCAAGATCAATATTACCCTGGTAAATGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCa                                                                                                     <  1:958054/97‑1 (MQ=255)
      aaCGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCt                                                                                               >  1:3552586/1‑100 (MQ=255)
         gacgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGa                                                                                            <  2:40965/100‑1 (MQ=255)
         gacgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGa                                                                                            <  2:581620/100‑1 (MQ=255)
           cgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACa                                                                                          <  2:1298777/100‑1 (MQ=255)
           cgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACa                                                                                          <  2:2290771/100‑1 (MQ=255)
           cgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACa                                                                                          <  2:303644/100‑1 (MQ=255)
           cgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATAGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACa                                                                                          <  1:1087759/100‑1 (MQ=255)
            gCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACAg                                                                                         <  1:2754590/100‑1 (MQ=255)
              cAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTATCCAGCAATGGCTGGACAGGc                                                                                       >  2:2184995/1‑100 (MQ=255)
              cAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCGGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACAgga                                                                                       <  1:3470388/100‑2 (MQ=255)
                aGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACAggata                                                                                     <  2:1622318/100‑4 (MQ=39)
                 gATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGAt                                                                                    <  1:3973588/100‑1 (MQ=255)
                  aTCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATg                                                                                   <  2:4027177/100‑1 (MQ=255)
                  aTCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATg                                                                                   <  2:1181316/100‑1 (MQ=255)
                  aTCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATg                                                                                   <  2:3715090/100‑1 (MQ=255)
                  aTCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATg                                                                                   <  2:2131523/100‑1 (MQ=255)
                   tCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACAggat                                                                                      <  1:1105584/96‑3 (MQ=255)
                   tCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGc                                                                                  <  2:2049938/100‑1 (MQ=255)
                    cAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCa                                                                                 >  1:4065678/1‑100 (MQ=255)
                     aaTATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCaa                                                                                >  2:1743687/1‑100 (MQ=255)
                      atatTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCaaa                                                                               <  2:2034091/100‑1 (MQ=255)
                      atatTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCaaa                                                                               >  2:442100/1‑100 (MQ=255)
                        atTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAAc                                                                             <  1:3021980/100‑1 (MQ=255)
                        atTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAAc                                                                             <  2:4228089/100‑1 (MQ=255)
                          ttctcTGGTTACTGACTCCGCTTTTGCAGGTGTTGTTTTTTACTTTGGTAACAATTTCAATATCCTGCCCAGCAATGGCTGGACAGGCGATGCAAAACAg                                                                           <  1:6921/96‑1 (MQ=255)
                          tACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAg                                                                           >  2:3291014/1‑100 (MQ=255)
                          tACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAg                                                                           <  2:813987/100‑1 (MQ=255)
                            cccTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTa                                                                         >  1:1347205/1‑100 (MQ=255)
                            cccTCGTTACTGATTCCGCTTTTGCAGGTGATGTTTTTTACATTGGCAACAAAGTCAATATACTGTCCACCAATGGCTGGACAGGCGATTCaaa                                                                               >  1:847866/1‑94 (MQ=255)
                               tGGCTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATg                                                                      <  2:1887619/100‑1 (MQ=255)
                                ggTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGa                                                                     >  2:1817718/1‑100 (MQ=255)
                                 gTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGac                                                                    <  2:385608/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                                                                 >  1:41365/1‑100 (MQ=255)
                                      tGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaa                                                               >  2:3699272/1‑100 (MQ=255)
                                      tGATTCCGCTTTTGCAGGTGATGTTTTTTACATTGGCAACAAGGTCAAGATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATCacagaa                                                                >  1:4352382/1‑99 (MQ=255)
                                        aTTCCGCTTTTGCAGGGGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaaca                                                             <  1:3605402/100‑1 (MQ=255)
                                         ttCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaaca                                                             <  1:3819367/99‑1 (MQ=255)
                                          tCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaacaga                                                           >  2:3984261/1‑100 (MQ=255)
                                          tCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaacaga                                                           >  2:1233825/1‑100 (MQ=255)
                                          tCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaacaga                                                           <  2:279159/100‑1 (MQ=255)
                                            cGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAat                                                         >  2:4278091/1‑100 (MQ=255)
                                            cGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAat                                                         >  2:1704428/1‑100 (MQ=255)
                                              cTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAatta                                                       <  2:3387048/100‑1 (MQ=255)
                                                 ttGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATtc                                                    <  1:2474269/100‑1 (MQ=255)
                                                 ttGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATtc                                                    <  1:1811779/100‑1 (MQ=255)
                                                    cAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCGGTCAAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTctct                                                 <  1:2047787/100‑1 (MQ=255)
                                                      ggTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtc                                               >  2:2433158/1‑100 (MQ=255)
                                                         gCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtcatc                                            >  2:889540/1‑100 (MQ=255)
                                                         gCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtcatc                                            >  1:1392258/1‑100 (MQ=255)
                                                          cTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtcatca                                           <  1:3358456/100‑1 (MQ=255)
                                                            gTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGc                                         <  2:3584351/100‑1 (MQ=255)
                                                            gTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGc                                         <  2:2322313/100‑1 (MQ=255)
                                                              tttttACATTGGCAACAAGGTCCATATCCTGTCCAGCCATGGCTGGATAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCtt                                       <  1:2975937/100‑1 (MQ=255)
                                                                   aCATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGtt                                  <  1:1260828/100‑1 (MQ=255)
                                                                   aCATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGtt                                  <  2:2977939/100‑1 (MQ=255)
                                                                    cATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTg                                 >  2:2550788/1‑100 (MQ=255)
                                                                    cATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTg                                 <  2:2424585/100‑1 (MQ=255)
                                                                     aTTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGa                                >  2:263021/1‑100 (MQ=255)
                                                                     aTTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCCGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGa                                <  2:3084409/100‑1 (MQ=255)
                                                                       tGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAg                               >  1:1966618/1‑99 (MQ=255)
                                                                        ggCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTg                             <  2:1388423/100‑1 (MQ=255)
                                                                        ggCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTg                             <  2:1347210/100‑1 (MQ=255)
                                                                        ggCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTg                             >  1:256370/1‑100 (MQ=255)
                                                                          aaacaaGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTggg                           >  1:3034382/2‑100 (MQ=255)
                                                                           aacaaGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTTTCTTCATCAGCTTACGTTGAGTGGGa                          <  2:2682156/100‑1 (MQ=255)
                                                                               aGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCa                      >  2:2722187/1‑100 (MQ=255)
                                                                               aGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCa                      <  1:4052589/100‑1 (MQ=255)
                                                                                tttgcatcgccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCag                     >  1:2110377/10‑99 (MQ=37)
                                                                                ggTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCat                     >  1:4002412/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:1777558/1‑100 (MQ=255)
                                                                                    catcgccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACACAATTATTCTCTTCATCAGCTTACGTTGCGTGGGATGCatatat                 >  1:1032030/6‑100 (MQ=38)
                                                                                     atatCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCAtatata                <  1:2257466/100‑1 (MQ=255)
                                                                                        tCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGtt             >  2:1336422/1‑100 (MQ=255)
                                                                                        gccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGtt             >  2:3203091/2‑100 (MQ=255)
                                                                                        gccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGtt             >  2:31418/2‑100 (MQ=255)
                                                                                         ccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTACTGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTc            >  2:2098405/1‑100 (MQ=255)
                                                                                         ccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTc            >  1:524444/1‑100 (MQ=255)
                                                                                         ccTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTc            <  1:2222973/100‑1 (MQ=255)
                                                                                         ccTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTc            >  1:3946022/1‑100 (MQ=255)
                                                                                          cTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTcc           >  2:3945631/1‑100 (MQ=255)
                                                                                          cTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTcc           >  2:1029180/1‑100 (MQ=255)
                                                                                          cTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTcc           >  1:3201428/1‑100 (MQ=255)
                                                                                          cTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTcc           <  2:1112293/100‑1 (MQ=255)
                                                                                          atgaccAGCAAGTGCTGTAAAGACGATGCAAAACAGTAAAGACACAAAAAGAAATATTATCTCCATCAGATTACGATGAGCAGGATGAATAAATAgatac           >  1:1875558/5‑96 (MQ=255)
                                                                                            gTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCt          >  1:4253809/1‑99 (MQ=255)
                                                                                             tCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGa        <  2:4081082/100‑1 (MQ=255)
                                                                                             tCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGa        <  2:3336951/100‑1 (MQ=255)
                                                                                              ccAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGaa       <  2:3959237/100‑1 (MQ=255)
                                                                                              ccAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGaa       <  1:2722172/100‑1 (MQ=255)
                                                                                               cAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAAt      <  2:2814922/100‑1 (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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: