Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,823,839 T→G 15.4% 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,8390TG15.4% 236.2 / 38.9 97intergenic (‑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 (40/42);  new base G (9/6);  total (49/48)
Fisher's exact test for biased strand distribution p-value = 5.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

CTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGC  >  NC_000913/2823742‑2823937
                                                                                                 |                                                                                                  
cTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCa                                                                                                  >  1:1045993/1‑100 (MQ=255)
cTTTCTGTTTGGTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATAAGTTCTGTTGAGGCAATTATACTGTATGCTCa                                                                                                  >  2:512456/1‑100 (MQ=255)
  ttCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATTCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATa                                                                                                <  1:442028/100‑1 (MQ=255)
  ttCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTAATa                                                                                                >  1:1747705/1‑100 (MQ=255)
   tCTGTTTGTTTTCGTCGATAGCCATTTTTACGCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATAc                                                                                               <  1:2129581/100‑1 (MQ=255)
    cTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACa                                                                                              <  1:222066/100‑1 (MQ=255)
    cTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACa                                                                                              <  2:550315/100‑1 (MQ=255)
     tgtttgttTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAg                                                                                             >  1:1648770/1‑100 (MQ=255)
      gtttgtttTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGt                                                                                            >  2:1498376/1‑100 (MQ=255)
      gtttgtttTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGt                                                                                            >  1:1281346/1‑100 (MQ=255)
      gtttgtttTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGt                                                                                            >  1:1625354/1‑100 (MQ=255)
         tgtttTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGAGCATACAGTAta                                                                                         <  1:1234184/100‑2 (MQ=255)
             ttCGTCGATAGCCATTTTTTCTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAgt                                                                                     <  2:145835/100‑1 (MQ=255)
             ttCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAgt                                                                                     <  2:2061715/100‑1 (MQ=255)
             ttCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGAGCATACAGTAtaattg                                                                                     <  1:246946/100‑6 (MQ=39)
             ttCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGAGCATACAGTAtaattg                                                                                     <  2:234337/100‑6 (MQ=39)
             ttCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGAGCATACAGTAtaattg                                                                                     <  2:1821639/100‑6 (MQ=39)
               cgtcgATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGAGCATACAGTAtaattgct                                                                                   <  2:1015381/100‑8 (MQ=38)
                 tcgATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAgtgt                                                                                   >  1:822428/1‑98 (MQ=255)
                 tcgATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGttt                                                                                 >  2:1360775/1‑100 (MQ=255)
                  cgATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGAGCATACAGTAtaattgcttca                                                                                <  1:1016897/100‑11 (MQ=37)
                   gATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTg                                                                               <  2:1769812/100‑1 (MQ=255)
                   gATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTg                                                                               >  2:866799/1‑100 (MQ=255)
                   gATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTg                                                                               >  1:1277566/1‑100 (MQ=255)
                     tAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTa                                                                             <  1:1978357/100‑1 (MQ=255)
                     tAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTa                                                                             <  2:1792548/100‑1 (MQ=255)
                      aGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAg                                                                            >  1:1386181/1‑100 (MQ=255)
                       gCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTa                                                                             >  1:295497/1‑98 (MQ=255)
                       gCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAg                                                                            <  1:481666/99‑1 (MQ=255)
                       gCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAATGTTTTGTAGa                                                                           >  1:214532/1‑100 (MQ=255)
                         cATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGaaa                                                                         <  1:1065908/100‑1 (MQ=255)
                               tACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAAttgttg                                                                   <  1:2004003/100‑1 (MQ=255)
                               tACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAAttgttg                                                                   <  1:1832983/100‑1 (MQ=255)
                               tACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAAttgttg                                                                   <  2:2108054/100‑1 (MQ=255)
                               tACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAAttgttg                                                                   <  1:450377/100‑1 (MQ=255)
                               tACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAAttgttg                                                                   <  1:687251/100‑1 (MQ=255)
                                aCTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGc                                                                  <  2:514733/100‑1 (MQ=255)
                                 cTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGcc                                                                 <  1:1977254/100‑1 (MQ=255)
                                 cTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGcc                                                                 >  2:888743/1‑100 (MQ=255)
                                 cTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGcc                                                                 >  1:1481005/1‑100 (MQ=255)
                                   ccTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCcac                                                               >  1:1237155/1‑100 (MQ=255)
                                    cTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCcaca                                                              >  2:379661/1‑100 (MQ=255)
                                    cTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCcaca                                                              >  2:1464234/1‑100 (MQ=255)
                                       tCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAgg                                                           <  1:119396/100‑1 (MQ=255)
                                       tCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAgg                                                           <  1:271620/100‑1 (MQ=255)
                                             cGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGAAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGATCTGCa                                                     >  2:90192/1‑100 (MQ=255)
                                              gggTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCaa                                                    >  2:766821/1‑100 (MQ=255)
                                              gggTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCaa                                                    <  2:247480/100‑1 (MQ=255)
                                              gggTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCa                                                     <  1:1725620/99‑1 (MQ=255)
                                               ggTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAAt                                                   <  2:1923118/100‑1 (MQ=255)
                                                gTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCa                                                     <  1:577887/97‑1 (MQ=255)
                                                gTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATg                                                  >  2:1524518/1‑100 (MQ=255)
                                                gTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATg                                                  <  1:415911/100‑1 (MQ=255)
                                                  aaTACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCa                                                <  2:31179/100‑1 (MQ=255)
                                                    tACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATa                                              <  1:614654/100‑1 (MQ=255)
                                                    tACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATa                                              <  1:1464244/100‑1 (MQ=255)
                                                    tACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATa                                              >  2:787010/1‑100 (MQ=255)
                                                     aCCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATAc                                             >  2:1607332/1‑100 (MQ=255)
                                                      ccGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACg                                            <  1:1380923/100‑1 (MQ=255)
                                                       cGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACg                                            >  1:94286/1‑99 (MQ=255)
                                                           tAGTCAATATGTTCTTGTGAATTAATTATAGTGTATGCTCATACAGTATCAATTGTTTTTTAGCAATTGTTGCCACAAGGTCTGCAATTCAGACGCAGtt                                       <  2:737997/100‑2 (MQ=255)
                                                             gTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGc                                     >  2:718762/1‑100 (MQ=255)
                                                             gTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGAAGTAg                                      >  2:1031081/1‑99 (MQ=255)
                                                                 atatGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTga                                 <  1:171195/100‑1 (MQ=255)
                                                                 atatGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTga                                 <  1:381667/100‑1 (MQ=255)
                                                                 atatGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTga                                 >  2:283222/1‑100 (MQ=255)
                                                                     gTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTgacgac                             <  1:876209/100‑1 (MQ=255)
                                                                     gTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTgacgac                             >  1:793214/1‑100 (MQ=255)
                                                                     gTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTaacga                              >  1:39642/1‑99 (MQ=255)
                                                                      ttCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGacgacg                            <  2:1237235/100‑1 (MQ=255)
                                                                       tCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGc                           <  1:825945/100‑1 (MQ=255)
                                                                         tGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCa                          <  1:638967/99‑1 (MQ=255)
                                                                         tGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCAc                         >  2:1948158/1‑100 (MQ=255)
                                                                          gTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCAcc                        >  1:1298018/1‑100 (MQ=255)
                                                                          gTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCAc                         >  1:1631417/1‑99 (MQ=255)
                                                                          gTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCAc                         >  2:1395696/1‑99 (MQ=255)
                                                                            tGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGc                      >  1:1098969/1‑100 (MQ=255)
                                                                            tGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGc                      <  1:866800/100‑1 (MQ=255)
                                                                             gAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCa                     <  1:169175/100‑1 (MQ=255)
                                                                               acacttgataCTGTATGAGCATACAGTATCAAGTGTAGATTAGAAAGGGTTCCGACAAATTCTGCAATGAATACGAAGTAGACTGACGACGCAACGCATc                   >  2:1045502/8‑100 (MQ=9)
                                                                               aGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTCGAAATTGTTGCCACAAGGTCTGCAATGCCTACGCAGTAGCCTGACGACGCACCGCATc                   >  1:2208393/1‑100 (MQ=255)
                                                                               aGCAATTATACTGTATGCTCATACAGTATCAAGTGTTCTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATc                   >  2:1196645/1‑100 (MQ=255)
                                                                                cacttgataCTGTATGAGCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCAt                    >  2:1086924/7‑98 (MQ=38)
                                                                                aacttgataCTGTATGAGCATACAGTATCAAGTGTTTTGTAGAAATTGTTGACACAAGGTCTGCAATGCATACGCAGTAGACTGACGACACACCGCATCa                  >  1:1875342/7‑100 (MQ=37)
                                                                                 acttgataCTGTATGAGCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTAACGACAAACCGCATCAc                 >  1:1799317/6‑100 (MQ=38)
                                                                                  cttgataCTGTATGAGCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACg                >  2:7254/5‑100 (MQ=39)
                                                                                  cttgataCTGTATGAGCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACg                >  1:1901343/5‑100 (MQ=39)
                                                                                  cttgataCTGTATGAGCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACg                >  1:1165414/5‑100 (MQ=39)
                                                                                   aTTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACgg               <  1:347264/100‑1 (MQ=255)
                                                                                    tgataCTGTATGAGCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGt              >  1:2233130/3‑100 (MQ=255)
                                                                                    tgataCTGTATGAGCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGt              >  1:426667/3‑100 (MQ=255)
                                                                                      ataCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTcg            >  2:577315/1‑100 (MQ=255)
                                                                                        aCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTcgcc          <  2:1480995/100‑1 (MQ=255)
                                                                                              tGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTg      >  1:916757/1‑98 (MQ=255)
                                                                                               gCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAg   <  1:463327/100‑1 (MQ=255)
                                                                                               gCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAg   >  1:985232/1‑100 (MQ=255)
                                                                                                cTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGc  <  1:1963018/100‑1 (MQ=255)
                                                                                                 |                                                                                                  
CTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGC  >  NC_000913/2823742‑2823937

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