Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,988,115 C→T D318N (GAT→AAT)  hemX ← PF04375 family protein HemX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,988,1150CT98.2% 197.2 / NA 54D318N (GAT→AAT) hemXPF04375 family protein HemX
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (26/27);  minor base G (0/1);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTTTTAGCATGGTTATTCTCCTTGCGGCGCGGCTGCCGGAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCATAAAGTTCTGCCAGCTTTTTTGCAGATTGATACGCCATTCGCTGATGGAACT  >  NC_000913/3987871‑3988358
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
cacTTTTAGCATGGTTATTCTCCTTCCGGCGCGGCTGCCGGAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTGGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGGCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTgtagta                                                                                                                                                                                                                                               <  1:337707/251‑1 (MQ=255)
    tttAGCATGGTTATTCTCCTTGCGGCGCGGCTGCCGGAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGc                                                                                                                                                                                                                                            >  1:57183/1‑250 (MQ=255)
            ggTTATTCTCCTTGCGGCGCGGCTGCCGGAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTAccc                                                                                                                                                                                                                                    >  2:347970/1‑250 (MQ=255)
                     ccTTGCGGCGCGGCTGCCGGAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGAc                                                                                                                                                                                                                          >  1:13018/1‑251 (MQ=255)
                          cggcgcggcTGCCGGAGTATCAGCTTGCGGTGCAGGGGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACTCACGCTAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCGCGTACCCAGGTGGAGACGtgct                                                                                                                                                                                                                     <  2:192821/251‑4 (MQ=255)
                             cgcggcTGCCGGAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACg                                                                                                                                                                                                                         >  2:207841/1‑244 (MQ=255)
                             cgcggcTGCCGGAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACg                                                                                                                                                                                                                         <  1:207841/244‑1 (MQ=255)
                                       ggAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCg                                                                                                                                                                                                        >  2:308900/1‑251 (MQ=255)
                                              cAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGt                                                                                                                                                                                                  <  1:53717/250‑1 (MQ=255)
                                                     cGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGCACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTg                                                                                                                                                                                         >  2:145335/1‑251 (MQ=255)
                                                                  ggtTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCCGCAGGTTACGCACGCGAGTCTGCCTCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGGGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGAGATAAGTCTCTTCCTGGTGACGCGGTa                                                                                                                                                                              >  1:343081/1‑250 (MQ=255)
                                                                         cTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTtgt                                                                                                                                                                      >  2:92281/1‑251 (MQ=255)
                                                                                         tGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGc                                                                                                                                                      >  2:122756/1‑251 (MQ=255)
                                                                                            cGGTGGTGCCATCAGGGTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCGGGCTTTGCAGGGGTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACTTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGGACAGCTTGTGCTGCGACCAGCAGGCgag                                                                                                                                                   <  2:274987/251‑1 (MQ=255)
                                                                                                     cagcagGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAatat                                                                                                                                          <  2:311204/251‑1 (MQ=255)
                                                                                                     cagcagGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAatat                                                                                                                                          >  1:99438/1‑251 (MQ=255)
                                                                                                         agGTTACGCGCGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTc                                                                                                                                      <  2:3251/251‑1 (MQ=255)
                                                                                                         agGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTc                                                                                                                                      >  1:153710/1‑251 (MQ=255)
                                                                                                             tACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTcgcg                                                                                                                                   <  2:294395/250‑1 (MQ=255)
                                                                                                             tACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTcgcg                                                                                                                                   <  2:57183/250‑1 (MQ=255)
                                                                                                                  acgcGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCagata                                                                                                                             <  1:206015/251‑1 (MQ=255)
                                                                                                                        gTCTGCATCAGTTTTTCCGGCATCGCTTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGCTAGATATc                                                                                                                       >  2:41471/1‑251 (MQ=255)
                                                                                                                            gCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCt                                                                                                                     <  1:360140/249‑1 (MQ=255)
                                                                                                                            gCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTGGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATcc                                                                                                                      <  2:205105/248‑1 (MQ=255)
                                                                                                                                cAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCATGATTTg                                                                                                               >  2:54412/1‑251 (MQ=255)
                                                                                                                                         cAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCATGATTTGGCGCTAACa                                                                                                      <  1:54412/251‑1 (MQ=255)
                                                                                                                                             aTCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCgg                                                                                                  <  1:278840/251‑1 (MQ=255)
                                                                                                                                                gCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTAc                                                                                               >  2:287973/1‑251 (MQ=255)
                                                                                                                                                 ccTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTAcg                                                                                              <  2:8236/251‑1 (MQ=255)
                                                                                                                                                     gCTTTGCAGGGTTTCCGGACGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCCGCAGGCGAGAGCGAATATTTTCGAGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTAcggcgg                                                                                          >  2:147016/1‑251 (MQ=255)
                                                                                                                                                         tGCAGGGCTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGtc                                                                                      >  1:315891/1‑251 (MQ=255)
                                                                                                                                                          gCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGtc                                                                                      <  2:23868/250‑1 (MQ=255)
                                                                                                                                                              ggTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGGCATCac                                                                                 >  1:307556/1‑251 (MQ=255)
                                                                                                                                                                   ccGGAAGATGCATCGAGATATTTTGCTGGCTTAACTGTTCCACCTCGGCGAGGAACGCTTTGGTGGGGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGTCGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCATGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCacgacg                                                                             <  2:226757/250‑1 (MQ=255)
                                                                                                                                                                   ccGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACgg                                                                            <  1:308900/251‑1 (MQ=255)
                                                                                                                                                                     ggAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGGCATCACGCCGGCg                                                                          >  2:278662/1‑251 (MQ=255)
                                                                                                                                                                               tCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATg                                                                >  1:313937/1‑251 (MQ=255)
                                                                                                                                                                                 gagaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGa                                                               >  1:177004/1‑250 (MQ=255)
                                                                                                                                                                                            cTGGCTTAACTGGTCCACATCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCCTCACGCCGGCGGATCGGCATGAAATTGTCCATa                                                    >  2:162673/1‑250 (MQ=255)
                                                                                                                                                                                              ggCTTTACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCGCGTTCCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTTCTGCGGCCAGCAGGTGAGAGCGAATATTTTCGCGCAGGTAGATATCCTGAGTTGGCGGTAACAGCAGTACGGCGGGGTGTTCACGACGGTGAATCTTATTGAAGTGGTCCATaaa                                                  <  2:211032/250‑1 (MQ=255)
                                                                                                                                                                                              ggCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCAGCACGACGGCGAATCGTAATGAAGTTGTCCATaaa                                                  <  2:238192/250‑1 (MQ=255)
                                                                                                                                                                                                    aCTGGTCCACCTCGTCGAGGCACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGGATCGTCATGAAGTTGTCCATAAAGTTCTGc                                           >  1:267162/1‑251 (MQ=255)
                                                                                                                                                                                                    aCTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCATAAAGTTCTGc                                           <  1:209164/251‑1 (MQ=255)
                                                                                                                                                                                                     cTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGGGAATCGTAATGAAGTTGTCCCTAAAGTTCttcc                                          >  2:117515/1‑248 (MQ=255)
                                                                                                                                                                                                       ggTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCATAAAGTTCTGCCAg                                        <  1:117515/251‑1 (MQ=255)
                                                                                                                                                                                                            cccTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCATAAAGTTCTGCCAGCtttt                                   <  2:204965/250‑1 (MQ=255)
                                                                                                                                                                                                            aCCTCGTCGAGGAACGCTTTGGTGGAGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGCCGGCGAATCGTAATGAAGTTGTCCATAAAGTTCTGCCAGCtttt                                   >  1:300052/1‑251 (MQ=255)
                                                                                                                                                                                                                   cgAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTGCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAAGGAAGTGGTCCATAAAGTTCTGCCAGCTTTTTTGCAg                             <  2:161372/250‑1 (MQ=255)
                                                                                                                                                                                                                            cTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGGTGTCCATAAAGTTCTGCCCGCTTTTTTGCAGATTGATACGc                   >  1:176334/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gtggtggCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTGGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCATAAAGTTCTGCCAGCTTTTTTGCAGATTGATACGCCATTc              <  2:231214/251‑1 (MQ=255)
                                                                                                                                                                                                                                        catcatCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCATAAAGTTCTGCCAGCTTTTTTGCAGATTGATACGCCATTCGCTGATg       <  1:322371/251‑1 (MQ=255)
                                                                                                                                                                                                                                         atcatcAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCATAAAGTTCTGCCAGCTTTTTTGCAGATTTATACGCCATTCGCTGATgg      >  1:294221/1‑251 (MQ=255)
                                                                                                                                                                                                                                              caGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCATAAAGTGCTGCCAGCTTTTTTGCAGATTGATACGCCATTCGCTGATGGAAcg  <  2:116800/250‑2 (MQ=255)
                                                                                                                                                                                                                                              caGTATGGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGAAGGCGAGAGCTAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCATAAAGTTCTGCCAGCTTTTTTGCAGATTGATACGCCATTCGCTGACGGAAcg  <  2:329423/250‑2 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
CACTTTTAGCATGGTTATTCTCCTTGCGGCGCGGCTGCCGGAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCATAAAGTTCTGCCAGCTTTTTTGCAGATTGATACGCCATTCGCTGATGGAACT  >  NC_000913/3987871‑3988358

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