Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,657,898 A→T 6.2% intergenic (+28/‑171) ynfD → / → ynfE DUF1161 family periplasmic protein/putative selenate reductase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,657,8980AT 6.2% 446.7 / 15.3 130intergenic (+28/‑171)ynfD/ynfEDUF1161 family periplasmic protein/putative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (56/66);  new base T (3/5);  total (59/71)
Fisher's exact test for biased strand distribution p-value = 7.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

Only 100 of 138 total aligned reads displayed.

CCAGGTCGTCGGCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAAC  >  NC_000913/1657750‑1658046
                                                                                                                                                    |                                                                                                                                                    
ccAGGTCGTCGGCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGa                                                                                                                                                      >  2:567314/1‑149 (MQ=255)
ccAGGTCGTCGGCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGa                                                                                                                                                      <  2:471273/149‑1 (MQ=255)
ccAGGTCGTCGGCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGa                                                                                                                                                      >  3:818411/1‑149 (MQ=255)
ccAGGTCGTCGGCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGa                                                                                                                                                      >  1:430615/1‑149 (MQ=255)
  aGGTCGTCGGCCATTGCGCTAATGATACGCATAATATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACa                                                                                                                                                    >  1:416871/1‑149 (MQ=255)
  aGGTCGTCGGCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACa                                                                                                                                                    >  2:839243/1‑149 (MQ=255)
       gtcgGCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATcc                                                                                                                                               <  3:602636/149‑1 (MQ=255)
         cgGCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCgg                                                                                                                                             >  3:356524/1‑149 (MQ=255)
          ggCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGc                                                                                                                                            <  2:204048/149‑1 (MQ=255)
            ccATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgt                                                                                                                                          <  1:296138/149‑2 (MQ=255)
             cATTGCACTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGAt                                                                                                                                         <  1:148362/149‑1 (MQ=255)
              aTTGCGCTAATGATACGCAGAAAATTATCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCCCGCAGTCTAGCCAGGATGCCGCGCCTGCCGAGCCGCAGTAACACATTATCGACTGAACGCCGGCTATGACAAATCCGGCGAtt                                                                                                                                        >  2:73442/1‑149 (MQ=255)
               ttGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgttca                                                                                                                                       <  4:648206/149‑5 (MQ=255)
               ttGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGAttt                                                                                                                                       <  3:239523/149‑1 (MQ=255)
                  cgcTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGaa                                                                                                                                    <  2:35171/149‑1 (MQ=255)
                  cgcTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGaa                                                                                                                                    <  4:171605/149‑1 (MQ=255)
                   gcTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGa                                                                                                                                                      >  4:110258/1‑130 (MQ=255)
                   gcTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAAc                                                                                                                                   <  4:443698/149‑1 (MQ=255)
                    cTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACa                                                                                                                                  >  2:569524/1‑149 (MQ=255)
                    cTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACa                                                                                                                                  >  1:878500/1‑149 (MQ=255)
                     tAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACAt                                                                                                                                 >  3:24932/1‑149 (MQ=255)
                         gATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaa                                                                                                                             <  3:517655/149‑1 (MQ=255)
                          aTACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaac                                                                                                                            >  1:395539/1‑149 (MQ=255)
                           tACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaaca                                                                                                                           <  1:839243/149‑1 (MQ=255)
                             cGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATa                                                                                                                         <  3:122742/149‑1 (MQ=255)
                               cATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAAt                                                                                                                       >  3:659040/1‑149 (MQ=255)
                                aTAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATc                                                                                                                      <  2:980175/149‑1 (MQ=255)
                                aTAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATc                                                                                                                      >  1:1065216/1‑149 (MQ=255)
                                 tAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATcc                                                                                                                     >  4:833496/1‑149 (MQ=255)
                                  aaaaTTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATccc                                                                                                                    <  1:682171/149‑1 (MQ=255)
                                         tctATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCttat                                                                                                             <  1:250175/149‑1 (MQ=255)
                                           tataCCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCttatta                                                                                                           >  1:674768/1‑149 (MQ=255)
                                              aCCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTc                                                                                                        >  2:1037344/1‑149 (MQ=255)
                                               cccGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCGGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATACCACCTTATTACTCa                                                                                                       >  3:734640/1‑149 (MQ=255)
                                                 cGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCGCCTTATTACTCATa                                                                                                     <  2:175675/149‑1 (MQ=255)
                                                 cGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATa                                                                                                     <  4:970858/149‑1 (MQ=255)
                                                 cGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATa                                                                                                     <  1:497672/149‑1 (MQ=255)
                                                    accaccAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAAc                                                                                                                                   >  4:212229/1‑116 (MQ=255)
                                                    accaccAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAAc                                                                                                                                   <  3:212229/116‑1 (MQ=255)
                                                       accaGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAAGATAATCCCACCTTATTACTCATACCCTTc                                                                                               <  3:877976/149‑1 (MQ=255)
                                                        ccaGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCt                                                                                              <  2:371470/149‑1 (MQ=255)
                                                        ccaGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCt                                                                                              <  3:1000882/149‑1 (MQ=255)
                                                         caGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTa                                                                                             >  4:182317/1‑149 (MQ=255)
                                                              gTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCg                                                                                                                                              <  2:536089/95‑1 (MQ=255)
                                                                aaCGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGatat                                                                                      <  1:536037/149‑1 (MQ=255)
                                                                aaCGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGatat                                                                                      <  4:624030/149‑1 (MQ=255)
                                                                 aCGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATg                                                                                     <  1:1024965/149‑1 (MQ=255)
                                                                   gTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGa                                                                                   >  2:893969/1‑149 (MQ=255)
                                                                   gTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGa                                                                                   <  4:417956/149‑1 (MQ=255)
                                                                   gTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGa                                                                                   <  3:270657/149‑1 (MQ=255)
                                                                    tctcTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACGTTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGAt                                                                                  >  1:568171/1‑149 (MQ=255)
                                                                    tctcTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTAtt                                                                                           >  2:816397/1‑140 (MQ=255)
                                                                    tctcTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCACAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTAtt                                                                                           <  1:816397/140‑1 (MQ=255)
                                                                        tGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATAACAAATCCTGCAATTTGAACATAGAAAATAATCCTTCCTTATTACTCATACCCTTCTATTGATATGAATtaat                                                                              <  2:802368/149‑1 (MQ=255)
                                                                         gCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATtaata                                                                             >  3:368203/1‑149 (MQ=255)
                                                                           tCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGata                                                                                       >  3:587883/1‑137 (MQ=255)
                                                                           tCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGata                                                                                       <  4:587883/137‑1 (MQ=255)
                                                                           tCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTaataat                                                                           >  4:627605/1‑149 (MQ=255)
                                                                              cgcgCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCt                                                                        <  1:888281/149‑1 (MQ=255)
                                                                                 gcAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTaa                                                                     <  3:992069/149‑1 (MQ=255)
                                                                                   aGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAAcc                                                                   >  1:1035540/1‑149 (MQ=255)
                                                                                    gTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAAccc                                                                  >  4:67713/1‑149 (MQ=255)
                                                                                     tCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCa                                                                 >  2:700591/1‑149 (MQ=255)
                                                                                       tAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCt                                                                                              >  4:716763/1‑118 (MQ=255)
                                                                                       tAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCt                                                                                              <  3:716763/118‑1 (MQ=255)
                                                                                         gCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAAt                                                             <  2:1092663/149‑1 (MQ=255)
                                                                                          ccAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATg                                                            >  4:481205/1‑149 (MQ=255)
                                                                                          ccAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATg                                                            >  4:27441/1‑149 (MQ=255)
                                                                                           cAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATgg                                                           <  3:120966/149‑1 (MQ=255)
                                                                                           cAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATgg                                                           <  2:1035540/149‑1 (MQ=255)
                                                                                           cAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATgg                                                           <  2:674768/149‑1 (MQ=255)
                                                                                           cAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATgg                                                           <  4:78894/149‑1 (MQ=255)
                                                                                            aGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATac                                                                                                                               <  2:328462/80‑1 (MQ=255)
                                                                                            aGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATggg                                                          >  1:161044/1‑149 (MQ=255)
                                                                                             ggATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGt                                                         <  3:18800/149‑1 (MQ=255)
                                                                                                     cgcCTGCCTAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTccc                                                 <  4:1062654/149‑1 (MQ=255)
                                                                                                      gcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCct                                                >  1:1011554/1‑149 (MQ=255)
                                                                                                         tGCCGAACCGCAGAAACACATTATCGACTGAACGCCGTATATGACAAATTCAGCCACTAGAACAAACAACATAATACCACCTTATTACTCATTCCCTTCTATTGATGGGGTTTAATATTTCTTAACCCAAAATTAATAAAATCCCTCTa                                             >  2:275569/1‑149 (MQ=255)
                                                                                                           ccGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTAtt                                           <  1:567314/149‑1 (MQ=255)
                                                                                                            cGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAATATACAACATAATCCCACCTTATTACTc                                                                                                        >  2:436520/1‑87 (MQ=255)
                                                                                                              aaCCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGtt                                        >  1:820651/1‑149 (MQ=255)
                                                                                                               aCCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTa                                       <  4:554739/149‑1 (MQ=255)
                                                                                                                 cGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAgc                                     <  4:875461/149‑1 (MQ=255)
                                                                                                                      tAACACATTATCGAATGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCt                                >  3:84040/1‑149 (MQ=255)
                                                                                                                          acaTTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAat                            >  3:552295/1‑149 (MQ=255)
                                                                                                                             ttATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATAtt                         <  4:958618/149‑1 (MQ=255)
                                                                                                                               aTCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTaaa                             <  4:173789/143‑1 (MQ=255)
                                                                                                                               aTCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTaaa                             >  3:173789/1‑143 (MQ=255)
                                                                                                                               aTCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCa                       >  3:953380/1‑149 (MQ=255)
                                                                                                                                  gACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAata                    <  2:883142/149‑1 (MQ=255)
                                                                                                                                  gACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAata                    >  1:437832/1‑149 (MQ=255)
                                                                                                                                   aCTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAatat                   >  2:224631/1‑149 (MQ=255)
                                                                                                                                       aaCGCCGGATACGATAAATCCAGCGATTTGAACCCACAACATAATCCCGTCATATTACACATACCCTTCTATTTATATGGATTAATAATTCTTGATCCAAAATGGGTAGACTCCATCTATTGTTAGGGCGCTAAATATACGATATATaa               <  2:177684/149‑1 (MQ=255)
                                                                                                                                              gATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTta        <  3:367792/149‑1 (MQ=255)
                                                                                                                                              gATATGACAAATACGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGTTTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGATAAATATTCAATATATATACTTTta        >  4:276504/1‑149 (MQ=255)
                                                                                                                                                   gACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTAtt                                           <  2:664780/109‑1 (MQ=255)
                                                                                                                                                   gACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTAtt                                           >  1:664780/1‑109 (MQ=255)
                                                                                                                                                   gACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATACGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCaa                      <  3:991323/130‑1 (MQ=255)
                                                                                                                                                    aCAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAAc  >  2:541601/1‑149 (MQ=255)
                                                                                                                                                    aCAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAAc  >  2:972740/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCAGGTCGTCGGCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAAC  >  NC_000913/1657750‑1658046

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