Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,823,839 T→G 7.1% 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 7.1% 210.6 / 5.2 79intergenic (‑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 (34/39);  new base G (4/2);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

ATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCAAACC  >  NC_000913/2823690‑2823985
                                                                                                                                                     |                                                                                                                                                  
ttgGAGACTTTACCAAATTGTTTCTCAATCTGGCGCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTGAATATGTTCTGTTGAAGCAATTATACTGTATGCt                                                                                                                                                    <  1:765554/149‑1 (MQ=255)
     gCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATAc                                                                                                                                               <  1:364127/150‑1 (MQ=255)
      ccTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCATTTATACTGTATGCTCATACa                                                                                                                                              >  1:178060/1‑150 (MQ=255)
            ccAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCa                                                                                                                                        <  2:981224/150‑1 (MQ=255)
             cAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGAGCATACAGTAtaat                                                                                                                                       <  2:320997/150‑4 (MQ=255)
             cAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGATTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGAAATTATAATGTATGCTCATACAGTATCaa                                                                                                                                       >  2:490758/1‑150 (MQ=255)
               aaTTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGAGCATACAGTAtaattg                                                                                                                                     <  1:1120645/150‑6 (MQ=255)
                  tGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGtt                                                                                                                                  <  2:705600/150‑1 (MQ=255)
                   gTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGttt                                                                                                                                 <  2:982991/150‑1 (MQ=255)
                       ctcAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTa                                                                                                                             <  1:250859/150‑1 (MQ=255)
                       ctcAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATAATGTATGCTCATACAGTCTCAAGTGTTTTGTa                                                                                                                             >  2:636523/1‑150 (MQ=255)
                            tCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACAGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAAt                                                                                                                        <  1:161736/150‑1 (MQ=255)
                             cTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAAtt                                                                                                                       >  1:1045244/1‑150 (MQ=255)
                             cTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAAtt                                                                                                                       <  1:1219874/150‑1 (MQ=255)
                             cTGGCCCAGTGCAGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAAtt                                                                                                                       <  1:107891/150‑1 (MQ=255)
                                 cccAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAAttgttg                                                                                                                   >  2:737094/1‑150 (MQ=255)
                                 cccAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAAttgttg                                                                                                                   >  2:860431/1‑150 (MQ=255)
                                 cccAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTGGAAGCAATTATACTGTTTGCTCATACAGTATCAAGTGTTTTGTAGAAAttgttg                                                                                                                   <  1:314758/150‑1 (MQ=255)
                                  ccAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGc                                                                                                                  <  1:436365/150‑1 (MQ=255)
                                   cAGTGCTGCCGCCAATGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGcc                                                                                                                 >  2:749821/1‑150 (MQ=255)
                                           ccgccAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGGGTTTTGTAGAAATTGTTGCCACAAGGTc                                                                                                         <  1:1158916/150‑1 (MQ=255)
                                               cAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCTAGATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAATTTGTTGCCACAAGGTCTGCa                                                                                                     >  2:850814/1‑150 (MQ=255)
                                                 aggCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATTGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGTTCTGCAAt                                                                                                   <  2:1286235/148‑1 (MQ=255)
                                                       tCTGTTTGTTTTCGTCGATAGCCATTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACg                                                                                            <  2:353886/150‑1 (MQ=255)
                                                        cTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGTTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAATGTCTGCAATGCATACg                                                                                            <  2:1072931/150‑1 (MQ=255)
                                                        cTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACg                                                                                            >  2:570995/1‑150 (MQ=255)
                                                          gtttgtttTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCTATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCa                                                                                          >  2:967837/1‑150 (MQ=255)
                                                             tgtttTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTa                                                                                       >  1:287348/1‑150 (MQ=255)
                                                                tttCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGcc                                                                                    <  1:615752/150‑1 (MQ=255)
                                                                 ttCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCt                                                                                   >  2:225880/1‑150 (MQ=255)
                                                                 ttCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGATCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCt                                                                                   <  1:812460/150‑1 (MQ=255)
                                                                    gtcgATAGCCATTTTTCCTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTgac                                                                                <  2:838528/150‑1 (MQ=255)
                                                                     tcgATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTgacg                                                                               <  1:334210/150‑1 (MQ=255)
                                                                       gATAGCCATTTTTACTCCTGTCAGGACGGGTAATCCAGGCTAGTCAATATGGTCTGATGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACACGGTCTGCACTGCATACGCAGTAGCCTgacgac                                                                             <  2:499329/150‑1 (MQ=255)
                                                                        aTAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTGCTGTTGAAGCAATTATACTGTATGCTCATCCAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGacgacg                                                                            <  1:935948/150‑1 (MQ=255)
                                                                           gCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCAc                                                                         >  2:821730/1‑150 (MQ=255)
                                                                           gCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCAc                                                                         >  2:291896/1‑150 (MQ=255)
                                                                            ccATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCAcc                                                                        <  2:360916/150‑1 (MQ=255)
                                                                             cATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGAGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTTACGACGCACCg                                                                       >  2:567501/1‑150 (MQ=255)
                                                                                  ttACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCa                                                                  <  1:1111886/150‑1 (MQ=255)
                                                                                   tACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCAc                                                                 >  2:1312033/1‑150 (MQ=255)
                                                                                     cTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACgg                                                               <  2:368984/150‑1 (MQ=255)
                                                                                     cTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACgg                                                               >  2:757356/1‑150 (MQ=255)
                                                                                     cTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACgg                                                               <  2:460074/150‑1 (MQ=255)
                                                                                       ccTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTc                                                             >  1:1240108/1‑150 (MQ=255)
                                                                                       ccTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTc                                                             >  2:287288/1‑150 (MQ=255)
                                                                                                ccGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGaa                                                    >  1:463630/1‑150 (MQ=255)
                                                                                                 cGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAg                                                   >  2:364380/1‑150 (MQ=255)
                                                                                                  gggTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGc                                                  >  2:282570/1‑150 (MQ=255)
                                                                                                   ggTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCa                                                 >  1:31164/1‑150 (MQ=255)
                                                                                                      aaTACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTc                                              >  2:1007230/1‑150 (MQ=255)
                                                                                                       aTACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTcc                                             >  1:407014/1‑150 (MQ=255)
                                                                                                            ggATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCccgccg                                        <  2:430295/150‑1 (MQ=255)
                                                                                                            ggATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCccgccg                                        <  1:41786/150‑1 (MQ=255)
                                                                                                             gATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCgg                                       >  2:1312955/1‑150 (MQ=255)
                                                                                                              aTAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCggg                                      <  2:473499/150‑1 (MQ=255)
                                                                                                              aTAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCggg                                      <  2:1015403/150‑1 (MQ=255)
                                                                                                              aTAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCggg                                      <  2:1293975/150‑1 (MQ=255)
                                                                                                              aTAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGGCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCggg                                      <  2:296346/150‑1 (MQ=255)
                                                                                                               tAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGt                                     <  1:859684/150‑1 (MQ=255)
                                                                                                                aGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTa                                    >  1:134191/1‑150 (MQ=255)
                                                                                                                    aaTATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGc                                <  1:182053/150‑1 (MQ=255)
                                                                                                                      tatGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACAGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGccg                              >  1:95066/1‑149 (MQ=255)
                                                                                                                            cTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACcg                        >  2:1314127/1‑150 (MQ=255)
                                                                                                                            atGTTGAAGCAATTATACTGAATGCTCCTACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGGCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACcg                        <  1:907932/149‑1 (MQ=255)
                                                                                                                             tGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACAGCATCACGGTCGCCGCTGAAGCATTCCCGACGGGTAATGCCTTCACcgc                       >  1:191508/1‑150 (MQ=255)
                                                                                                                                   acacttgataCTGTATGAGCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAg                 >  2:618517/8‑150 (MQ=255)
                                                                                                                                   aGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAg                 >  1:39872/1‑150 (MQ=255)
                                                                                                                                      aaTTATCCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACCAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGGCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTgg              <  1:863906/150‑1 (MQ=255)
                                                                                                                                        tgataCTGTATGAGCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCa            >  1:919476/3‑150 (MQ=255)
                                                                                                                                        tgataCTGTATGAGCATACAGTATCAAGTGTTTTGTAGAAAATGTTGCCACAAGGTCTGCAATGCATACTAAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGTGCACTGGCa            >  1:185050/3‑150 (MQ=255)
                                                                                                                                         tataCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCaa           >  2:781665/1‑150 (MQ=255)
                                                                                                                                         tataCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAAAGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCCCCGCGGGCAGTGGCaa           >  1:236286/1‑150 (MQ=255)
                                                                                                                                          ataCTGTATGAGCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTCGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCaaa          >  2:1012535/1‑150 (MQ=255)
                                                                                                                                            aCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAg        <  1:425663/150‑1 (MQ=255)
                                                                                                                                            aCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAg        <  1:225880/150‑1 (MQ=255)
                                                                                                                                                  tGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCAAAcc  >  1:140443/1‑150 (MQ=255)
                                                                                                                                                  tGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCAAAcc  <  2:1270365/150‑1 (MQ=255)
                                                                                                                                                  tGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCAAAcc  <  1:7381/150‑1 (MQ=255)
                                                                                                                                                  tGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCAAAcc  <  2:866333/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCAAACC  >  NC_000913/2823690‑2823985

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