Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,823,839 T→G 6.7% intergenic (‑70/+10) recA ← / ← pncC DNA recombination and repair protein; ssDNA‑dependent ATPase; synaptase; ssDNA and dsDNA binding protein; ATP‑dependent homologous DNA strand exchanger; recombinase A; LexA autocleavage cofactor/nicotinamide‑nucleotide amidohydrolase; NMN amidohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,823,8390TG 6.7% 363.7 / 17.0 104intergenic (‑70/+10)recA/pncCDNA recombination and repair protein; ssDNA‑dependent ATPase; synaptase; ssDNA and dsDNA binding protein; ATP‑dependent homologous DNA strand exchanger; recombinase A; LexA autocleavage cofactor/nicotinamide‑nucleotide amidohydrolase; NMN amidohydrolase
Reads supporting (aligned to +/- strand):  ref base T (49/48);  new base G (3/4);  total (52/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

Only 100 of 105 total aligned reads displayed.

GGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCAAACCAG  >  NC_000913/2823692‑2823987
                                                                                                                                                   |                                                                                                                                                     
ggAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTc                                                                                                                                                      >  2:646548/1‑149 (MQ=255)
     cTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACa                                                                                                                                                 >  2:535661/1‑149 (MQ=255)
      tttACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAg                                                                                                                                                >  1:653401/1‑149 (MQ=255)
        tACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGAGCATACAGTa                                                                                                                                              <  4:332498/149‑1 (MQ=255)
          ccAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATc                                                                                                                                            <  3:254838/149‑1 (MQ=255)
           cAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCa                                                                                                                                           >  3:125960/1‑149 (MQ=255)
           cAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCa                                                                                                                                           >  2:127437/1‑149 (MQ=255)
           cAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCa                                                                                                                                           >  3:611805/1‑149 (MQ=255)
            aaaTTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCaa                                                                                                                                          >  4:68799/1‑149 (MQ=255)
            aaaTTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCaa                                                                                                                                          <  4:655701/149‑1 (MQ=255)
              aTTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAgt                                                                                                                                        >  4:364019/1‑149 (MQ=255)
               ttGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGAGCATACAGTAtaattgc                                                                                                                                       <  3:383132/149‑7 (MQ=255)
                      tcAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTa                                                                                                                                <  2:591675/149‑1 (MQ=255)
                           cTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAAt                                                                                                                           <  4:537664/149‑1 (MQ=255)
                              gCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAAttgt                                                                                                                        <  2:292075/149‑1 (MQ=255)
                              gCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAAttgt                                                                                                                        <  2:237862/149‑1 (MQ=255)
                               cccAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAAttgtt                                                                                                                       <  1:127437/149‑1 (MQ=255)
                                ccAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAAttgttg                                                                                                                      >  3:411205/1‑149 (MQ=255)
                                   gTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCaa                                                                                                                                          >  2:635023/1‑126 (MQ=255)
                                   gTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCaa                                                                                                                                          <  1:635023/126‑1 (MQ=255)
                                        gccgccAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAgg                                                                                                              <  3:68799/149‑1 (MQ=255)
                                             cAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGc                                                                                                         >  1:183172/1‑149 (MQ=255)
                                             cAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGc                                                                                                         <  1:181727/149‑1 (MQ=255)
                                             cAACGCTTTCTGTTTGTTTTCGTCAATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGc                                                                                                         >  3:570001/1‑149 (MQ=255)
                                              aaCGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCa                                                                                                        >  2:517052/1‑149 (MQ=255)
                                                   tttCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACa                                                                                                                                                 >  2:134731/1‑103 (MQ=255)
                                                   tttCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACa                                                                                                                                                 <  1:134731/103‑1 (MQ=255)
                                                   tttCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCa                                                                                                   <  1:535661/149‑1 (MQ=255)
                                                   tttCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCa                                                                                                   >  4:646948/1‑149 (MQ=255)
                                                     tCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATa                                                                                                 <  1:310045/149‑1 (MQ=255)
                                                         tttgtttTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGttt                                                                                                                                    <  1:47296/110‑1 (MQ=255)
                                                         tttgtttTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCa                                                                                             >  2:67805/1‑149 (MQ=255)
                                                             ttttCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTa                                                                                          <  1:578443/148‑1 (MQ=255)
                                                             ttttCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTa                                                                                          >  2:578443/1‑148 (MQ=255)
                                                                tcgtcgATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCt                                                                                      >  3:280887/1‑149 (MQ=255)
                                                                tcgtcgATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCt                                                                                      <  4:333510/149‑1 (MQ=255)
                                                                tcgtcgATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCt                                                                                      <  4:611805/149‑1 (MQ=255)
                                                                tcgtcgATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCt                                                                                      >  4:9000/1‑149 (MQ=255)
                                                                     gATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTa                                                                                                                                >  3:94209/1‑102 (MQ=255)
                                                                     gATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATAc                                                                                                <  4:243791/134‑1 (MQ=255)
                                                                     gATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATAc                                                                                                >  3:243791/1‑134 (MQ=255)
                                                                        aGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGc                                                                              >  1:374965/1‑149 (MQ=255)
                                                                         gCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCa                                                                             <  1:418904/149‑1 (MQ=255)
                                                                          ccATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATa                                                                                                 <  1:669317/128‑1 (MQ=255)
                                                                          ccATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATa                                                                                                 >  2:669317/1‑128 (MQ=255)
                                                                          ccATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCAc                                                                            >  2:478088/1‑149 (MQ=255)
                                                                           cATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCAcc                                                                           >  4:388634/1‑149 (MQ=255)
                                                                            aTTTTTACTCCTGTCATGCCGGGTAATACCGGATAGACAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAACTGTTGCCACAAGGTCTGCAATGCATACGCAGCAGCCTGACGACGCACCg                                                                          >  4:583612/1‑149 (MQ=255)
                                                                                  aCTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCAc                                                                    >  1:32880/1‑149 (MQ=255)
                                                                                   cTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACg                                                                   >  2:257994/1‑149 (MQ=255)
                                                                                    tCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACgg                                                                  >  3:64883/1‑149 (MQ=255)
                                                                                      cTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTc                                                                >  4:300277/1‑149 (MQ=255)
                                                                                         tCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTcgcc                                                             <  3:388634/149‑1 (MQ=255)
                                                                                         tCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTAAACTGTATGCTCAGACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTcgcc                                                             >  4:658097/1‑149 (MQ=255)
                                                                                               cGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTcgc                                                              <  2:392240/142‑1 (MQ=255)
                                                                                               cGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTcgc                                                              >  1:392240/1‑142 (MQ=255)
                                                                                               cGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGaa                                                       <  1:67805/149‑1 (MQ=255)
                                                                                                     aTACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAATGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTc                                                <  1:473940/149‑1 (MQ=255)
                                                                                                      tACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTcc                                               >  3:246907/1‑149 (MQ=255)
                                                                                                            aTAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCgg                                         <  3:658097/149‑1 (MQ=255)
                                                                                                              aGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTATTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGt                                       >  4:48145/1‑149 (MQ=255)
                                                                                                                tCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTaa                                     >  4:691719/1‑149 (MQ=255)
                                                                                                                 cAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTACCCGGCgggaat                                     <  1:478088/149‑4 (MQ=255)
                                                                                                                 cAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAAt                                    <  3:85166/149‑1 (MQ=255)
                                                                                                                   atatGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGc                                  <  4:481736/149‑1 (MQ=255)
                                                                                                                   atatGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGc                                  >  1:476741/1‑149 (MQ=255)
                                                                                                                    tatGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGcc                                 <  2:573592/149‑1 (MQ=255)
                                                                                                                        ttCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCTCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCa                             >  1:375588/1‑149 (MQ=255)
                                                                                                                        ttCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCa                             >  1:76495/1‑149 (MQ=255)
                                                                                                                            gTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACg                                                                   >  1:394272/1‑108 (MQ=255)
                                                                                                                            gTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACg                                                                   <  2:394272/108‑1 (MQ=255)
                                                                                                                             ttGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACcgcg                        >  1:17903/1‑149 (MQ=255)
                                                                                                                                aaGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGc                     >  1:430161/1‑149 (MQ=255)
                                                                                                                                   cAATTATACTTTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGt                  <  4:514518/149‑1 (MQ=255)
                                                                                                                                   cAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGt                  <  4:64883/149‑1 (MQ=255)
                                                                                                                                    aaTTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCTCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTg                 <  2:375588/149‑1 (MQ=255)
                                                                                                                                    aaTTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTg                 <  2:76495/149‑1 (MQ=255)
                                                                                                                                    aaTTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTg                 <  1:646548/149‑1 (MQ=255)
                                                                                                                                    aaTTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTg                 <  1:517052/149‑1 (MQ=255)
                                                                                                                                      ttATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCggg                      <  2:671881/142‑1 (MQ=255)
                                                                                                                                      ttATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCggg                      >  1:671881/1‑142 (MQ=255)
                                                                                                                                      ttATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGc               >  3:71901/1‑149 (MQ=255)
                                                                                                                                        ataCTGTATGAGCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCaa             >  3:565963/1‑149 (MQ=255)
                                                                                                                                         taCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGGCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCaaa            >  4:232249/1‑149 (MQ=255)
                                                                                                                                         taCTGTATGAGCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGaa                                                       <  3:307455/107‑1 (MQ=255)
                                                                                                                                         taCTGTATGAGCATACAGTATCAAGTGTTCTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCaaa            >  1:543300/1‑149 (MQ=255)
                                                                                                                                          aCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCaaaa           <  4:480890/149‑1 (MQ=255)
                                                                                                                                          aCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCaaaa           <  1:257994/149‑1 (MQ=255)
                                                                                                                                          aCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCaaaa           <  4:246907/149‑1 (MQ=255)
                                                                                                                                          aCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCaaaa           <  4:125960/149‑1 (MQ=255)
                                                                                                                                           cTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGaa                                                       <  3:415263/105‑1 (MQ=255)
                                                                                                                                           cTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGaa                                                       >  4:415263/1‑105 (MQ=255)
                                                                                                                                           cTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAg          >  1:137416/1‑149 (MQ=255)
                                                                                                                                              tATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCaa       >  3:273198/1‑149 (MQ=255)
                                                                                                                                              tATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCaa       <  4:280887/149‑1 (MQ=255)
                                                                                                                                              tATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCaa       <  3:9000/149‑1 (MQ=255)
                                                                                                                                               aTGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCaaa      <  2:653401/149‑1 (MQ=255)
                                                                                                                                                tGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCAAAc     >  4:517623/1‑149 (MQ=255)
                                                                                                                                                  cTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCAAACCa   <  2:183172/149‑1 (MQ=255)
                                                                                                                                                   tCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCAAACCAg  >  1:430449/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCAAACCAG  >  NC_000913/2823692‑2823987

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