Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,428,671 T→C 5.1% H16R (CAC→CGC)  ytfB ← cell division protein YtfB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,428,6710TC 5.1% 303.4 / 3.6 97H16R (CAC→CGC) ytfBcell division protein YtfB
Reads supporting (aligned to +/- strand):  ref base T (49/43);  new base C (3/2);  total (52/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGCGCTTCACGCGTCACCACCGGCGCGTTGGGCGTATCATCAGATGGAAGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATA  >  NC_000913/4428525‑4428820
                                                                                                                                                  |                                                                                                                                                     
aCTGCGCTTCACGCGTCACCACCGGCGCGTTGGGCGTATCATCAGATGGAAGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGcc                                                                                                                                                    >  1:635823/1‑150 (MQ=255)
 cTGCGCTTCACGCGTCACCACCGGCGCGTTGGGCGTATCATCAGATGGAAGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCa                                                                                                                                                   >  1:909874/1‑150 (MQ=255)
    cgcTTCACGCGTCACCACCGGCGCGTTGGGCGTATCATCAGATGGAAGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGAc                                                                                                                                                >  2:577743/1‑150 (MQ=255)
      cTTCACGCGTCACCACCGGCGCGTTGGGCGTATCATCAGATGGAAGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACtt                                                                                                                                              >  1:1667158/1‑150 (MQ=255)
      cTTCACGCGTCACCACCGGCGCGTTGGGCGTATCATCAGATGGAAGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACtt                                                                                                                                              >  1:1233184/1‑150 (MQ=255)
                  ccaccGGCGCGTTGGGCGTATCATCAGATGGAAGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGggttg                                                                                                                                  >  1:112274/1‑150 (MQ=255)
                       ggCGCGTTGGGCGTATCATCAGATGGAAGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGtttt                                                                                                                             <  1:111104/150‑1 (MQ=255)
                         cgcgTTGGGCGTATCATCAGATGGAAGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGCGCCAGACTTTCTCCAGGGTTGGTTTTAg                                                                                                                           >  1:274740/1‑150 (MQ=255)
                          gcgTTGGGCGTATCATCAGATGGAAGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGt                                                                                                                          <  2:491402/150‑1 (MQ=255)
                                ggCGTATCATCAGATGGAAGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGCCTTTCTCCAGGGTTGGTTTTAGTTCaaac                                                                                                                    >  1:1008708/1‑149 (MQ=255)
                                ggCGTATCATCAGATGGAAGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAg                                                                                                                    <  2:901742/150‑1 (MQ=255)
                                  cGTATCATCAGATGGAAGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAgcg                                                                                                                  >  2:414558/1‑150 (MQ=255)
                                    tATCATCAGATGGAAGCAGAAATCCGACCCCCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTGGGTTTTAGTTCAAAGCGcc                                                                                                                >  2:553088/1‑150 (MQ=255)
                                    tATCATCAGATGGAAGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGcc                                                                                                                >  1:1666822/1‑150 (MQ=255)
                                    tATCATCAGATGGAAGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGcc                                                                                                                >  1:941703/1‑150 (MQ=255)
                                    tATCATCAGATGGAAGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGcc                                                                                                                >  2:1173760/1‑150 (MQ=255)
                                     atcatcAGATGGAAGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGccc                                                                                                               >  2:58812/1‑150 (MQ=255)
                                      tcatcaGATGGAAGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTc                                                                                                                        >  1:190300/1‑140 (MQ=255)
                                      tcatcaGATGGAAGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTc                                                                                                                        <  2:190300/140‑1 (MQ=255)
                                        atcaGATGGAAGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCggg                                                                                                            <  1:577743/150‑1 (MQ=255)
                                         tcaGATGGAAGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGc                                                                                                           >  1:120468/1‑150 (MQ=255)
                                          caGATGGAAGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCa                                                                                                          >  2:240211/1‑150 (MQ=255)
                                            gATGGAAGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATg                                                                                                        <  1:916093/150‑1 (MQ=255)
                                                gAAGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGttt                                                                                                    <  2:1632456/150‑1 (MQ=255)
                                                gAAGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGttt                                                                                                    >  1:1308202/1‑150 (MQ=255)
                                                 aaGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTtc                                                                                                   >  1:1644064/1‑150 (MQ=255)
                                                 aaGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTtc                                                                                                   >  1:753522/1‑150 (MQ=255)
                                                  aGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGGGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTtcc                                                                                                  <  1:1022308/150‑1 (MQ=255)
                                                    cAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTtcctc                                                                                                >  2:1428840/1‑150 (MQ=255)
                                                    cAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTtcctc                                                                                                <  2:967861/150‑1 (MQ=255)
                                                      gAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTtcctc                                                                                                <  2:1557077/148‑1 (MQ=255)
                                                      gAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTtcctc                                                                                                >  1:1557077/1‑148 (MQ=255)
                                                      gAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTcctcct                                                                                              >  1:462434/1‑150 (MQ=255)
                                                        aaTCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGc                                                                                            <  1:617561/150‑1 (MQ=255)
                                                         atacGACCACCAATACAATGGCGGCAAAAATGATCCCCCGACGATGCATCTGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCTTGCCATACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCt                                                                                           <  1:1598062/147‑1 (MQ=255)
                                                            cGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCcg                                                                                        <  2:417933/150‑1 (MQ=255)
                                                            cGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCcg                                                                                        >  2:477610/1‑150 (MQ=255)
                                                             gACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCcgc                                                                                       >  2:594274/1‑150 (MQ=255)
                                                              accaccAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCGGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCcgcg                                                                                      <  1:1588686/150‑1 (MQ=255)
                                                              accaccAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCcgcg                                                                                      <  2:1233184/150‑1 (MQ=255)
                                                                caccaATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTc                                                                                    >  1:110377/1‑150 (MQ=255)
                                                                 accaATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCt                                                                                   >  1:5042/1‑150 (MQ=255)
                                                                 accaATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCt                                                                                   <  2:1308202/150‑1 (MQ=255)
                                                                       aCAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTcc                                                                             >  2:1256918/1‑150 (MQ=255)
                                                                       aCAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTcc                                                                             <  2:112274/150‑1 (MQ=255)
                                                                         aaTGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTc                                                                           <  2:635823/150‑1 (MQ=255)
                                                                             gcggcAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATc                                                                       >  1:1514062/1‑150 (MQ=255)
                                                                             gcggcAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATc                                                                       >  1:675593/1‑150 (MQ=255)
                                                                               ggcAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCAt                                                                     <  2:1667158/150‑1 (MQ=255)
                                                                                 cAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAg                                                                   >  1:805921/1‑150 (MQ=255)
                                                                                      atGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTAt                                                              >  1:269599/1‑150 (MQ=255)
                                                                                       tGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATg                                                             >  1:692353/1‑150 (MQ=255)
                                                                                         aTCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGGCTTGTTCCTCAATCATAGCCTATGaa                                                           <  2:1644064/150‑1 (MQ=255)
                                                                                           cccccGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATa                                                         <  1:1428840/150‑1 (MQ=255)
                                                                                            ccccGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATaa                                                        >  2:173300/1‑150 (MQ=255)
                                                                                              ccGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGc                                                      <  1:971448/150‑1 (MQ=255)
                                                                                              ccGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGc                                                      <  1:477610/150‑1 (MQ=255)
                                                                                                    aTGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTGTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGc                                                >  1:1486948/1‑150 (MQ=255)
                                                                                                        aTCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCAAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTatg                                            >  2:582786/1‑150 (MQ=255)
                                                                                                         tcgtgggcATCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTTGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTatga                                           <  2:307576/145‑1 (MQ=255)
                                                                                                         tCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTatga                                           >  1:1602998/1‑150 (MQ=255)
                                                                                                         tCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTatga                                           <  1:1256918/150‑1 (MQ=255)
                                                                                                         tCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTatga                                           <  2:5042/150‑1 (MQ=255)
                                                                                                           ggcggcAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTatgatg                                         >  2:288410/1‑150 (MQ=255)
                                                                                                            gcggcAGCGGGTCCATGAAGTGAAAATTGTCGGGCGCGTGGGAGACTTTCTCCTGGGTTGGTTTTTGTTCATAGCGCCCGGGCATGGTTTGCTCCTGGTGGGCGTCTTGTTCCTCAATCATAGGCTATGAATAAGCTAACGCTATGATGt                                        <  2:1598645/150‑1 (MQ=255)
                                                                                                             cggcAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTc                                       >  1:484047/1‑150 (MQ=255)
                                                                                                                 aGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTg                                   <  1:594274/150‑1 (MQ=255)
                                                                                                                  gCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTgg                                  <  1:1522751/150‑1 (MQ=255)
                                                                                                                  gCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTgg                                  <  1:1040172/150‑1 (MQ=255)
                                                                                                                  gCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTgg                                  <  1:240211/150‑1 (MQ=255)
                                                                                                                   cGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGt                                 >  1:1007297/1‑150 (MQ=255)
                                                                                                                   cGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGt                                 >  2:1632907/1‑150 (MQ=255)
                                                                                                                      gTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTaaa                              >  1:696965/1‑150 (MQ=255)
                                                                                                                         cATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAAccc                           >  2:217309/1‑150 (MQ=255)
                                                                                                                              aGCGAAAATTGTCCGGCGCGCGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCtt                      <  2:274740/150‑1 (MQ=255)
                                                                                                                                cGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTtta                    >  2:668905/1‑150 (MQ=255)
                                                                                                                                cGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTtta                    <  2:511476/150‑1 (MQ=255)
                                                                                                                                cGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTtta                    >  1:913750/1‑150 (MQ=255)
                                                                                                                                 taaaaTTGTCCGGCTCGTGCAAGACTTTCTCCAGGTTTGGTTTTAGTTCAAAGCGCCCGGGAATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTttat                   <  1:404811/149‑1 (MQ=255)
                                                                                                                                 taaaaTTGTCCGGCGCGTGCCAGACTTTCTCCATGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATTAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTttat                   <  1:410826/149‑1 (MQ=255)
                                                                                                                                 taaaaTTGGCGGGCGCGTGCCAGACTTTCTCGAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGATTAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTttat                   <  2:922752/149‑1 (MQ=255)
                                                                                                                                 gAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTttat                   <  2:110377/150‑1 (MQ=255)
                                                                                                                                  gagaaagTCTGGCACGCGCCGGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTttatt                  >  2:420669/7‑150 (MQ=255)
                                                                                                                                  aaaaTTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGc                                                      <  2:144032/114‑1 (MQ=255)
                                                                                                                                  aaaaTTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGc                                                      >  1:144032/1‑114 (MQ=255)
                                                                                                                                  aaaaTTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTtcatt                  <  2:675593/150‑1 (MQ=255)
                                                                                                                                   aaaTTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCATGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTttatta                 >  1:629754/1‑150 (MQ=255)
                                                                                                                                    gaaagTCTGGCACGCGCCGGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTtattat                >  2:907923/5‑150 (MQ=255)
                                                                                                                                     tttGTCCGGGGCGTGCCTGACTTTCTCCAGGGGTGGTTTTAGTGCATAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTGTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATtata               <  2:1661101/149‑1 (MQ=255)
                                                                                                                                       tGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATAtt             <  2:941703/150‑1 (MQ=255)
                                                                                                                                        gTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATAtta            >  2:815987/1‑149 (MQ=255)
                                                                                                                                         tCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCa           <  1:414558/150‑1 (MQ=255)
                                                                                                                                          ccGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCAt          >  1:1411668/1‑150 (MQ=255)
                                                                                                                                            ggCACGCGCCGGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATAc        <  1:907923/150‑1 (MQ=255)
                                                                                                                                             gcgcGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACg       >  2:390479/1‑150 (MQ=255)
                                                                                                                                               gcgTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACggg     <  1:1566732/150‑1 (MQ=255)
                                                                                                                                                 gTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGat   <  2:1514062/150‑1 (MQ=255)
                                                                                                                                                  tGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGata  >  1:486206/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
ACTGCGCTTCACGCGTCACCACCGGCGCGTTGGGCGTATCATCAGATGGAAGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATA  >  NC_000913/4428525‑4428820

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