Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 569,161 A→C 18.2% 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,1610AC18.2% 294.9 / 68.8 132D87A (GAC→GCC) ybcKDLP12 prophage; putative recombinase
Reads supporting (aligned to +/- strand):  ref base A (56/52);  new base C (10/14);  total (66/66)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

Only 100 of 132 total aligned reads displayed.

AGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCT  >  NC_000913/569063‑569259
                                                                                                  |                                                                                                  
 ggAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACa                                                                                                  <  1:2559123/100‑1 (MQ=255)
  gAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTGAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAg                                                                                                 <  2:2041174/100‑1 (MQ=255)
   tccgatctcacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGa                                                                                                <  1:48636/92‑1 (MQ=255)
     tgCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACt                                                                                              <  1:2375339/99‑1 (MQ=255)
      gCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACta                                                                                             >  2:2375350/1‑99 (MQ=255)
        aTGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTc                                                                                           <  2:1176919/100‑1 (MQ=255)
        aTGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTc                                                                                           <  1:416120/100‑1 (MQ=255)
        aTGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTc                                                                                           <  2:2355051/100‑1 (MQ=255)
         tGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCa                                                                                          <  1:2865000/100‑1 (MQ=255)
         tGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCa                                                                                          <  2:182509/100‑1 (MQ=255)
         tGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCa                                                                                          >  1:287723/1‑100 (MQ=255)
         tGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCa                                                                                          >  1:569539/1‑100 (MQ=255)
         tGCACAATCAGGAGATTTTTCCGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCa                                                                                          <  2:1583747/100‑1 (MQ=255)
          gCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCaa                                                                                         <  2:2719004/100‑1 (MQ=255)
          gCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCaa                                                                                         <  2:608548/100‑1 (MQ=255)
          gCACAATCAGGAGCTTTGTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTATTTGAAAGTCTGTCCAGACTTTCaa                                                                                         <  1:77498/100‑1 (MQ=255)
           cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCatc                                                                                        <  1:506881/100‑3 (MQ=255)
           cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaac                                                                                        <  1:2356687/100‑2 (MQ=255)
           cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaac                                                                                        <  1:2360736/100‑2 (MQ=255)
           cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaac                                                                                        <  2:786211/100‑2 (MQ=255)
           cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaac                                                                                        <  2:1461490/100‑2 (MQ=255)
           cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGaagatcgga                                                                                        >  2:48637/1‑92 (MQ=255)
           cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAag                                                                                        >  2:181073/1‑100 (MQ=255)
           cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAag                                                                                        <  1:2497139/100‑1 (MQ=255)
           cacaATCAGGAGATTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAag                                                                                        >  2:1913360/1‑100 (MQ=255)
            acaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAaga                                                                                       >  1:1924924/1‑100 (MQ=255)
             caATCAGGAGCTTTTTCGGAATTTTTCGATGCTATAGAGGATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGCCAGACTTTCaac                                                                                        <  1:761546/98‑2 (MQ=255)
                  aGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGaaaaa                                                                                 >  2:1659925/1‑100 (MQ=255)
                   ggAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAg                                                                                >  2:2498317/1‑100 (MQ=255)
                   ggAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAg                                                                                >  2:641860/1‑100 (MQ=255)
                   ggAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAg                                                                                <  2:2476281/100‑1 (MQ=255)
                      gCTTTCTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCg                                                                             >  2:2661776/1‑100 (MQ=255)
                        tttttCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGt                                                                           <  2:276420/100‑1 (MQ=255)
                        tttttCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGt                                                                           <  1:2829103/100‑1 (MQ=255)
                         ttttCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGta                                                                          >  2:2829113/1‑99 (MQ=255)
                         ttttCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTg                                                                          >  1:2004208/1‑100 (MQ=255)
                            tCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAg                                                                       >  2:777086/1‑100 (MQ=255)
                              ggAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCg                                                                     >  1:243143/1‑100 (MQ=255)
                              ggAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTAAATTGTTAGTTGAAAGTCTGGAAAGACTTTCAAGAGAAAAAGTCGGTGAATCg                                                                     >  1:277848/1‑100 (MQ=255)
                               naaTTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGa                                                                    >  1:1807431/2‑100 (MQ=255)
                                aaTTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGAt                                                                   <  2:1267853/100‑1 (MQ=255)
                                 aTTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGAtt                                                                  >  2:1132961/1‑100 (MQ=255)
                                  tttttAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTg                                                                 <  1:2096358/100‑1 (MQ=255)
                                  tttttAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTg                                                                 <  2:814107/100‑1 (MQ=255)
                                   ttttAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGa                                                                <  1:2539894/100‑1 (MQ=255)
                                    tttAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGaa                                                               >  1:2463091/1‑100 (MQ=255)
                                      tAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACg                                                             <  2:860843/100‑1 (MQ=255)
                                       aGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGt                                                            >  2:369294/1‑100 (MQ=255)
                                            cTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGaa                                                       >  1:1063360/1‑100 (MQ=255)
                                              ataGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAAt                                                     >  1:424904/1‑100 (MQ=255)
                                              ataGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAAt                                                     >  1:1290355/1‑100 (MQ=255)
                                              ataGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAAt                                                     >  1:1893774/1‑100 (MQ=255)
                                                  agCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAt                                                 >  2:2181747/1‑100 (MQ=255)
                                                  agCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAt                                                 <  1:1168165/100‑1 (MQ=255)
                                                  agCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAt                                                 <  1:1132959/100‑1 (MQ=255)
                                                   gCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAtt                                                >  2:42414/1‑100 (MQ=255)
                                                   gCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAtt                                                <  1:2234734/100‑1 (MQ=255)
                                                      tGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAttttn                                             <  1:2373636/100‑2 (MQ=255)
                                                      tGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTg                                             <  2:287725/100‑1 (MQ=255)
                                                       ggTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGa                                            <  1:1024281/100‑1 (MQ=255)
                                                         ttATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAAt                                          <  2:2260841/100‑1 (MQ=255)
                                                         ttATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAAt                                          <  2:2259862/100‑1 (MQ=255)
                                                          tatataTTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATc                                         <  2:1323999/100‑1 (MQ=255)
                                                            tatatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCAc                                       >  2:1358232/1‑100 (MQ=255)
                                                            tatatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCAc                                       <  2:2077209/100‑1 (MQ=255)
                                                             atatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACg                                      >  2:250693/1‑100 (MQ=255)
                                                             atatTGCCAGGAACTACATTGTTAGTTGAAAGTCCGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAGTTTGAATCACg                                      <  1:1525989/100‑1 (MQ=255)
                                                               atTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGt                                    >  2:1216676/1‑100 (MQ=255)
                                                               atTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGt                                    >  1:1201884/1‑100 (MQ=255)
                                                               atTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGt                                    >  2:1843724/1‑100 (MQ=255)
                                                                ttGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGAATATAAAAAAAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTa                                   <  1:899769/100‑1 (MQ=255)
                                                                  gCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTAtt                                 >  1:2826729/1‑100 (MQ=255)
                                                                   ccAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTg                                <  1:181071/100‑1 (MQ=255)
                                                                    cAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCAGTGAAGCGATTGAACATCTTCGATTTATTTTTCATCCCGGTATTGa                               <  1:409289/100‑1 (MQ=255)
                                                                           tACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTtata                        <  2:770646/100‑1 (MQ=255)
                                                                           tACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTtata                        >  2:2776801/1‑100 (MQ=255)
                                                                             cATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAAc                      <  1:369292/100‑1 (MQ=255)
                                                                                 ttctcttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCtt                  >  1:335770/6‑100 (MQ=38)
                                                                                   tAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTg                >  2:392601/1‑100 (MQ=255)
                                                                                   ctcttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTg                >  2:631303/4‑100 (MQ=39)
                                                                                   ctcttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTg                >  2:2227962/4‑100 (MQ=39)
                                                                                    tcttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGc               >  1:1993138/3‑100 (MQ=255)
                                                                                     gTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATTTTATAACTCTTTGCg              >  2:2861763/1‑100 (MQ=255)
                                                                                     gTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCg              >  1:2866829/1‑100 (MQ=255)
                                                                                     gTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCg              >  2:2831000/1‑100 (MQ=255)
                                                                                       tGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGa             >  1:2621736/1‑99 (MQ=255)
                                                                                       tGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGAc            >  1:1148816/1‑100 (MQ=255)
                                                                                       tGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGAc            >  2:883407/1‑100 (MQ=255)
                                                                                       tGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGAc            <  1:250691/100‑1 (MQ=255)
                                                                                       tGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGAGGTTATAACTCTTTGCGAc            >  1:2784759/1‑100 (MQ=255)
                                                                                        gAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACa           >  2:154687/1‑100 (MQ=255)
                                                                                         aaaGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAACGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACaa          >  1:1788394/1‑100 (MQ=255)
                                                                                          aaGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAAt         <  1:1358230/100‑1 (MQ=255)
                                                                                            gTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGATGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTCTGCGACAATAc       <  1:1532114/100‑1 (MQ=255)
                                                                                            gTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATAc       <  1:2852747/100‑1 (MQ=255)
                                                                                               tGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGt    >  2:2453993/1‑100 (MQ=255)
                                                                                               tGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGt    <  2:1893782/100‑1 (MQ=255)
                                                                                               tGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGt    <  2:137128/100‑1 (MQ=255)
                                                                                                ggACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTc   >  1:1897404/1‑100 (MQ=255)
                                                                                                 gACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCt  >  2:2902388/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
AGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCT  >  NC_000913/569063‑569259

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