Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,532,409 T→G intergenic (‑758/+28) sgcX ← / ← yjhP putative endoglucanase with Zn‑dependent exopeptidase domain/putative methyltransferase YjhP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,532,4090TG98.6% 261.5 / NA 71intergenic (‑758/+28)sgcX/yjhPputative endoglucanase with Zn‑dependent exopeptidase domain/putative methyltransferase YjhP
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (49/21);  total (49/22)
Fisher's exact test for biased strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCG  >  NC_000913/4532170‑4532653
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    
aCCTGATAACGCTTCTTCCCAGGGTTTGAGATTGCCGTTTCTGTCGCTGCATTACGCACCATATTGATCAATGGCCCTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGAAATCTTGTTTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGAACAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGaa                                                                                                                                                                                                                                            >  2:44706/1‑250 (MQ=255)
        ccGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCggg                                                                                                                                                                                                                                    >  2:299374/1‑250 (MQ=255)
             tcagcccAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTCCGCACAATATTGATCACTGGCACTGTGAGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCCTGTTGAGACATATAATGCTCCCTATGGTGACATCTTGATTAATAAACATCCTGAATATGTATTGTGCAGGTGATACATCAGTTCAAAGCAAAGGTAAACCCTCGATATTGCCAGCCTGAATTACAGGGAACCtg                                                                                                                                                                                                                              >  2:336028/5‑251 (MQ=255)
                   cAGGGTTTGAAATTGCAGTTGCTGTCGCTGCATTACGCGCAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGCCAGCACTGTTCATGTTGTGACATATGATGCTACTTATGGTGACATCTTGATTAATAAACATACTGCATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGAGACCTgttgtt                                                                                                                                                                                                                         >  2:184299/1‑250 (MQ=255)
                    aGGGTTTGAGATTGCCGTTGCTGTAGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGAAAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGAATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTAc                                                                                                                                                                                                                       >  2:71330/1‑251 (MQ=255)
                     gggTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTAcc                                                                                                                                                                                                                      <  1:60045/251‑1 (MQ=255)
                      ggTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCg                                                                                                                                                                                                                     <  1:17651/251‑1 (MQ=255)
                                 gCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTg                                                                                                                                                                                                          >  2:101614/1‑251 (MQ=255)
                                  ccGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGCTCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGc                                                                                                                                                                                                         >  2:339638/1‑251 (MQ=255)
                                       gCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGACTATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAac                                                                                                                                                                                                     >  1:197141/1‑250 (MQ=255)
                                          gTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAacac                                                                                                                                                                                                   <  2:335445/249‑1 (MQ=255)
                                             gCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGcccc                                                                                                                                                                                              >  2:88499/1‑251 (MQ=255)
                                                 cATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGcc                                                                                                                                                                                          >  1:337459/1‑251 (MQ=255)
                                                              tttGATCCCTGGCACTGTGTGCCTCTGGTGGCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTCCCCTATGATGCTCCTTATGGTGACATCTTGCTTACTCAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCAAGCCAAAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCACCCAAAGCATTCCcgc                                                                                                                                                                             >  2:176745/2‑251 (MQ=255)
                                                              atTGATCACTGGCACTGTGTGCATCTGGAGCAGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTCTGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGACATACCGCGACCCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCACAAGCATTACcgc                                                                                                                                                                             >  2:94977/1‑251 (MQ=255)
                                                                tGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCcgcgc                                                                                                                                                                           >  1:391905/1‑251 (MQ=255)
                                                                 gATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCcgcgcg                                                                                                                                                                          >  1:82736/1‑251 (MQ=255)
                                                                        ggCACTGTGTGCATCTGGTGCCGACATCGGCATTCAGCAAGCCTGACAGAAATGTTCATGTTGTGACATATGATGCTCCTTATGGTGCCATCTTGATTACTAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGACCACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTcc                                                                                                                                                                   >  2:203897/1‑250 (MQ=255)
                                                                                           gCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGAt                                                                                                                                                >  1:6192/1‑251 (MQ=255)
                                                                                                 aTCGGCATTCCGGAAGCCTGACAGAAATGTTCATGTTGTGACATAAGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGAGACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTTTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTTGCGTTTCGGCGTCATGTTTAg                                                                                                                                          <  2:194011/251‑1 (MQ=255)
                                                                                                          ccGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCt                                                                                                                                 >  2:74347/1‑251 (MQ=255)
                                                                                                              aaGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACt                                                                                                                             >  1:172275/1‑251 (MQ=255)
                                                                                                                                cATGTTGTTACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATATATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCg                                                                                                           <  2:244010/251‑1 (MQ=255)
                                                                                                                                cATGTTGTGACATAAGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGAAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGGAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCg                                                                                                           >  2:287319/1‑251 (MQ=255)
                                                                                                                                   gttgtGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTc                                                                                                         <  1:66096/250‑1 (MQ=255)
                                                                                                                                     tgtgACATATGATGCTCCTTATGGTGACATCTTGCTTAATAAACATACTGAATCTGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAgg                                                                                                      >  1:123800/1‑251 (MQ=255)
                                                                                                                                          cATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGTCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCg                                                                                                 <  2:243990/251‑1 (MQ=255)
                                                                                                                                            tatGATGCTCCTTATGGTGACATATTGATTAATAAACATATTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGACGTCATCGTCAGGATTCGCt                                                                                               >  2:31748/1‑251 (MQ=255)
                                                                                                                                             atgatgCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGACTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGACACGAAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCtt                                                                                              >  2:342306/1‑251 (MQ=255)
                                                                                                                                               gatgCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCAAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCAGCGGCGAAGTCATCGTCAGGATTCGCTTcc                                                                                            >  1:168333/1‑251 (MQ=255)
                                                                                                                                                atgCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCCTTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTccc                                                                                           >  2:401497/1‑250 (MQ=255)
                                                                                                                                                     ccTTATGGTGAAATCTTGATTAATAAACATCCTGCCTACGTATTGTGCAGGAGATACATCAGATCAAAGCACAGGAACAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGATACCGGGCGATGAGAGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCCTCGTCAGGATTCGCTTccagaca                                                                                       >  2:175889/1‑247 (MQ=255)
                                                                                                                                                                    ttGATTAATAAACATACTGAATATGTATTGTGCGGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGc                                                                       >  2:274014/1‑251 (MQ=255)
                                                                                                                                                                    ttGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATAAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGc                                                                       >  2:407335/1‑251 (MQ=255)
                                                                                                                                                                     tGATTAATAAACAAACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGTTATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGc                                                                       <  2:92210/250‑1 (MQ=255)
                                                                                                                                                                      gATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGcc                                                                      >  1:69317/1‑250 (MQ=255)
                                                                                                                                                                         taataaACATACTGAATAGGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGTCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGGCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGTTCAGCCAtt                                                                   <  1:165973/250‑1 (MQ=255)
                                                                                                                                                                          aataaACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCAtttt                                                                 >  2:100713/1‑251 (MQ=255)
                                                                                                                                                                          aataaACATACTGAATATGAATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATAACGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGAATTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGACAtttt                                                                 >  1:329750/1‑251 (MQ=255)
                                                                                                                                                                             aaaCATACTGAATATGTATTGTGCGGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGc                                                               <  1:274014/250‑1 (MQ=255)
                                                                                                                                                                               aCATACTGAATATGTATTGTGCAGGTGATCCATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGCTGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGc                                                            >  2:82191/1‑251 (MQ=255)
                                                                                                                                                                                cATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCAAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGACCACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGAGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGACGTCATCGACAGGATTCGCTTCCAGCCAGCGGCGACTGGTCAGCCATTTTGCGGCt                                                           >  1:410024/1‑251 (MQ=255)
                                                                                                                                                                                  tACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTc                                                         >  2:132354/1‑251 (MQ=255)
                                                                                                                                                                                    cTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCCCGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCAt                                                       >  2:271044/1‑251 (MQ=255)
                                                                                                                                                                                      gAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATAc                                                     >  2:90226/1‑251 (MQ=255)
                                                                                                                                                                                                tgtgCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACAGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGAAGCTCACGTAGCGTTTCGGCGCGCTGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGCTTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCATCATACCGGTCCCAGc                                           >  1:219151/1‑251 (MQ=255)
                                                                                                                                                                                                tgtgCAGGTGATACATCAGACCAAAGCAAAGGTAAAGCCTCGATAATGCCAGCCAGAAATACCGGGACCCTGTTGTTACCGGGCGATGAGTGCCAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAg                                            <  2:62410/250‑1 (MQ=255)
                                                                                                                                                                                                   gCAGGTGCTACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGccc                                                                                                                                  >  2:212990/1‑161 (MQ=255)
                                                                                                                                                                                                   gCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGccc                                                                                                                                  <  1:212990/161‑1 (MQ=255)
                                                                                                                                                                                                      ggTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATCCCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCt                                     <  2:337080/251‑1 (MQ=255)
                                                                                                                                                                                                         gATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCCGCAATACCGGGAACCTGTTGTTACCGAGCGCTGAGCGAGAACAAGCCACAGCCCAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGACCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCAAGCCAGCGACGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCCGCCTTCCTGGt                                  >  1:31096/1‑251 (MQ=255)
                                                                                                                                                                                                          aTACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTc                                 <  2:117458/251‑1 (MQ=255)
                                                                                                                                                                                                             cATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCAGCGAAGTCATCGTCAGGCTTCGCTTCCAGCCAGCGGCGCCTGGTCAGCCATTTTACGGCTTCATACCGGTCCCCGCCTTCCTGGTCtcc                              >  2:319857/1‑249 (MQ=255)
                                                                                                                                                                                                                cAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATgg                                                                            >  2:8588/1‑202 (MQ=255)
                                                                                                                                                                                                                cAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATgg                                                                            <  1:8588/202‑1 (MQ=255)
                                                                                                                                                                                                                cAGATCAAAGCAAAGGTAAAGACTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAACCTTCCTgtgccgcc                             >  2:286675/1‑242 (MQ=255)
                                                                                                                                                                                                                   aTCAAAGCAAAGGAAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTAc                        >  1:132278/1‑251 (MQ=255)
                                                                                                                                                                                                                     cacagcaacgGTAAAGCCTCGATATTGCCAGCCAGAACTACCGGGCACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCACAGCCAAAGAATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCAAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGCCCCAGCCTTCCTGGTCTGCCAGTAccc                      >  2:246594/4‑250 (MQ=255)
                                                                                                                                                                                                                       aagcaacgGTAAAGCCTCGATATTGCCAGACAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAgacccatt                    >  2:399332/1‑244 (MQ=255)
                                                                                                                                                                                                                       aagcaaagGTAAAGCCTCGATATTGACAGCCAGAAATACCGGGAACCTGTTGTAACCGGGCGAGGAGTGCGAACACGCCCCAGCCATAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCAtt                    <  2:217524/251‑1 (MQ=255)
                                                                                                                                                                                                                        agcaaagGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCAttt                   <  2:383415/251‑1 (MQ=255)
                                                                                                                                                                                                                         gcaaagGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTc                  >  1:71443/1‑251 (MQ=255)
                                                                                                                                                                                                                            aagGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTTTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTccac               <  1:193911/251‑1 (MQ=255)
                                                                                                                                                                                                                            aagGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATAACGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTccac               >  1:237497/1‑251 (MQ=255)
                                                                                                                                                                                                                              ggTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTccacca             <  2:337459/251‑1 (MQ=255)
                                                                                                                                                                                                                               gTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCcaccac            <  1:339638/251‑1 (MQ=255)
                                                                                                                                                                                                                                        cGATATTGCCAGCCATAAATTCCGGGAACCTGTTGTTCCCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCCCGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCTCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGcc   <  2:305770/251‑1 (MQ=255)
                                                                                                                                                                                                                                        cGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTTCCAGTACCATTTCCACCACGTCATAGcc   >  2:156866/1‑251 (MQ=255)
                                                                                                                                                                                                                                         gATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCa                                             >  2:406211/1‑208 (MQ=255)
                                                                                                                                                                                                                                         gATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCa                                             <  1:406211/208‑1 (MQ=255)
                                                                                                                                                                                                                                          atatTGCCCGCCAGAAATAACGGGCACCTGTTGTTACCGGGCGATGCGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCCGTACCATTTCCACCACGTCATAGCCg  >  2:43528/1‑250 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    
ACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCG  >  NC_000913/4532170‑4532653

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