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,1150CT95.6% 219.0 / 1.3 68D318N (GAT→AAT) hemXPF04375 family protein HemX
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (29/36);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

TAGCATGGTTATTCTCCTTGCGGCGCGGCTGCCGGAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCATAAAGTTCTGCCAGCTTTTTTGCAGATTGATACGCCATTCGCTGATGGAA  >  NC_000913/3987877‑3988356
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    
tAGCAGGGTGATTCTCCTTGCGGCGCGGCTGCCGGAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCAc                                                                                                                                                                                                                                           <  2:553292/250‑1 (MQ=255)
           ttCTCCTTGCGGCGCGGCTGCCGGAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCCCGTACCCAGGTgg                                                                                                                                                                                                                               >  2:102996/1‑251 (MQ=255)
                  tGCGGCGCGGCTGCCGGAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACg                                                                                                                                                                                                                          <  1:76713/249‑1 (MQ=255)
                  tGCGGCGCGGCTGCCGGAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGtt                                                                                                                                                                                                                        >  1:332442/1‑251 (MQ=255)
                            cTGCCGGAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGcc                                                                                                                                                                                                              >  1:499965/1‑251 (MQ=255)
                                cGGAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGc                                                                                                                                                                                                          <  1:69764/251‑1 (MQ=255)
                                        cAGCTTGCGGGCCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGGCCACCTCGTCGAGGAACGCTTTGTTGATGGCATCATCCGTATTGTAGTAAGCACGTACCCAGGAGGAGAAGTTCTCCAGCGCCTGGCGATAAGtc                                                                                                                                                                                                  >  2:506748/1‑251 (MQ=255)
                                           cTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTctct                                                                                                                                                                                               >  1:536589/1‑251 (MQ=255)
                                                ggtgcaggtgcaGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTg                                                                                                                                                                                          >  2:179653/1‑251 (MQ=255)
                                                     aggtgcaggtTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGAc                                                                                                                                                                                     >  2:331706/1‑251 (MQ=255)
                                                       gtgcaggtTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTACGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGAcgc                                                                                                                                                                                   >  2:427305/1‑251 (MQ=255)
                                                             gtTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGAc                                                                                                                                                                                                                           >  1:523650/1‑205 (MQ=255)
                                                             gtTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGAc                                                                                                                                                                                                                           <  2:523650/205‑1 (MQ=255)
                                                              tttAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTAc                                                                                                                                                                              >  1:202317/1‑249 (MQ=255)
                                                                 aGCTTCCGTTGTCCCCGCGGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCtt                                                                                                                                                                         <  2:442819/251‑1 (MQ=255)
                                                                   cTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTtgt                                                                                                                                                                       >  2:238634/1‑251 (MQ=255)
                                                                       cGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGtgctg                                                                                                                                                                  <  1:303594/251‑1 (MQ=255)
                                                                        gTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGtgctgc                                                                                                                                                                 <  1:129573/251‑1 (MQ=255)
                                                                             ccccGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGAcc                                                                                                                                                             >  1:36437/1‑250 (MQ=255)
                                                                             ccccGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGAcc                                                                                                                                                             <  2:348662/250‑1 (MQ=255)
                                                                                      cgggtggtccAGCAGGGTACGCACGCGAGTCTGCATCAGTGTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCgag                                                                                                                                                   <  2:170267/243‑1 (MQ=255)
                                                                                      cGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCgag                                                                                                                                                   <  1:379195/251‑1 (MQ=255)
                                                                                              ccAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAat                                                                                                                                            >  1:24183/1‑250 (MQ=255)
                                                                                                agcagGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGATGCGACCAGCAGGCGAGAGCGAATAtt                                                                                                                                         >  2:386756/1‑251 (MQ=255)
                                                                                                  cagGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATAtttt                                                                                                                                       <  1:285584/251‑1 (MQ=255)
                                                                                                  cagGATACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGAACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTATTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATAtttt                                                                                                                                       >  2:499584/1‑251 (MQ=255)
                                                                                                         cgctggcgAGTCTGCACCAGTTTTCCCAGCAGCGCCTGCCGTTGCCGGTGTCCCGTAATATCCATCGAGTGATGTTGCTGGCTTAACGTGTTCACCTCGTCGAGGATTGCTTTGGTGGTGGCCTCATCAGTATTGTAGTAATCACGTACCCAGGTGTAGACGTTCTCCAGTGCCTGGCGATAAGTCTCTGTCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGATCGAATATTTTCGCGCag                                                                                                                                <  2:462003/246‑1 (MQ=255)
                                                                                                              gcgAGTCTGCCTAAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGCAGCTCCATCGCGATCTTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGACCGCTTTGGTGGTGGCATCCTCAGTATTGTAGTAAGCACGTACCAAGGCGGAGACGTTCTCCAGCGCCTGCCGTTAAGTCTCTTCCTGGTGACGCGGTACAACTTTGTCTGCGACCCGCAGGCGCGAGCGCTTATTTTCGCGCagttagc                                                                                                                           >  2:241988/1‑250 (MQ=255)
                                                                                                                          cAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATAATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTg                                                                                                               >  1:40484/1‑251 (MQ=255)
                                                                                                                           aGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTgg                                                                                                              <  1:209049/251‑1 (MQ=255)
                                                                                                                                  ccAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTCTTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGTTTTGGCGCTAAc                                                                                                       >  2:554220/1‑251 (MQ=255)
                                                                                                                                     gCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGc                                                                                                    >  1:206046/1‑251 (MQ=255)
                                                                                                                                      cATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCg                                                                                                   >  1:476328/1‑251 (MQ=255)
                                                                                                                                        tCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCGGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTCCCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCgg                                                                                                  <  2:151033/250‑1 (MQ=255)
                                                                                                                                                  ttGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGgt                                                                                        <  2:562764/249‑1 (MQ=255)
                                                                                                                                                        ggTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCac                                                                                <  1:238634/251‑1 (MQ=255)
                                                                                                                                                         gTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCacg                                                                               <  1:102996/251‑1 (MQ=255)
                                                                                                                                                              cGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGAGATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGc                                                                         <  2:206046/251‑1 (MQ=255)
                                                                                                                                                               ggAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCg                                                                        <  2:536589/251‑1 (MQ=255)
                                                                                                                                                                gAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGa                                                                       <  1:482569/251‑1 (MQ=255)
                                                                                                                                                                   gATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTTGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGTGAATc                                                                    <  1:118880/251‑1 (MQ=255)
                                                                                                                                                                   gATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGGCATCACGACGGCGAATc                                                                    >  1:354690/1‑251 (MQ=255)
                                                                                                                                                                    aTCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCg                                                                   <  2:354690/251‑1 (MQ=255)
                                                                                                                                                                    aTCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGAGTCATCACGACggggatgcg                                                                    >  1:513569/1‑244 (MQ=255)
                                                                                                                                                                    aTCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCGTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCg                                                                   <  2:364738/251‑1 (MQ=255)
                                                                                                                                                                      ccATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGATCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACTGCGGTGTTCTCACCACGGAGAATCgg                                                                  >  2:394064/1‑249 (MQ=255)
                                                                                                                                                                           gagaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCAAGACGGGGAATAGTAATGa                                                             >  1:399220/1‑250 (MQ=255)
                                                                                                                                                                            agaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTcg                                                                                                                                     <  2:149408/179‑1 (MQ=255)
                                                                                                                                                                            agaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTcg                                                                                                                                     >  1:149408/1‑179 (MQ=255)
                                                                                                                                                                                   tggCTGGCTTAACTGGTCCTCCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGGAAGCACGTTCCCAGGTGGAGACGTGCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGAGAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTcc                                                     <  2:483388/248‑1 (MQ=255)
                                                                                                                                                                                     gCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCAt                                                   <  2:36437/250‑1 (MQ=255)
                                                                                                                                                                                        ggCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCATaa                                                 <  1:201740/249‑1 (MQ=255)
                                                                                                                                                                                          cTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGAACGGCGGTGGCGTCACGACGGCGAGACGTAAttaaagtttacaaaacagt                                             >  2:224338/1‑233 (MQ=255)
                                                                                                                                                                                             aaCTGGTCCTCCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCTGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCATAAAGTTCTg                                          <  1:206357/251‑1 (MQ=255)
                                                                                                                                                                                             aaCTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGGGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCATAAAGTTCTg                                          <  2:332442/251‑1 (MQ=255)
                                                                                                                                                                                              aCTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGAACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCAt                                                   <  1:310396/241‑1 (MQ=255)
                                                                                                                                                                                                   tCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGGGTCATCCCCACGGCGGATAGTAATGAAGATGTCCCTAAAATTCTGcccacc                                    >  2:133555/1‑247 (MQ=255)
                                                                                                                                                                                                        cTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGATGCATCACGACGGGCGATCGTAATGAAGTTGTCCCAAAAGTTCTGCCAGAtttttt                               >  1:453041/1‑251 (MQ=255)
                                                                                                                                                                                                           gtcgAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCATAAAGTTCTGCCAGCTTTTTTGCa                            <  1:216002/251‑1 (MQ=255)
                                                                                                                                                                                                              gAGGAACGCGTTGGTGGTGGCCTCAGCAGTATTGTAGGAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATGTGGCGCTAACAGCGGGACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCATAAAGTTCTGCCAGCTTTTTTGCAGAt                         <  1:546105/251‑1 (MQ=255)
                                                                                                                                                                                                               aGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGCATCGTAATGAAGGTGTCCATAAAGTTCTGCCAGCTTTTTTGCAGAtt                        >  2:292825/1‑251 (MQ=255)
                                                                                                                                                                                                                         tggtggtggCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCt                                                                                                                     <  2:420142/150‑1 (MQ=255)
                                                                                                                                                                                                                         tggtggtggCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCt                                                                                                                     >  1:420142/1‑150 (MQ=255)
                                                                                                                                                                                                                          ggtggtggCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCATAAAGTTCTGCCAGCTTTTTTGCAGATTGATACGCCAtt             <  2:513569/251‑1 (MQ=255)
                                                                                                                                                                                                                          ggtggtggCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGAGATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCATAAAGTTCTGCCAGCTTTTTTGCAGATTGATACGCCAt              <  2:436945/250‑1 (MQ=255)
                                                                                                                                                                                                                              gtggCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCATAAAGTTCTGCCAGCTTTTTTGCAGATTGATACGCCATTCGCt         <  1:447253/251‑1 (MQ=255)
                                                                                                                                                                                                                                     catcaGTATTGTAGTAAGCACGTTCCCAGGTGGATACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGTCCAGCAGGCGAGAGCGAATATTTTCGCGTAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCATAAAGTTCTGCCAGCTTTTTTGCAGATTGATACGCCATTCGCTGATGGa   <  2:412613/250‑1 (MQ=255)
                                                                                                                                                                                                                                     catcaGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCATAAAGTTCTGCCAGCTTTTTTGCAGATTGATACGCCATTCGCTGATGGaa  <  1:443478/251‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    
TAGCATGGTTATTCTCCTTGCGGCGCGGCTGCCGGAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCATAAAGTTCTGCCAGCTTTTTTGCAGATTGATACGCCATTCGCTGATGGAA  >  NC_000913/3987877‑3988356

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