Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,778,124 A→C 5.3% 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.3% 639.3 / 17.0 184intergenic (+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 (80/92);  new base C (7/4);  total (87/97)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 184 total aligned reads displayed.

AGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATACCTGCCTGGTA  >  NC_000913/2777977‑2778271
                                                                                                                                                   |                                                                                                                                                   
aGCCCACCCATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAATGGAGCCGCGGCTCCTGGTAACATGAAAGCCCCCAGAtg                                                                                                                                                   >  1:188333/1‑150 (MQ=255)
   ccaccaATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACGGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACAGGAAAGCCCACAGATGTgg                                                                                                                                                >  2:2375770/1‑150 (MQ=255)
    caccaATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAATTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTggg                                                                                                                                               <  1:277608/150‑1 (MQ=255)
          tCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATAACAGCCCACATCTGTGGGCtttca                                                                                                                                         <  2:1710775/150‑3 (MQ=255)
            aTGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCg                                                                                                                                       >  2:411206/1‑150 (MQ=255)
            aTGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCg                                                                                                                                       >  2:2627131/1‑150 (MQ=255)
            aTGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCg                                                                                                                                       >  1:1749033/1‑150 (MQ=255)
            aTGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCg                                                                                                                                       >  1:1749024/1‑150 (MQ=255)
            aTGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCg                                                                                                                                       >  2:1216535/1‑150 (MQ=255)
            aTGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCg                                                                                                                                       >  2:1216534/1‑150 (MQ=255)
              gATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGtt                                                                                                                                     <  2:1119316/150‑1 (MQ=255)
                     cGGTGTAAGGACAACACCAACAAAACCAGGAATTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTc                                                                                                                              <  1:1563204/150‑1 (MQ=255)
                     cGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTc                                                                                                                              <  1:1563212/150‑1 (MQ=255)
                      ggTGTACGGACACCACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAACCACCGACCAGCCGGTGAGGAGACATAAAAGGATCCGCAACGGCGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCa                                                                                                                             <  1:1291801/150‑1 (MQ=255)
                      ggTGTAAGGACAACACCAACAACCACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCa                                                                                                                             <  2:711685/150‑1 (MQ=255)
                      ggTGTAAGGACAACACCAACAAAAACCGGACGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCa                                                                                                                             <  2:253565/150‑1 (MQ=255)
                         gtAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACCCCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGaa                                                                                                                          <  1:1483622/150‑1 (MQ=255)
                              gACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCCGAACGAcc                                                                                                                     >  2:1945685/1‑150 (MQ=255)
                                    acCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTca                                                                                                               >  1:484039/1‑150 (MQ=255)
                                       aacaaAAACAGTAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAATGGATTCGCAAAGGAGCCGCGGCTCCTGGTAACATGATAGCCCACAGTTGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCacac                                                                                                            <  2:1411176/150‑1 (MQ=255)
                                        acaaAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACAcc                                                                                                           >  1:1336259/1‑150 (MQ=255)
                                         caaAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCa                                                                                                          <  2:2616908/150‑1 (MQ=255)
                                          aaaaaCAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAg                                                                                                         >  1:2215242/1‑150 (MQ=255)
                                            aaaCAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTGTCGTTGATGGTCAGAACGACCAGTTCACACCAGCt                                                                                                       >  2:1312535/1‑150 (MQ=255)
                                                agcaaGTTAGAAGTCTCACCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACcg                                                                                                   <  1:1664630/147‑1 (MQ=255)
                                                    acttAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGAAACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACCGATGTGGGCTATTTCGTTGAGGGGCAGAACGACCAGTTCACACCAGCTACCGCGtt                                                                                               <  2:2386526/148‑1 (MQ=255)
                                                       tAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCa                                                                                            >  1:1630948/1‑150 (MQ=255)
                                                         gAAGTCTCAGCAAAACCCCGACCCGACGGTGAGGAGACATAACAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACg                                                                                          <  2:440517/150‑1 (MQ=255)
                                                            gTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACgggg                                                                                       <  2:1987245/150‑1 (MQ=255)
                                                               tcAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGAtt                                                                                    <  2:587313/150‑1 (MQ=255)
                                                                 aGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTcc                                                                                  >  1:1182902/1‑150 (MQ=255)
                                                                  gCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCca                                                                                 >  1:2868528/1‑150 (MQ=255)
                                                                   cAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGATCACACCAGCTACCGCGTTCCACGGGGATTCcac                                                                                >  2:2498249/1‑150 (MQ=255)
                                                                    aaaaCACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCcaca                                                                               <  2:629829/150‑1 (MQ=255)
                                                                    aaaaCACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCcaca                                                                               <  1:311699/150‑1 (MQ=255)
                                                                     aaaCACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCacac                                                                              >  2:152667/1‑150 (MQ=255)
                                                                       acacCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCg                                                                            <  1:1954545/150‑1 (MQ=255)
                                                                         acCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGc                                                                          >  2:1617548/1‑150 (MQ=255)
                                                                           cGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCAc                                                                        >  1:1000400/1‑150 (MQ=255)
                                                                                aGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATgg                                                                   >  2:2516956/1‑150 (MQ=255)
                                                                                aGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATgg                                                                   <  1:2613733/150‑1 (MQ=255)
                                                                                aGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATgg                                                                   <  1:1440855/150‑1 (MQ=255)
                                                                                 gACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGc                                                                  <  2:1354130/150‑1 (MQ=255)
                                                                                 gACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGc                                                                  <  1:2649487/150‑1 (MQ=255)
                                                                                 gACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGc                                                                  <  2:322746/150‑1 (MQ=255)
                                                                                  aCGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCt                                                                 <  1:1550675/150‑1 (MQ=255)
                                                                                      tGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATac                                                             <  1:2727497/150‑1 (MQ=255)
                                                                                      tGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATac                                                             >  2:242086/1‑150 (MQ=255)
                                                                                       gaggagACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATacc                                                            <  1:2125488/150‑1 (MQ=255)
                                                                                        aggagACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCAACAGATGTGGGCTTTTTAGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACCGGGAATCCACACCGGCACTATGGCTCTaccc                                                           >  1:318717/1‑150 (MQ=37)
                                                                                           agaCATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATaccccac                                                        <  1:1235657/150‑1 (MQ=255)
                                                                                           agaCATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATaccccac                                                        <  1:2211058/150‑1 (MQ=255)
                                                                                              cATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATaccccacccc                                                     <  1:2135263/150‑1 (MQ=255)
                                                                                                tAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCaa                                                   <  1:642082/150‑1 (MQ=255)
                                                                                                tAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCaa                                                   <  1:2658443/150‑1 (MQ=255)
                                                                                                tAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACAAGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCaa                                                   <  1:279881/150‑1 (MQ=255)
                                                                                                  aaaGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCACGa                                                 >  2:349040/1‑150 (MQ=255)
                                                                                                  aaaGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGa                                                 >  2:391139/1‑150 (MQ=255)
                                                                                                       aTACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGAAGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGccc                                            >  1:1512808/1‑150 (MQ=255)
                                                                                                          cGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCt                                         <  1:233452/150‑1 (MQ=255)
                                                                                                           gCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGACCGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCtt                                        >  2:148860/1‑150 (MQ=255)
                                                                                                           gCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCtt                                        >  1:431044/1‑150 (MQ=255)
                                                                                                           gCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTATCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCtt                                        >  1:2281498/1‑150 (MQ=255)
                                                                                                            cAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTa                                       <  1:438199/150‑1 (MQ=255)
                                                                                                              caGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAAc                                     >  2:77522/2‑150 (MQ=255)
                                                                                                              aaGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAAc                                     <  2:1033666/150‑1 (MQ=255)
                                                                                                              aaGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAAc                                     <  1:754255/150‑1 (MQ=255)
                                                                                                                 gAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTa                                  <  2:617213/150‑1 (MQ=255)
                                                                                                                 gAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCAATACCCCCCCCCAAGATGCCCGCTTAACGTa                                  >  1:1941272/1‑150 (MQ=255)
                                                                                                                  aGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTAc                                 <  1:2569411/150‑1 (MQ=255)
                                                                                                                   gCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTAct                                >  2:1714456/1‑150 (MQ=255)
                                                                                                                   gCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTAct                                >  1:982020/1‑150 (MQ=255)
                                                                                                                       cgGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTActgctg                            >  1:805748/1‑150 (MQ=255)
                                                                                                                        ggCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGa                           <  2:672347/150‑1 (MQ=255)
                                                                                                                         gCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGaa                          <  2:436607/150‑1 (MQ=255)
                                                                                                                         gCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGaa                          <  2:777952/150‑1 (MQ=255)
                                                                                                                         gCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGaa                          <  1:2508434/150‑1 (MQ=255)
                                                                                                                         gCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGGTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCGATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGaa                          >  2:1874154/1‑150 (MQ=255)
                                                                                                                         gCTCCGGGTAACATGAAAGCCCACAGATGGGGGCTCCTTCGATGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCACGATGCCCGCTTAACGTACTGCTGaa                          <  2:2019581/150‑1 (MQ=18)
                                                                                                                           tCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGGACTGCTGAATg                        >  2:2167067/1‑150 (MQ=255)
                                                                                                                            ccTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGa                       <  1:1160886/150‑1 (MQ=255)
                                                                                                                            ccTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGa                       <  2:439064/150‑1 (MQ=255)
                                                                                                                             cTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAg                      >  1:936475/1‑150 (MQ=255)
                                                                                                                              tGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAgg                     >  1:2300151/1‑150 (MQ=255)
                                                                                                                                 tAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCtt                  >  2:2485989/1‑150 (MQ=255)
                                                                                                                                 tAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCtt                  <  1:578552/150‑1 (MQ=255)
                                                                                                                                   aCATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACACCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTa                >  2:2859487/1‑150 (MQ=255)
                                                                                                                                    cATGAAATCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTaa               <  1:55643/150‑1 (MQ=255)
                                                                                                                                      tGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAAc                                     <  1:535862/126‑1 (MQ=255)
                                                                                                                                      tGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAAc                                     >  2:535862/1‑126 (MQ=255)
                                                                                                                                      tGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATa             >  1:1122669/1‑150 (MQ=255)
                                                                                                                                      tGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATa             >  2:2083915/1‑150 (MQ=255)
                                                                                                                                      tGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATa             >  1:2473951/1‑150 (MQ=255)
                                                                                                                                       gAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATAc            >  2:2797298/1‑150 (MQ=255)
                                                                                                                                        aaaGCCCACATCTGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATAcc           >  1:1334744/1‑150 (MQ=255)
                                                                                                                                        aaaGCCCACATCTGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATAcc           >  2:855389/1‑150 (MQ=255)
                                                                                                                                         aaGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATAcct          <  1:2557904/150‑1 (MQ=255)
                                                                                                                                             ccACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATAcctgcct      >  2:557390/1‑150 (MQ=255)
                                                                                                                                                 aGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATACCTGCCTGGta  <  1:2024596/150‑1 (MQ=255)
                                                                                                                                                 aGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATACCTGCCTGGta  <  1:1991609/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATACCTGCCTGGTA  >  NC_000913/2777977‑2778271

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