Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 774,356 C→G 7.4% P16P (CCC→CCG)  ybgE → putative inner membrane protein in cydABX‑ybgE operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913774,3560CG 7.4% 601.2 / 31.5 175P16P (CCC→CCG) ybgEputative inner membrane protein in cydABX‑ybgE operon
Reads supporting (aligned to +/- strand):  ref base C (66/96);  new base G (7/6);  total (73/102)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

Only 100 of 185 total aligned reads displayed.

GCATGGATTCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGC  >  NC_000913/774208‑774503
                                                                                                                                                    |                                                                                                                                                    
gCATGGATTCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGcccc                                                                                                                                                      >  3:1087831/1‑149 (MQ=255)
       ttCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCggg                                                                                                                                               >  2:1006489/1‑149 (MQ=255)
       ttCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCggg                                                                                                                                               >  1:1326935/1‑149 (MQ=255)
          tGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGgcgc                                                                                                                                            <  4:955853/149‑1 (MQ=255)
           gggAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGTTTACCGCGACTTTGTATGCGGTAATGGACAAGGGCCCCCTGCGGGCGCt                                                                                                                                           >  1:465222/1‑149 (MQ=255)
            ggAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCtt                                                                                                                                          >  3:302806/1‑149 (MQ=255)
             gAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCGCAGGGGGCGCttg                                                                                                                                         <  3:1287392/149‑2 (MQ=255)
              aaCGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCGCAGGGGGCGCttgt                                                                                                                                        <  4:1427934/149‑3 (MQ=255)
              aaCGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑Tc                                                                                                                                       <  1:275166/149‑1 (MQ=255)
               aCGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑Tcc                                                                                                                                      >  4:1360764/1‑149 (MQ=255)
               aCGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAAAACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑Tcc                                                                                                                                      <  2:1117078/149‑1 (MQ=255)
                 gCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCtt                                                                                                                                    >  2:30488/1‑149 (MQ=255)
                  cttcttGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTc                                                                                                                                   >  4:224336/1‑149 (MQ=255)
                  cttcttGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTc                                                                                                                                   <  2:114932/149‑1 (MQ=255)
                  cttcttGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTc                                                                                                                                   <  1:791297/149‑1 (MQ=255)
                    tcttGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGt                                                                                                                                 <  2:1130234/149‑1 (MQ=255)
                         ccTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATgg                                                                                                                            <  2:620441/149‑1 (MQ=255)
                            gTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGgcgc                                                                                                                         >  2:251376/1‑149 (MQ=255)
                                gTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCt                                                                                                                     <  3:224336/149‑1 (MQ=255)
                                gTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCt                                                                                                                     <  2:300770/149‑1 (MQ=255)
                                 tttGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTg                                                                                                                    >  4:521977/1‑149 (MQ=255)
                                  ttGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGt                                                                                                                   <  2:56878/149‑1 (MQ=255)
                                       gTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCa                                                                                                              >  3:25384/1‑149 (MQ=255)
                                        tAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATCAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAg                                                                                                             >  1:330650/1‑149 (MQ=255)
                                           tCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGat                                                                                                          <  3:1355933/149‑1 (MQ=255)
                                           tCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATCAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGat                                                                                                          <  2:330650/149‑1 (MQ=255)
                                               cgcgCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGatgtat                                                                                                      <  4:874733/149‑1 (MQ=255)
                                                 cgcTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGatgtatgt                                                                                                    <  3:360229/149‑1 (MQ=255)
                                                    tGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGtttt                                                                                                 <  1:1238477/149‑1 (MQ=255)
                                                    tGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGtttt                                                                                                 <  4:302806/149‑1 (MQ=255)
                                                     ggCGCTTGAACACGTCGAATCAGACAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTg                                                                                                <  1:461103/149‑1 (MQ=255)
                                                       cgcTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTggg                                                                                              <  4:1285483/149‑1 (MQ=255)
                                                        gcTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAACAGGATGTATGTTTTGGGa                                                                                             <  4:1238052/149‑1 (MQ=255)
                                                        acTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGAAGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGa                                                                                             >  2:1255318/2‑149 (MQ=255)
                                                           tGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGAccc                                                                                          <  2:883293/149‑1 (MQ=255)
                                                              acacGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCAtc                                                                                       >  1:1287822/1‑149 (MQ=255)
                                                              acacGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCAtc                                                                                       >  3:1224474/1‑149 (MQ=255)
                                                               cacGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCAtct                                                                                      >  2:1214254/1‑149 (MQ=255)
                                                                acGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCAtctc                                                                                     >  2:1403838/1‑149 (MQ=255)
                                                                  gTCGAATCAGGGAAAGGCGGGCAAAAAGACAACTGATGGGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGt                                                                                   <  2:465222/149‑1 (MQ=255)
                                                                  gTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATgg                                                                                                                            >  3:1003590/1‑108 (MQ=255)
                                                                  gTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATgg                                                                                                                            <  4:1003590/108‑1 (MQ=255)
                                                                  gTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGt                                                                                   >  3:332067/1‑149 (MQ=255)
                                                                  gTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGt                                                                                   <  2:250699/149‑1 (MQ=255)
                                                                     gAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGtttt                                                                                <  4:1279094/149‑1 (MQ=255)
                                                                      aaTCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTg                                                                               >  2:533777/1‑149 (MQ=255)
                                                                        tCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTg                                                                                                >  3:758664/1‑130 (MQ=255)
                                                                        tCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTg                                                                                                <  4:758664/130‑1 (MQ=255)
                                                                           ggCGAAGCCGGTCAAGAAGACATCTGAGGCGTAAGATTATCGCGACATTGCACGCGGTCATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCTTGATGGCGCTTCTGTTAGCAGGATGTCTGTTTTGGGACCTATCTCGGTTTGCcggg                                                                          >  4:562186/1‑147 (MQ=255)
                                                                           ggCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCcgcg                                                                          >  1:688834/1‑149 (MQ=255)
                                                                            gCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGa                                                                         >  3:182582/1‑149 (MQ=255)
                                                                             cAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGaa                                                                        <  1:261820/149‑1 (MQ=255)
                                                                             cAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGaa                                                                        <  1:673415/149‑1 (MQ=255)
                                                                                   cGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAg                                                                  >  1:236102/1‑149 (MQ=255)
                                                                                    ggTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGt                                                                 <  4:556003/149‑1 (MQ=255)
                                                                                     gTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTg                                                                >  3:388723/1‑149 (MQ=255)
                                                                                      tCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGa                                                               <  2:108227/149‑1 (MQ=255)
                                                                                       cAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGaa                                                              <  2:434082/149‑1 (MQ=255)
                                                                                        aagaagACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAAc                                                             <  3:958104/149‑1 (MQ=255)
                                                                                         agaagaCATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACt                                                            >  1:536131/1‑149 (MQ=255)
                                                                                          gaagaCATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTg                                                           <  3:535498/149‑1 (MQ=255)
                                                                                          gaagaCATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTg                                                           >  4:199599/1‑149 (MQ=255)
                                                                                           aagaCATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTgg                                                          >  3:15076/1‑149 (MQ=255)
                                                                                             gaCATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGaa                                                        <  1:1006489/149‑1 (MQ=255)
                                                                                                    gatgaGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGc                                                 >  2:1299013/1‑149 (MQ=255)
                                                                                                       gaGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATg                                              <  1:367615/149‑1 (MQ=255)
                                                                                                         gTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGt                                            <  2:918637/149‑1 (MQ=255)
                                                                                                          tAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTg                                                                                                <  1:1405279/96‑1 (MQ=255)
                                                                                                          tAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGtt                                           >  1:447297/1‑149 (MQ=255)
                                                                                                             gATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTAt                                        >  3:1314532/1‑149 (MQ=255)
                                                                                                              aTTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTAtt                                       <  4:15076/149‑1 (MQ=255)
                                                                                                                   cgcgACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGa                                  <  1:30488/149‑1 (MQ=255)
                                                                                                                   cgcgACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGa                                  <  1:488180/149‑1 (MQ=255)
                                                                                                                    gcgACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTgg                                                          >  4:968134/1‑124 (MQ=255)
                                                                                                                    gcgACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGAt                                 <  4:25384/149‑1 (MQ=255)
                                                                                                                       aCTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGAtgtg                              >  4:349796/1‑149 (MQ=255)
                                                                                                                        cTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTgg                             <  3:490657/149‑1 (MQ=255)
                                                                                                                        cTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTgg                             <  1:800775/149‑1 (MQ=255)
                                                                                                                               tGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCt                      >  3:620341/1‑149 (MQ=255)
                                                                                                                                   gTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGGAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGc                  >  3:293560/1‑149 (MQ=255)
                                                                                                                                   gTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGc                  <  2:1242401/149‑1 (MQ=255)
                                                                                                                                       tGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGgt              <  1:1118895/149‑1 (MQ=255)
                                                                                                                                        ggACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGgtg             >  1:216712/1‑149 (MQ=255)
                                                                                                                                         gACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGgtgt            >  1:166302/1‑149 (MQ=255)
                                                                                                                                          gaaaGCGCCCGCAGGGGGCGCTTGTCCTTCGCGATGGCGCTTATGTTAGCAGGAGGTGTGGTTTGGGTGCCATCTCGTTTTGCCGCGAATACCAGTGAAATGGAAATCTGGCATGGTTTTTTGCTGATGTGGGCCGTCTGTTCTGgtgt            >  2:588697/3‑149 (MQ=255)
                                                                                                                                          gaaaGCGCCCGCAGGGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGt                                            >  4:1340901/3‑116 (MQ=255)
                                                                                                                                           cAAGCGCCCCCTGCGGGCGCTC‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGAGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGGGCTGGTGTGa          <  3:898104/149‑1 (MQ=255)
                                                                                                                                           aaaGCGCCCGCAGGGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGa          >  4:687100/2‑149 (MQ=255)
                                                                                                                                            aaGCGCCCGCAGGGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGAtgt                               <  2:1037321/127‑1 (MQ=255)
                                                                                                                                            aaGCGCCCGCAGGGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGAtgt                               >  1:1037321/1‑127 (MQ=255)
                                                                                                                                            aaGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCtgt                    >  2:385587/1‑138 (MQ=255)
                                                                                                                                            aaGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCtgt                    <  1:385587/138‑1 (MQ=255)
                                                                                                                                              gcgcCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTc       <  3:1358832/149‑1 (MQ=255)
                                                                                                                                              gcgcCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTc       >  4:477272/1‑149 (MQ=255)
                                                                                                                                               cgcCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCa      <  1:1374394/149‑1 (MQ=255)
                                                                                                                                                 cccccTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACg    >  4:1152966/1‑149 (MQ=255)
                                                                                                                                                  ccccTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTg                                                                                                >  2:1179617/1‑56 (MQ=255)
                                                                                                                                                  ccccTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTg                                                                                                <  1:1179617/56‑1 (MQ=255)
                                                                                                                                                  ccccTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGAt         >  1:1041722/1‑143 (MQ=255)
                                                                                                                                                   cccTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGc  >  2:791748/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCATGGATTCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGC  >  NC_000913/774208‑774503

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