Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,778,124 A→C 5.4% intergenic (+342/+22) ypjF → / ← ypjA CP4‑57 prophage; toxin of the YpjF‑YfjZ toxin‑antitoxin system/adhesin‑like autotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,778,1240AC 5.4% 1097.8 / 43.3 331intergenic (+342/+22)ypjF/ypjACP4‑57 prophage; toxin of the YpjF‑YfjZ toxin‑antitoxin system/adhesin‑like autotransporter
Reads supporting (aligned to +/- strand):  ref base A (155/154);  new base C (8/10);  total (165/166)
Fisher's exact test for biased strand distribution p-value = 8.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

Only 100 of 331 total aligned reads displayed.

AGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATACCTGCCTGGTA  >  NC_000913/2777977‑2778271
                                                                                                                                                   |                                                                                                                                                    
  cccACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGAtgt                                                                                                                                                   <  3:2210556/149‑1 (MQ=255)
    caccaATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTgg                                                                                                                                                 <  3:1201544/149‑1 (MQ=255)
    caccaATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTgg                                                                                                                                                 >  2:1002153/1‑149 (MQ=255)
    caccaATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTgg                                                                                                                                                 <  3:338958/149‑1 (MQ=255)
     accaATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTggg                                                                                                                                                >  4:499003/1‑149 (MQ=255)
        aaTCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACATCTGTGGGCtt                                                                                                                                             <  2:175864/149‑1 (MQ=255)
        aaTCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACATCTGTGGGCtt                                                                                                                                             <  1:1672802/149‑1 (MQ=255)
         aTCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCttt                                                                                                                                            >  1:772459/1‑149 (MQ=255)
           cATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCttt                                                                                                                                            <  2:274494/147‑1 (MQ=255)
            aTGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACATCTGTGGGCtttcat                                                                                                                                         <  2:1654156/149‑4 (MQ=255)
            aTGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTc                                                                                                                                         >  3:587470/1‑149 (MQ=255)
            aTGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTc                                                                                                                                         <  3:997719/149‑1 (MQ=255)
              gATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGt                                                                                                                                       <  1:409937/149‑1 (MQ=255)
                   gACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATg                                                                                                                                  >  4:1657253/1‑149 (MQ=255)
                   gACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATg                                                                                                                                  <  2:1436040/149‑1 (MQ=255)
                   gACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATg                                                                                                                                  <  4:2225851/149‑1 (MQ=255)
                      ggTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATAAGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTc                                                                                                                               <  3:2108925/149‑1 (MQ=255)
                            aGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACg                                                                                                                         >  3:46927/1‑149 (MQ=255)
                            aGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACg                                                                                                                         >  1:1141552/1‑149 (MQ=255)
                              gACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGAc                                                                                                                       <  2:405118/149‑1 (MQ=255)
                                caacaCCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCa                                                                                                                     <  1:1002153/149‑1 (MQ=255)
                                 aacaCCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑g                                                                                                                   >  3:800501/1‑149 (MQ=255)
                                 aacaCCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑g                                                                                                                   >  1:2243830/1‑149 (MQ=255)
                                  acacCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑Gt                                                                                                                  <  4:663778/149‑1 (MQ=255)
                                       aacaaAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑Gtt                                                                                                                 <  4:2040685/145‑1 (MQ=255)
                                       aacaaAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTcaca                                                                                                             >  1:308888/1‑149 (MQ=255)
                                       aacaaAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTcaca                                                                                                             >  4:1700199/1‑149 (MQ=255)
                                       aacaaAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTcaca                                                                                                             >  3:2442214/1‑149 (MQ=255)
                                       aacaaAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTcaca                                                                                                             <  1:1901413/149‑1 (MQ=255)
                                        acaaAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCacac                                                                                                            <  2:1953365/149‑1 (MQ=255)
                                          aaaaaCAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCa                                                                                                          <  3:1642109/149‑1 (MQ=255)
                                           aaaaCAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAg                                                                                                         >  1:2474921/1‑149 (MQ=255)
                                                      ttAGAAGTCTCAGCAAAACCCCGACCAGACGGTGAGGTGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCGTTTTCGTTGATGGTCAGAACGACCA‑GTTCACATCAGCTACCGCGTTc                                                                                              <  2:1671794/149‑1 (MQ=255)
                                                      ttAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTc                                                                                              >  3:2058460/1‑149 (MQ=255)
                                                      ttAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTc                                                                                              >  3:2010840/1‑149 (MQ=255)
                                                      ttAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTc                                                                                              >  3:1785020/1‑149 (MQ=255)
                                                      ttAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGAGACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTc                                                                                              <  1:596611/149‑1 (MQ=255)
                                                       tAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTcc                                                                                             >  2:1110589/1‑149 (MQ=255)
                                                           aGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACgg                                                                                         <  4:1867692/149‑1 (MQ=255)
                                                           aGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACgg                                                                                         <  2:576032/149‑1 (MQ=255)
                                                           aGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACAGGAAAGCCCACAGATGTGGGCTTTTTCGTGGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACgg                                                                                         <  3:1303292/149‑1 (MQ=255)
                                                            gTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCGGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACggg                                                                                        >  2:2059611/1‑149 (MQ=255)
                                                            gTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACggg                                                                                        <  4:587470/149‑1 (MQ=255)
                                                             tctcAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACgggg                                                                                       >  1:2132661/1‑149 (MQ=255)
                                                               tcAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGTATAAGCAAAGGAGCCGTGGCTCATGGTAACATGAAAGCACACAGAGATGGGCTTTTTCGTTGATGGTCAGAACGACCG‑GTTCACACCAGCTACCACGTTCGACGGGGAt                                                                                     >  3:1532875/1‑149 (MQ=37)
                                                               tcAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGtt                                                                                               >  4:849729/1‑139 (MQ=255)
                                                               tcAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGAt                                                                                     >  2:200159/1‑149 (MQ=255)
                                                                    aaaaCACCGACCAGCCGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGT‑AGATGGAAGCCCCCAAGTTTTGGCCTTTTCCTTTGTGGTTAAAAAGATCAGGTTCACA‑CAGCCTCCGCCTTTCACGGGGAATcccaa                                                                               >  2:2381799/1‑146 (MQ=255)
                                                                    aaaaCACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTTCTGGTAACATGAAAGCCCACAGCTGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCcac                                                                                >  3:644618/1‑149 (MQ=255)
                                                                        cacCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCg                                                                            >  1:856322/1‑149 (MQ=255)
                                                                                aGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGTGGATTCCACACCGGCACCATg                                                                    >  4:335833/1‑149 (MQ=255)
                                                                                 gACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATgg                                                                   <  4:295981/149‑1 (MQ=255)
                                                                                  aCGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGc                                                                  <  3:1700199/149‑1 (MQ=255)
                                                                                    ggTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCta                                                                >  1:195349/1‑149 (MQ=255)
                                                                                       gaggagACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATac                                                             <  4:1988450/149‑1 (MQ=255)
                                                                                       gaggagACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATac                                                             >  4:2408332/1‑149 (MQ=255)
                                                                                              cATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATaccccaccc                                                      <  1:2023749/149‑1 (MQ=255)
                                                                                                tAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGa                                                                                                                        <  3:695359/82‑1 (MQ=255)
                                                                                                 aaaaGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCaa                                                   >  4:1402046/1‑149 (MQ=255)
                                                                                                  aaaGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAg                                                  <  4:800501/149‑1 (MQ=255)
                                                                                                   aaGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAgc                                                 <  3:849239/149‑2 (MQ=255)
                                                                                                   aaGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGa                                                 <  2:2335041/149‑1 (MQ=255)
                                                                                                   aaGGATACGCAAAGCAGCCGGGGCTCATGGTAACATGAAAGCCCACAGATGTGCGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTACACGGGGATACCACACCGGCACTATGGCTATACCCCACCCCaata                                                 <  4:1303749/149‑3 (MQ=21)
                                                                                                    aGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGATGATGGTCAGAACGACCA‑GTTCATACCAGCTACCGCGTTCCACGGGGATTCCACGCCGGCACTCTGGCTATACCCCACCCCAAGAt                                                <  2:2199236/149‑1 (MQ=255)
                                                                                                     ggATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGAt                                                                                                                                                     <  3:1404158/48‑1 (MQ=255)
                                                                                                     ggATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTACACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATg                                               >  1:1024759/1‑149 (MQ=255)
                                                                                                      gATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGc                                              >  4:1227268/1‑149 (MQ=255)
                                                                                                          cgaaaaGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAATCTGCCCGc                                          <  2:227316/146‑1 (MQ=255)
                                                                                                           gCAAAGGAGCCGCGGCTCCTGGTAACATGTGAGCCCACAGTTGTGGGCTTTTTCGTTGATGGTCAGAACGACGA‑GTTCACACCAGCTACCGCGATCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCt                                         <  4:644618/149‑1 (MQ=21)
                                                                                                           gCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCt                                         >  4:2116713/1‑149 (MQ=255)
                                                                                                           gCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCt                                         >  1:1544091/1‑149 (MQ=255)
                                                                                                              caGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTc                                                                                   <  4:1874623/103‑1 (MQ=255)
                                                                                                                 gAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑Gtt                                                                                                                 >  3:1493255/1‑71 (MQ=255)
                                                                                                                  aGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTa                                  <  3:2301396/149‑1 (MQ=255)
                                                                                                                  aGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTa                                  <  1:1778313/149‑1 (MQ=255)
                                                                                                                   gCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTAc                                 >  1:1233002/1‑149 (MQ=255)
                                                                                                                    ccGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCacac                                                                              <  4:1326048/103‑1 (MQ=255)
                                                                                                                    ccGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCacac                                                                              >  3:1326048/1‑103 (MQ=255)
                                                                                                                    ccGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTAct                                >  2:713521/1‑149 (MQ=255)
                                                                                                                          cTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCAAAGATGCCCGCTTAACGTACTGCTGaa                          >  1:1492385/1‑149 (MQ=255)
                                                                                                                          cTCCTGGTAACATGAAAGCCCACAGATGTCGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGGACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGaa                          <  2:920034/149‑1 (MQ=255)
                                                                                                                              tGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTGAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAg                      >  4:1527660/1‑149 (MQ=255)
                                                                                                                              tGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCtgagtgag                      <  1:536066/149‑1 (MQ=255)
                                                                                                                              tGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAg                      <  1:2353432/149‑1 (MQ=255)
                                                                                                                               ggTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAgg                     <  3:1035740/149‑1 (MQ=255)
                                                                                                                               ggTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAgg                     <  1:939653/149‑1 (MQ=255)
                                                                                                                               ggTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAgg                     >  1:1398200/1‑149 (MQ=255)
                                                                                                                               ggTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAgg                     <  2:856322/149‑1 (MQ=255)
                                                                                                                                 tAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCt                   >  4:377629/1‑149 (MQ=255)
                                                                                                                                  aaCATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCtt                  <  2:782171/149‑1 (MQ=255)
                                                                                                                                  aaCATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCtt                  <  1:2314808/149‑1 (MQ=255)
                                                                                                                                  aaCATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCtt                  <  2:1331247/149‑1 (MQ=255)
                                                                                                                                   aCATGAAAGCACACAGATGTGGGCTTTTTCGTGGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTCCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCTCTTAACGTACGGCTGAATGATGCttt                 <  4:1581246/149‑1 (MQ=18)
                                                                                                                                      tGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAAt              >  3:2069004/1‑149 (MQ=255)
                                                                                                                                        aaaGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATAc            <  2:527191/149‑1 (MQ=255)
                                                                                                                                            cccACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATAcctgc        <  2:1181383/149‑1 (MQ=255)
                                                                                                                                               acaGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATAcctgcctg     >  3:1939350/1‑149 (MQ=255)
                                                                                                                                               acaGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATAcctgcctg     >  3:285187/1‑149 (MQ=255)
                                                                                                                                                taGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATACCTGCCTgg    <  4:1706733/148‑1 (MQ=255)
                                                                                                                                                 aGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATACCTGCCTGGt   >  1:1723554/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCA‑GTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATACCTGCCTGGTA  >  NC_000913/2777977‑2778271

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑