Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,285,321 T→A 7.3% intergenic (‑108/+92) yjcH ← / ← acs DUF485 family inner membrane protein/acetyl‑CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,285,3210TA 7.3% 378.3 / 16.0 109intergenic (‑108/+92)yjcH/acsDUF485 family inner membrane protein/acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base T (51/50);  new base A (5/3);  total (56/53)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

Only 100 of 109 total aligned reads displayed.

GCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCC  >  NC_000913/4285174‑4285466
                                                                                                                                                   |                                                                                                                                                 
gcgcATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCAGAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTGATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTa                                                                                                                                                  <  1:762747/149‑1 (MQ=255)
 cgcATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTAc                                                                                                                                                 <  2:256484/149‑1 (MQ=255)
   cATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATGCGGCGTAAAcgc                                                                                                                                               >  1:316436/1‑149 (MQ=255)
    aTTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                                                              >  1:864615/1‑149 (MQ=255)
        tcttctATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATGCGGCGTAAAcgccttat                                                                                                                                          <  2:316436/149‑5 (MQ=255)
          ttctATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATGCGGCGTAAAcgccttat                                                                                                                                          >  2:67341/1‑143 (MQ=255)
          ttctATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATGCGGCGTAAAcgccttat                                                                                                                                          <  1:67341/147‑5 (MQ=255)
          ttctATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATcc                                                                                                                                        >  4:866334/1‑149 (MQ=255)
          ttctATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATcc                                                                                                                                        >  2:602176/1‑149 (MQ=255)
              aTCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                    <  4:433811/149‑2 (MQ=255)
              aTCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGca                                                                                                                                    >  4:247249/1‑149 (MQ=255)
                  gCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGcaatca                                                                                                                                >  4:828732/1‑149 (MQ=255)
                    tGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGcaatcaat                                                                                                                              >  3:300899/1‑149 (MQ=255)
                      aTAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGc                                                                                                                            >  4:255893/1‑149 (MQ=255)
                       tAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCg                                                                                                                                       <  1:590399/137‑1 (MQ=255)
                          aTAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGa                                                                                                                        <  3:8509/149‑1 (MQ=255)
                                  aTTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCt                                                                                                                          >  4:393773/1‑139 (MQ=255)
                                  aTTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCt                                                                                                                          <  3:807374/139‑1 (MQ=255)
                                  aTTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCt                                                                                                                          >  4:807374/1‑139 (MQ=255)
                                  aTTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCt                                                                                                                          <  3:393773/139‑1 (MQ=255)
                                         acaGATTCTCCAGAGGTTATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgcg                                                                                                         >  1:403421/1‑149 (MQ=255)
                                         acaGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgcg                                                                                                         >  4:338710/1‑149 (MQ=255)
                                         acaGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgcg                                                                                                         >  4:791190/1‑149 (MQ=255)
                                         acaGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgcg                                                                                                         >  1:27642/1‑149 (MQ=255)
                                         acaGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgcg                                                                                                         >  4:768696/1‑149 (MQ=255)
                                           aGATTTTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTc                                                                                                       <  3:729543/149‑1 (MQ=255)
                                           aGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTc                                                                                                       <  1:174024/149‑1 (MQ=255)
                                           aGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTc                                                                                                       >  2:316101/1‑149 (MQ=255)
                                            gATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCt                                                                                                      >  4:144405/1‑149 (MQ=255)
                                            gATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCt                                                                                                      >  4:20832/1‑149 (MQ=255)
                                             aTTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGc                                                                                                                >  1:22407/1‑138 (MQ=255)
                                              ttCTCCAGAGGGAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTa                                                                                                    >  2:24902/1‑149 (MQ=255)
                                               tctcCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTAt                                                                                                   <  4:711361/149‑1 (MQ=255)
                                               tctcCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTAt                                                                                                   <  3:530239/149‑1 (MQ=255)
                                               tctcCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTAt                                                                                                   <  1:289063/149‑1 (MQ=255)
                                                 tcCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTGAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCa                                                                                                 >  2:392085/1‑149 (MQ=255)
                                                  ccAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAg                                                                                                >  4:269492/1‑149 (MQ=255)
                                                    agagGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGc                                                                                              <  2:728767/149‑1 (MQ=255)
                                                     gagGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGcc                                                                                             <  4:396891/149‑1 (MQ=255)
                                                              tagcgcttGTTTTAATTCCCGCTCCCTAATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAAcct                                                                                    <  4:798558/143‑2 (MQ=255)
                                                               aGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGt                                                                                   <  3:526392/149‑1 (MQ=255)
                                                               aGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGt                                                                                   <  4:306859/149‑1 (MQ=255)
                                                                 ggATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTa                                                                                 >  1:230888/1‑149 (MQ=255)
                                                                  gATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTAc                                                                                <  2:729286/149‑1 (MQ=255)
                                                                  gATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTAc                                                                                <  2:864615/149‑1 (MQ=255)
                                                                   aTTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTAcc                                                                               >  3:703088/1‑149 (MQ=255)
                                                                   aTTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTAcc                                                                               >  4:547940/1‑149 (MQ=255)
                                                                    ttGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGc                                                                                              >  3:202603/1‑133 (MQ=255)
                                                                    ttGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGc                                                                                              <  4:202603/133‑1 (MQ=255)
                                                                     tGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGa                                                                             <  3:13583/149‑1 (MQ=255)
                                                                      gTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGAc                                                                            >  3:583028/1‑149 (MQ=255)
                                                                        tttAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTc                                                                          >  2:876241/1‑149 (MQ=255)
                                                                         ttAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCg                                                                         >  2:88869/1‑149 (MQ=255)
                                                                           aaTTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCa                                                                       >  2:30656/1‑149 (MQ=255)
                                                                            aTTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCAt                                                                      <  3:247249/149‑1 (MQ=255)
                                                                             ttCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGACTTATCAGGCCTACAAACCGTTACCGACTCGCATc                                                                     >  2:417396/1‑149 (MQ=255)
                                                                               cccGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCgg                                                                   <  4:49246/149‑1 (MQ=255)
                                                                               cccGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCgg                                                                   <  3:496436/149‑1 (MQ=255)
                                                                               cccGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCgg                                                                   <  3:866334/149‑1 (MQ=255)
                                                                                   cTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCaa                                                               <  3:255893/149‑1 (MQ=255)
                                                                                    tCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAAt                                                              >  3:423581/1‑149 (MQ=255)
                                                                                        ttATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATtgtg                                                          >  4:868966/1‑149 (MQ=255)
                                                                                        ttATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATtgtg                                                          <  3:868966/149‑1 (MQ=255)
                                                                                          aTGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTggg                                                        >  4:742686/1‑149 (MQ=255)
                                                                                          aTGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTggg                                                        >  2:820768/1‑149 (MQ=255)
                                                                                              gagaAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTAc                                                    >  1:839165/1‑149 (MQ=255)
                                                                                               agaAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACg                                                   >  1:656946/1‑149 (MQ=255)
                                                                                               agaAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACg                                                   >  3:871230/1‑149 (MQ=255)
                                                                                                 aaGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTgg                                                         >  3:855410/1‑141 (MQ=255)
                                                                                                 aaGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTgg                                                         <  4:855410/141‑1 (MQ=255)
                                                                                                  aGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATg                                                <  3:541838/149‑1 (MQ=255)
                                                                                                    gTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGc                                              <  3:742686/149‑1 (MQ=255)
                                                                                                    gTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGc                                              <  1:820768/149‑1 (MQ=255)
                                                                                                            tCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGAt                                      <  1:316101/149‑1 (MQ=255)
                                                                                                                 tatcccTTGCCGAATGTAAGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTCCCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCt                                 >  4:71353/4‑149 (MQ=255)
                                                                                                                  aaCCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTg                                >  1:143440/1‑149 (MQ=255)
                                                                                                                   aCCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGc                               <  4:583028/149‑1 (MQ=255)
                                                                                                                     ccTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCtt                             >  1:750029/1‑149 (MQ=255)
                                                                                                                     ccTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCtt                             >  2:815180/1‑149 (MQ=255)
                                                                                                                        tGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTtctc                          <  3:71353/149‑1 (MQ=255)
                                                                                                                        tGCCGAATGTAGGCCGGATAAGCCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTtctc                          <  4:64603/149‑1 (MQ=255)
                                                                                                                            gAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTAc                                                                                >  2:544466/1‑91 (MQ=255)
                                                                                                                            gAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTAc                                                                                <  1:544466/91‑1 (MQ=255)
                                                                                                                            gAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCa                      >  1:280737/1‑149 (MQ=255)
                                                                                                                             aaTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCaa                     <  4:871230/149‑1 (MQ=255)
                                                                                                                             aaTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCaa                     <  3:835220/149‑1 (MQ=255)
                                                                                                                             aaTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCaa                     <  3:20832/149‑1 (MQ=255)
                                                                                                                             aaTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCaa                     >  2:327532/1‑149 (MQ=255)
                                                                                                                             aaTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCaa                     <  1:75184/149‑1 (MQ=255)
                                                                                                                             aaTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCaa                     <  1:392085/149‑1 (MQ=255)
                                                                                                                              ttgattgCCGGATGCGGCGTAAACGCCTTATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTATCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAag                    >  2:3015/7‑149 (MQ=37)
                                                                                                                                  ttgCCGGATGCGGCGTAAACGCCTTATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAagcagc                >  4:236733/3‑149 (MQ=37)
                                                                                                                                  aGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAagcagc                <  3:828732/149‑1 (MQ=255)
                                                                                                                                      cGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTtct            <  1:876241/149‑1 (MQ=255)
                                                                                                                                       ggATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTtctc           <  2:750029/149‑1 (MQ=255)
                                                                                                                                        gATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCg          >  2:173915/1‑149 (MQ=255)
                                                                                                                                         aTAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGa         <  4:806406/149‑1 (MQ=255)
                                                                                                                                               cGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGc   <  3:768696/149‑1 (MQ=255)
                                                                                                                                               cGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGc   <  3:338710/149‑1 (MQ=255)
                                                                                                                                                gTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGcc  <  2:360165/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCC  >  NC_000913/4285174‑4285466

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