Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 569,161 A→C 15.8% D87A (GAC→GCC) ‡ ybcK → DLP12 prophage; putative recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913569,1610AC15.8% 303.4 / 49.8 129D87A (GAC→GCC) ybcKDLP12 prophage; putative recombinase
Reads supporting (aligned to +/- strand):  ref base A (59/49);  new base C (9/12);  total (68/61)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01

Only 100 of 130 total aligned reads displayed.

AAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTC  >  NC_000913/569062‑569258
                                                                                                   |                                                                                                 
aaGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGa                                                                                                   >  2:1808670/1‑100 (MQ=255)
 aGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGAc                                                                                                  <  1:293287/100‑1 (MQ=255)
  ggAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACa                                                                                                 <  2:1811900/100‑1 (MQ=255)
      aGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACt                                                                                             >  1:1602310/1‑100 (MQ=255)
         aTGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTc                                                                                          >  2:1226339/1‑100 (MQ=255)
          tGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCa                                                                                         <  1:2414873/100‑1 (MQ=255)
           gCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCaa                                                                                        <  1:826489/100‑1 (MQ=255)
            cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGGATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTTaac                                                                                       <  1:2585967/100‑2 (MQ=255)
            cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaac                                                                                       <  2:1252254/100‑2 (MQ=255)
            cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaac                                                                                       <  1:179059/100‑2 (MQ=255)
            cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaac                                                                                       <  2:1036303/100‑2 (MQ=255)
            cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAag                                                                                       <  2:468696/100‑1 (MQ=255)
             accaTCAGGAGCTTTTTCGGAATTTTTAGCTGCTATAGAGCATGGTTATATATTGCCCACAACTACATTGGTAGTTAAAAGTCTGGACACACTTTCaa                                                                                        >  1:2706901/4‑98 (MQ=255)
             acaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAAATACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaac                                                                                       <  1:2053541/99‑2 (MQ=255)
              caATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAagag                                                                                     >  2:764512/1‑100 (MQ=255)
                 tCAGGAGCTTTTTCGGAATTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTTaat                                                                                       <  1:2708725/94‑2 (MQ=255)
                  cAGGGGCTCTTTCGGGATTTTGAGAAGGTATGGTGAATGGTTATAGGTTGCCAGGAACTACATCGTGAGTGGAAAGCCTGTACGGGATTTCagctaa                                                                                    <  1:1003328/97‑6 (MQ=255)
                   aGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaactaacaat                                                                                <  1:2609515/100‑9 (MQ=38)
                    ggAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAg                                                                               >  1:7344/1‑100 (MQ=255)
                    ggAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAg                                                                               >  1:1711807/1‑100 (MQ=255)
                    ggAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAg                                                                               <  1:133433/100‑1 (MQ=255)
                    ggAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAg                                                                               <  2:2060503/100‑1 (MQ=255)
                          ttttCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTg                                                                         >  2:2695978/1‑100 (MQ=255)
                             tCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAg                                                                      <  1:1999940/100‑1 (MQ=255)
                              cGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGc                                                                     >  1:797755/1‑100 (MQ=255)
                              cGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGc                                                                     >  2:2714900/1‑100 (MQ=255)
                               ggAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCg                                                                    >  1:1607179/1‑100 (MQ=255)
                                 aaTTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGAt                                                                  <  1:2695962/100‑1 (MQ=255)
                                  aTTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGAtt                                                                 >  2:2313468/1‑100 (MQ=255)
                                  aTTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGAtt                                                                 >  2:1996800/1‑100 (MQ=255)
                                   tttttAGATGCTAGAGAGCATGGTTATATATGGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTg                                                                >  1:1543739/1‑100 (MQ=255)
                                   tttttAAATGCTATAGATCATTGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTg                                                                <  1:2854821/100‑1 (MQ=255)
                                     tttAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGaa                                                              >  2:2360864/1‑100 (MQ=255)
                                         gATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTc                                                          >  2:1710094/1‑100 (MQ=255)
                                         gATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTc                                                          <  1:889692/100‑1 (MQ=255)
                                         gATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTc                                                          >  2:1710105/1‑100 (MQ=255)
                                         gATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTc                                                          <  1:889657/100‑1 (MQ=255)
                                           tGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTg                                                        >  2:1040939/1‑100 (MQ=255)
                                           tGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTg                                                        <  1:1884535/100‑1 (MQ=255)
                                               ataGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAAt                                                    <  2:2239246/100‑1 (MQ=255)
                                               agagagCATGGTTATATATTGCCAGGAACTAAATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAAt                                                    <  1:1709770/98‑1 (MQ=255)
                                                taGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAAtt                                                   >  1:2701874/1‑100 (MQ=255)
                                                 agagCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTg                                                  <  1:1725340/100‑1 (MQ=255)
                                                  gagCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGa                                                 <  2:2069726/100‑1 (MQ=255)
                                                   agCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAt                                                <  2:1803175/100‑1 (MQ=255)
                                                    gCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAtt                                               <  2:912706/100‑1 (MQ=255)
                                                    gCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAtt                                               <  2:325515/100‑1 (MQ=255)
                                                     cATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAttt                                              <  1:775724/100‑1 (MQ=255)
                                                      aTGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAtttt                                             <  2:1024956/100‑1 (MQ=255)
                                                          ttATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAAt                                         <  2:1557564/100‑1 (MQ=255)
                                                           tatataTTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATc                                        >  2:1697885/1‑100 (MQ=255)
                                                           tatataTTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATc                                        <  1:93881/100‑1 (MQ=255)
                                                             tatatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCAc                                      >  2:1812698/1‑100 (MQ=255)
                                                               tatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACgg                                    >  2:314061/1‑100 (MQ=255)
                                                               tatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACgg                                    <  2:2482891/100‑1 (MQ=255)
                                                                 ttGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTa                                  >  1:844722/1‑100 (MQ=255)
                                                                   gCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTAtt                                >  2:2912400/1‑100 (MQ=255)
                                                                   gCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTAtt                                <  1:1098887/100‑1 (MQ=255)
                                                                     cAGGAACAACATTGTTAGGTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGAGAAGCGATTGAACGTCTGAAAGTGACTGTGAATAACGGt                                   <  1:877194/95‑1 (MQ=255)
                                                                       ggAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATg                            <  2:291580/100‑1 (MQ=255)
                                                                        gAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGt                           <  1:1996789/100‑1 (MQ=255)
                                                                         aaCTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGtt                          >  2:1729943/1‑100 (MQ=255)
                                                                         aaCTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGtt                          <  1:62951/100‑1 (MQ=255)
                                                                         aaCTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGt                           >  1:464920/1‑99 (MQ=255)
                                                                          ncTACATTGTTAGTTAAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGAGAAGCGATTGAACGTCTGAAATTGATTTTGAATAACGGTATTGATGTta                         >  1:2807823/2‑100 (MQ=255)
                                                                          aCTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTta                         >  1:1938432/1‑100 (MQ=255)
                                                                            tACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTtata                       <  1:2582438/100‑1 (MQ=255)
                                                                              cATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAAc                     >  2:1622139/1‑100 (MQ=255)
                                                                              cATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAAc                     <  1:1654610/100‑1 (MQ=255)
                                                                                ttGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAActc                   >  2:2233624/1‑100 (MQ=255)
                                                                                 tGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAActct                  >  1:148913/1‑100 (MQ=255)
                                                                                 tGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAActct                  >  2:2099177/1‑100 (MQ=255)
                                                                                  nttAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCtt                 >  2:691525/2‑100 (MQ=255)
                                                                                  gTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCtt                 >  1:2470972/1‑100 (MQ=255)
                                                                                   ttAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAGTCACGGTATTGATATTATAACTCtt                 >  1:2818762/1‑99 (MQ=255)
                                                                                    ctcttGAAAGTCTGTCCAGCCTTTCAATAGAAAAAGTCGGTGAAGCGATTGAAAGTATGTAATTGATTTTGAATGACGGTATTGATGTTATAACtttt                 >  1:1131616/4‑98 (MQ=255)
                                                                                    ctcttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTg               >  1:362506/4‑100 (MQ=39)
                                                                                    ctcttGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTg               <  2:590589/97‑1 (MQ=255)
                                                                                     tcttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAActct                  >  1:2120411/3‑96 (MQ=255)
                                                                                     aGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGc              >  2:1551138/1‑100 (MQ=255)
                                                                                      nttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGc              >  1:2216104/2‑99 (MQ=255)
                                                                                      cttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCg             >  1:2551557/2‑100 (MQ=255)
                                                                                      cttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCg             >  2:2422905/2‑100 (MQ=255)
                                                                                      cttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCg             >  2:285278/2‑100 (MQ=255)
                                                                                      cttGAAAGTCTGTCCAGACTTTAAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCg             >  2:417343/2‑100 (MQ=255)
                                                                                      cttGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCg             >  2:288117/2‑100 (MQ=255)
                                                                                        tGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGAc           <  1:2099164/100‑1 (MQ=255)
                                                                                        tGAAAGTATGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGAc           >  2:1526563/1‑100 (MQ=255)
                                                                                         gAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGAcn          <  1:292511/100‑2 (MQ=255)
                                                                                         gAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACa          >  2:2329069/1‑100 (MQ=255)
                                                                                         gAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACa          <  2:1240912/100‑1 (MQ=255)
                                                                                           aaGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACaan        <  1:1615163/100‑2 (MQ=255)
                                                                                            aGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGagatcg       >  1:1081183/1‑95 (MQ=255)
                                                                                            aGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATa       <  2:2578247/100‑1 (MQ=255)
                                                                                               cTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACan    <  1:691522/100‑2 (MQ=255)
                                                                                                tGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGt   >  1:1885678/1‑100 (MQ=255)
                                                                                                tGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGt   >  2:1878323/1‑100 (MQ=255)
                                                                                                 ggACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTc  >  2:275835/1‑100 (MQ=255)
                                                                                                 ggACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTc  >  1:556355/1‑100 (MQ=255)
                                                                                                 ggACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTc  >  1:561398/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
AAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTC  >  NC_000913/569062‑569258

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: