Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,285,320 T→A 7.3% intergenic (‑107/+93) 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,3200TA 7.3% 375.8 / 15.8 109intergenic (‑107/+93)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 = 9.49e-01

Only 100 of 110 total aligned reads displayed.

ATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCC  >  NC_000913/4285172‑4285466
                                                                                                                                                    |                                                                                                                                                  
aTGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGtt                                                                                                                                                    >  4:526392/1‑149 (MQ=255)
  gcgcATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCAGAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTGATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTa                                                                                                                                                  <  1:762747/149‑1 (MQ=255)
   cgcATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTAc                                                                                                                                                 <  2:256484/149‑1 (MQ=255)
     cATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATGCGGCGTAAAcgc                                                                                                                                               >  1:316436/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)
                         tAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCg                                                                                                                                       >  2:590399/1‑137 (MQ=255)
                            aTAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGa                                                                                                                        <  3:8509/149‑1 (MQ=255)
                                    aTTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCt                                                                                                                          <  3:807374/139‑1 (MQ=255)
                                    aTTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCt                                                                                                                          <  3:393773/139‑1 (MQ=255)
                                    aTTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCt                                                                                                                          >  4:807374/1‑139 (MQ=255)
                                           acaGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgcg                                                                                                         >  4:338710/1‑149 (MQ=255)
                                           acaGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgcg                                                                                                         >  4:768696/1‑149 (MQ=255)
                                           acaGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgcg                                                                                                         >  4:791190/1‑149 (MQ=255)
                                           acaGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgcg                                                                                                         >  1:27642/1‑149 (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                                                                                                                <  2:22407/138‑1 (MQ=255)
                                               aTTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGc                                                                                                                >  1:22407/1‑138 (MQ=255)
                                                ttCTCCAGAGGGAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTa                                                                                                    >  2:24902/1‑149 (MQ=255)
                                                 tctcCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTAt                                                                                                   <  1:289063/149‑1 (MQ=255)
                                                 tctcCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTAt                                                                                                   <  4:711361/149‑1 (MQ=255)
                                                 tctcCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTAt                                                                                                   <  3:530239/149‑1 (MQ=255)
                                                 tctcCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTAt                                                                                                   <  4:300899/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)
                                                                tAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCg                                                                                    >  4:734646/1‑149 (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:864615/149‑1 (MQ=255)
                                                                    gATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTAc                                                                                <  2:729286/149‑1 (MQ=255)
                                                                     aTTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTAcc                                                                               >  4:547940/1‑149 (MQ=255)
                                                                     aTTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTAcc                                                                               >  3:703088/1‑149 (MQ=255)
                                                                      ttGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGc                                                                                              >  3:202603/1‑133 (MQ=255)
                                                                      ttGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGc                                                                                              <  4:202603/133‑1 (MQ=255)
                                                                      ttGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCg                                                                              >  3:852998/1‑149 (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)
                                                                                                gagaAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTAc                                                    >  1:839165/1‑149 (MQ=255)
                                                                                                 agaAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACg                                                   >  1:656946/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                                              <  1:820768/149‑1 (MQ=255)
                                                                                                      gTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGc                                              <  3:742686/149‑1 (MQ=255)
                                                                                                              tCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGAt                                      <  1:316101/149‑1 (MQ=255)
                                                                                                                  gTAACCCTTGCTGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGcc                                  >  3:541238/1‑149 (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                                                                                <  1:544466/91‑1 (MQ=255)
                                                                                                                              gAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTAc                                                                                >  2:544466/1‑91 (MQ=255)
                                                                                                                              gAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCa                      >  1:280737/1‑149 (MQ=255)
                                                                                                                               aaTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCaa                     <  3:835220/149‑1 (MQ=255)
                                                                                                                               aaTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCaa                     >  2:327532/1‑149 (MQ=255)
                                                                                                                               aaTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCaa                     <  1:75184/149‑1 (MQ=255)
                                                                                                                               aaTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCaa                     <  3:20832/149‑1 (MQ=255)
                                                                                                                               aaTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCaa                     <  1:392085/149‑1 (MQ=255)
                                                                                                                               aaTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCaa                     <  4:871230/149‑1 (MQ=255)
                                                                                                                                ttgattgCCGGATGCGGCGTAAACGCCTTATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTATCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAag                    >  2:3015/7‑149 (MQ=37)
                                                                                                                                   tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAagcag                 <  3:144405/149‑1 (MQ=255)
                                                                                                                                    ttgCCGGATGCGGCGTAAACGCCTTATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAagcagc                >  1:544849/3‑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)
                                                                                                                                                    |                                                                                                                                                  
ATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCC  >  NC_000913/4285172‑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: