Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 569,161 A→C 16.0% 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,1610AC16.0% 483.0 / 85.3 202D87A (GAC→GCC) ybcKDLP12 prophage; putative recombinase
Reads supporting (aligned to +/- strand):  ref base A (87/82);  new base C (20/13);  total (107/95)
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.42e-01

Only 100 of 207 total aligned reads displayed.

AAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTA  >  NC_000913/569062‑569260
                                                                                                   |                                                                                                   
aaGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGa                                                                                                     <  1:344839/100‑1 (MQ=255)
 aGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGAc                                                                                                    <  1:2325563/100‑1 (MQ=255)
  ggAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACa                                                                                                   <  1:4169365/100‑1 (MQ=255)
  ggAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACa                                                                                                   <  2:455501/100‑1 (MQ=255)
   gAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAg                                                                                                  <  2:138172/100‑1 (MQ=255)
         aTGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTc                                                                                            <  2:3990940/100‑1 (MQ=255)
          tGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCa                                                                                           <  2:2019470/100‑1 (MQ=255)
          tGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCa                                                                                           <  2:727436/100‑1 (MQ=255)
          tGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCa                                                                                           >  1:1360505/1‑100 (MQ=255)
          tGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCa                                                                                           <  1:535583/100‑1 (MQ=255)
          tGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCa                                                                                           >  2:325010/1‑100 (MQ=255)
            cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaac                                                                                         <  2:958993/100‑2 (MQ=255)
            cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaac                                                                                         <  2:2132848/100‑2 (MQ=255)
            cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaac                                                                                         <  2:1617467/100‑2 (MQ=255)
            cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaac                                                                                         <  2:3678581/100‑2 (MQ=255)
              caATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAagag                                                                                       >  2:1016723/1‑100 (MQ=255)
               aaTCAGGAGCTTTTTCGGAATTTTTAGATGCTAGAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAgtga                                                                                      <  1:259001/100‑1 (MQ=255)
                 tCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGaaa                                                                                    <  2:3183080/100‑1 (MQ=255)
                 tCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGaaa                                                                                    <  1:3390844/100‑1 (MQ=255)
                  cAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGaaaa                                                                                   >  1:3675070/1‑100 (MQ=255)
                   aGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaactaacaat                                                                                  <  1:3061912/100‑9 (MQ=38)
                    ggAGCTTTTTCGGAATTTTTAGATGCTATAGAGAATGGTTATATATTGCCAGGAACTACATTGTTAGTTGGAAGTCTGTCCAGACTTTCAAcagatcgga                                                                                 >  2:2521605/1‑95 (MQ=255)
                    cgAGCTTTTTGGGAATTTTTAGAAGCAATAGAGCATGGTTATATATTGCCAGGAACTAGGTTGTTAGTGGAAGATCTGGAGAGACTTTCAAGAGaaaaa                                                                                  >  2:1030068/2‑99 (MQ=255)
                       gCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCg                                                                              <  1:135320/100‑1 (MQ=255)
                       gCTTTTTCGGAATTTTTAGATGCAATAGGGCATGGTTATATAGTGCCAGGAACTACATTGTTAGTTGAAAGCCTGGACAGACTTTCAAGAGCAAAAGTc                                                                               <  2:839033/99‑1 (MQ=255)
                          ttttCGGAATTTTTAGATGCTATAGAGCATGGTTATATATCGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTg                                                                           >  2:1926480/1‑100 (MQ=255)
                             tCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCgagatcg                                                                        >  2:211914/1‑94 (MQ=255)
                              cGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGc                                                                       <  2:3831522/100‑1 (MQ=255)
                                  aTTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGAtt                                                                   <  2:3888320/100‑1 (MQ=255)
                                   tttttAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTg                                                                  <  1:495547/100‑1 (MQ=255)
                                    ttttAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGa                                                                 >  2:2188296/1‑100 (MQ=255)
                                     tttAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGaa                                                                <  1:2188300/100‑1 (MQ=255)
                                     tttAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGaa                                                                >  1:2672189/1‑100 (MQ=255)
                                     tttAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGaa                                                                <  1:1536344/100‑1 (MQ=255)
                                       tAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACg                                                              >  2:3231202/1‑100 (MQ=255)
                                       tAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACg                                                              <  2:259889/100‑1 (MQ=255)
                                       tAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACg                                                              >  1:759749/1‑100 (MQ=255)
                                          aTGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCt                                                           >  1:2858478/1‑100 (MQ=255)
                                            gCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGa                                                         >  1:2990811/1‑100 (MQ=255)
                                             cTATAGAGGATGGTTATATATTGCAAGGAACTGCATTGTTAGGTGAAAGTCTGGACAGACTTGCAAGAGAAAAAGTCGGTGAAGCGATTTAACGTCTGaa                                                        <  2:2907032/100‑1 (MQ=255)
                                             cTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTGGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGaa                                                        >  1:2368766/1‑100 (MQ=255)
                                                taGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAAtt                                                     <  1:105857/100‑1 (MQ=255)
                                                 agagCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTg                                                    >  1:697409/1‑100 (MQ=255)
                                                  gagCATGGTTATATATTGCAAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAATAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAAttaa                                                   >  2:4098608/1‑98 (MQ=255)
                                                    gCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAtt                                                 <  1:543096/100‑1 (MQ=255)
                                                    gCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAt                                                  >  2:2848225/1‑99 (MQ=255)
                                                     cATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAttt                                                >  1:2796798/1‑100 (MQ=255)
                                                       tGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTg                                              >  1:2768895/1‑100 (MQ=255)
                                                        ggTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGa                                             >  1:254432/1‑100 (MQ=255)
                                                           tatataTTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACTGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAAt                                           <  2:2677194/99‑1 (MQ=255)
                                                             tatatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCAc                                        >  1:2490981/1‑100 (MQ=255)
                                                             tatatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCAc                                        >  1:257268/1‑100 (MQ=255)
                                                              atatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACg                                       <  1:938668/100‑1 (MQ=255)
                                                               tatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACgg                                      <  1:1772455/100‑1 (MQ=255)
                                                               tatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACgg                                      >  1:3841613/1‑100 (MQ=255)
                                                               tatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACgg                                      >  1:225647/1‑100 (MQ=255)
                                                                   gCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTAtt                                  >  2:409370/1‑100 (MQ=255)
                                                                    ccAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTg                                 >  2:2414014/1‑100 (MQ=255)
                                                                      aGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGAt                               >  2:2394041/1‑100 (MQ=255)
                                                                          aCTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTta                           <  1:3200929/100‑1 (MQ=255)
                                                                           cTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTtat                          >  1:3758903/1‑100 (MQ=255)
                                                                           cTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTtat                          >  2:3012601/1‑100 (MQ=255)
                                                                            tACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTtata                         <  2:3062614/100‑1 (MQ=255)
                                                                            tACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTtata                         <  2:1325934/100‑1 (MQ=255)
                                                                              cATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGATAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATAGATGTTATAAc                       >  2:1157974/1‑100 (MQ=255)
                                                                              cATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAAc                       <  1:2397357/100‑1 (MQ=255)
                                                                                ttGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAActc                     <  2:2929128/100‑1 (MQ=255)
                                                                                ttGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAActc                     <  2:2614129/100‑1 (MQ=255)
                                                                                 tGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAActct                    >  2:182393/1‑100 (MQ=255)
                                                                                 tGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAActct                    <  1:3165483/100‑1 (MQ=255)
                                                                                    tAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTg                 <  2:4133149/100‑1 (MQ=255)
                                                                                    ctcttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTg                 >  1:2017422/4‑100 (MQ=39)
                                                                                     tcttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGc                >  2:3443840/3‑100 (MQ=255)
                                                                                      gatGAAAGACTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGATATTGATGTTATAACTCTTTGCg               >  2:2319947/3‑100 (MQ=255)
                                                                                      gTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCg               >  1:1915586/1‑100 (MQ=255)
                                                                                      gTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCg               <  2:2368761/100‑1 (MQ=255)
                                                                                      cttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCg               >  1:3608482/2‑100 (MQ=255)
                                                                                      cttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCg               >  1:251513/2‑100 (MQ=255)
                                                                                      cttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCg               >  2:485542/2‑100 (MQ=255)
                                                                                      cttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCg               >  2:4137061/2‑100 (MQ=255)
                                                                                      cttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATAGATTTTGAATCACGGTATTGATGTAATAACTCTTTGCg               >  2:3546806/2‑100 (MQ=255)
                                                                                        tGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGAc             >  1:2302979/1‑100 (MQ=255)
                                                                                        tGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGAc             >  1:3880067/1‑100 (MQ=255)
                                                                                        tGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGAc             <  1:1657991/100‑1 (MQ=255)
                                                                                        tGAAAGTCTGGACAGAAAAAAAAAAAAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGAc             <  2:614440/100‑1 (MQ=255)
                                                                                         gAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACa            >  2:175813/1‑100 (MQ=255)
                                                                                          aaaGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACaa           <  2:2888957/100‑1 (MQ=255)
                                                                                           aaGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAAt          <  1:610191/100‑1 (MQ=255)
                                                                                             gTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATAc        >  1:711684/1‑100 (MQ=255)
                                                                                               cTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAg      <  1:1818366/100‑1 (MQ=255)
                                                                                               cTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAg      <  1:2197683/100‑1 (MQ=255)
                                                                                                tGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGt     >  1:2109900/1‑100 (MQ=255)
                                                                                                 ggACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAATTTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTc    >  1:3945397/1‑100 (MQ=255)
                                                                                                 ggACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTc    >  1:393327/1‑100 (MQ=255)
                                                                                                 ggACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTc    >  2:1562622/1‑100 (MQ=255)
                                                                                                 ggACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTc    >  1:1149874/1‑100 (MQ=255)
                                                                                                  gACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCt   <  1:2089898/100‑1 (MQ=255)
                                                                                                  gACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCt   <  1:2728810/100‑1 (MQ=255)
                                                                                                   aCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCGCGGTATTAATGTTTTAACTCTTTTCGACAATACAGTCt   >  2:127695/1‑99 (MQ=255)
                                                                                                   aCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCta  <  1:1293478/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
AAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTA  >  NC_000913/569062‑569260

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: