Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,510,609 C→T P327L (CCG→CTG)  yfeO → putative transport protein YfeO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,510,6090CT80.6% 78.6 / ‑5.4 40P327L (CCG→CTG) yfeOputative transport protein YfeO
Reads supporting (aligned to +/- strand):  ref base C (4/0);  major base T (22/7);  minor base G (0/7);  total (26/14)
Fisher's exact test for biased strand distribution p-value = 4.11e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCCGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAA  >  NC_000913/2510371‑2510831
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                              
cACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGt                                                                                                                                                                                                                    >  1:467489/1‑251 (MQ=255)
  cgcgATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCTAATCTTCCTGGCAGTGTTTGTCg                                                                                                                                                                                                                  <  1:201765/251‑1 (MQ=255)
           gCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCAt                                                                                                                                                                                                         <  1:14100/251‑1 (MQ=255)
           gCATCAAATGAAAAATCCGGGGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGTACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCAt                                                                                                                                                                                                         <  2:509957/251‑1 (MQ=255)
           gCATCAAATGAAAAATCCGGGGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGTGTGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTGTCAGCACGAGCGATTCCGTTTTGCTGGCGGTAAGTAAACTTGCGGCGCTGGTCGTAGCTGCCGCCAGGGGCTTTCGCGGTGGGCGATTCTTCCTGGCAGCGTTAGTCGGCGTGGCa                                                                                                                                                                                                          <  2:682585/250‑1 (MQ=255)
                     aaaaaTCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAAATTGCCGCCCCGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGGGAATCTTTCTGGCCGTGTTTGTCGGCGGGGCATTAggggtggt                                                                                                                                                                                               >  1:316420/1‑246 (MQ=255)
                           ccGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGAtgctgta                                                                                                                                                                                         <  1:38198/251‑3 (MQ=255)
                               tGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTTCCGCCCTGGTCGTTGCTGCCGCCCGTGGCTTTCGCGGTGGGCGGATCTTCCTGGCAGTGTTTGTCGGCGGGGGCTTAGGGGTGATGCTGCATg                                                                                                                                                                                       >  1:260283/1‑249 (MQ=255)
                                          ggTATTGGCGGATTTATTCGCGGTATTCGGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCccgc                                                                                                                                                                          <  2:593745/251‑1 (MQ=255)
                                             aTTGGCGGTTTTTTTCTGGGTATGGTGGGGGTTATTGTTGGACAGGTTGGGTTTTTAAAAGGGGTGGATTAGGTTCAGAAGATGGTGGAAATTCAGGCTTTCAGCACCAGCGGTTACTTTTTGCTGGCGGTAATTAAACTTGCGGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAACTTTCCTGGCGGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCTCGTGCccgccg                                                                                                                                                                        <  2:171736/250‑1 (MQ=255)
                                                  cGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGCATCTTCCTGGCAGTGTTTGTCCGGGGGGCCTTTGGGTTGTTGCTGCATGAGCACGTTCCCGCCGTAcc                                                                                                                                                                    >  2:97779/1‑249 (MQ=255)
                                                                 tATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGCTGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGTCAGTGGCTTTTGCGGTGGCGGAATCTTCCCGGCAGTGGTTGGCGGCGGGGGATTAGGGTTGATGCTGCATGAGCACGTACCCGCCGTGACAGCGGACATAACCggtt                                                                                                                                                   >  2:591299/1‑248 (MQ=255)
                                                                       gggggTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGGGGGCGAATCTTCCTGGCAGGGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGGACCAGCGGCAATAACCGTTTCTTgcg                                                                                                                                             >  1:496514/1‑251 (MQ=255)
                                                                                  gtggACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGGGGCCTTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGc                                                                                                                                  >  1:288409/1‑251 (MQ=255)
                                                                                              tCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCAGGCGGTCATTAAACTTGACGCCCTGGTCGCTGCTGCCGCCAGTGGATTTCGCGGTGGGCGAATCCTCTTGGCCGGGTTTGTCGGCGTGGCATTAGGGTTGATGCTGAATGAGCACGGTCCCGCCGTTCCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATAGTGCtggt                                                                                                                       >  1:685565/1‑251 (MQ=255)
                                                                                                 cTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGATCATTCCACTTGCCCCCCTGGTCCTTTCTCCCGCCCGTGGCCTTCCCGGTGGGCGCATCTTCCCGGCAGTGTTTGTCGGCGTGGCATTAGGGGTGTTGATGCCTGAGCACGTTCCCGCCGTAACAGCGGCAATAACCGTTTCTTGCCCTATTTTCGGCATCGTGCTggtggt                                                                                                                    >  1:688145/1‑251 (MQ=255)
                                                                                                      tAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCCGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGGGGCATTTGGGTTTGTGCTGCATGGGCACGTTCCCGCCGGTCCCGCGGCAATAACCGTTTCTTGCGCTTTTTTCGGCATCggtctggggggaacac                                                                                                               >  2:404812/1‑236 (MQ=255)
                                                                                                         aggggCTGGAGAAGATGCAGCAGATGGTGGAAAATCAGTCTTTCACCCCCAGTTTTTTCTTTTGGGTGCCGTTAATAAACTTTGCCGCCCGGGTGGTTGCTGCCGCTAGTGGTTTTCGGGGGGGGTTATTTTTCCGGGTAGTGTTGGTCGGCGTGGCATTGGGGTTGATGTTGATTGATCACGTCTCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACAcgc                                                                                                             <  1:737440/249‑1 (MQ=255)
                                                                                                                ggagaagaggcagaaGAGGGTGGCAAATCAGTTTTTCACCACCAGCGTTTTTTTTTTGGTGGCGTTAATTAACTTCGCCCCCCTGGTTGTTGCTGCCGCCAGGGGTTTGCGGGGTGGGCGAATCTTCCGGGCGGTGTTTGGCGGCGTGGCATTGGGGTTGATGCTGCATGAGCCCTTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCt                                                                                                      <  1:529737/242‑1 (MQ=255)
                                                                                                                ggATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGcg                                                                                                                                 >  1:247479/1‑223 (MQ=255)
                                                                                                                        tGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTAACAGCGGCAATAACCGTTTTTTGCGCTATTCTCGGCATCGTGCTGGGGGTAACACGCGGTGGCTGGTTAAgg                                                                                              >  2:13411/1‑250 (MQ=255)
                                                                                                                          cagcagATGGTGGCAAATCCGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGAAATTAAAATAGCCCCCCTGGTCGTTGCAGCCGCCAGGGGCTTTCGCGGTGGGGGAATCTTTCTGGCAGGGTTTGTCGGCGTTGCATTAGGGTTGATGATGCATGACCACATTCCCGCCATACCACGGGCCATAACCGTTTCTTTCGCCTTTCCCGACATCGCGCTGGTGGTAACACGCGATGGCTAGTTAAGTCt                                                                                            >  2:443649/1‑251 (MQ=255)
                                                                                                                                      gCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGAAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGCTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGGGCTGGGGGTAACACGCGATGGCTGGGTAAGTCTTTTTTTggcggg                                                                                >  1:268985/1‑250 (MQ=255)
                                                                                                                                         aaTCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCCGTGGCTTTCGCGGTGGGCGAATCTTCCCGGCAGCGTTTGTCGGCGTGGCATTAGGGGTGATGCTGCATGAGCACGTTCCCGGCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGGGGTAACACGCGATGGCTGGGTAAGTCTTTTTATTGCGGCa                                                                               >  1:369297/1‑249 (MQ=255)
                                                                                                                                            cagtctttcagccccatggtttcttttttGCGGGCGTTATTTAACCTTGCCGCCGTGGTTGTGGTTCCCGCTGGTGGCTTTCGGGGTGGGCGATTCTTCCGGGCTGTGTTTGGCGGCGGGGCATTAGGTTGGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATTGCTGGTTAAGTCTTTTTATGGCGGCAGTCGt                                                                          <  1:76011/227‑1 (MQ=255)
                                                                                                                                                aatctagccccAGCGATTCTTTTTTGCGGGGGGAAATAAAACTTGCCGCCCTGGTCGTTGCTCCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCGGGCAGTGTTGGTCGGCGTGTCATTAGGGTTAATGCTCCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTGTa                                                                      <  1:378605/242‑1 (MQ=255)
                                                                                                                                                     aGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGATTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCCTTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTTGGCTATTCTCGGCATCGTGCTGGGGGTAAAAAGCGGTGGCTGGTTTAGTCTTTTTTTggggggagggggtgttaagtt                                                                 >  2:378605/1‑232 (MQ=255)
                                                                                                                                                              gATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCCGTGTTTGTCGGCGTGGCATTAGGGGTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGTTGGGAACACGCGATGGCTGGGTAAATCTTTTTATTGCGGGCGTCGttttttcccaaaaccacct                                                        >  2:25036/1‑234 (MQ=255)
                                                                                                                                                              gATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGGCGTTGCTGCCGCCCGTGGCTTTCGCGGTGGGCGCATCTTCCTCGGCGTGTTTGTCGGCGTGGGCTTAGGGGTGATGCTGCATGGGCACGTTCCCGCCGTAACAGCGGCAATAACCGTTTATTGCGCTATTCTCCGCCTCGTGCTGGGGGTAAAACGCGATGGCTGGGTAAATCTTTTTTTGGCGGCCGTCGTTGTTCCCCATACcaccat                                                        >  1:368559/1‑248 (MQ=255)
                                                                                                                                                               atttttttttgtgggtggTAATAAACCTTGCCGCCCGGGTCGTTGCTGCCGCCGGTGGCTTTCGGGGGGGGCAATTCTTCCGGGCGGTGTTTGGCGGCGTGGCATTGGGGTTGATGCTGCGTGAGCACTTCCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTGTACCCAATACCACATTg                                                       <  1:146240/235‑1 (MQ=255)
                                                                                                                                                                ttgtttttgggggggggTAATAAACCTTGCCGCCGTTGTTGTGGCTGCCGCTAGGGGTTTTCGGGGTGGGGAAATCTTCCGGGCAGTGTTGGGCGGCGGGGCATTAGGGTTGATGCTGCATGAGCGCGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTGTACCCAATACCACATTg                                                       <  1:393837/235‑1 (MQ=255)
                                                                                                                                                                        tGCTGGCGGTAACTAAACTTGCCCCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGCTGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGGCGTTGTACCCCATAccccaatgttacccccg                                              >  2:44797/1‑236 (MQ=255)
                                                                                                                                                                               ggTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCATGCGTTGGACCCAATACCCCATTTCTACCGCTGCCCTGt                                        >  1:411007/1‑250 (MQ=255)
                                                                                                                                                                                  aaTTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTGTACCCAAAACCACATTTCTACCGCTGctcctgatct                                    >  2:486238/1‑244 (MQ=255)
                                                                                                                                                                                   aTTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCg                                                                                                                                    >  2:250049/1‑154 (MQ=255)
                                                                                                                                                                                     tAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGATGATGCTGAATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTATCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGACGTTTTAACCCATACCACATTGCTTCCGCTGCTCTGTATAGTc                                  >  1:217335/1‑250 (MQ=255)
                                                                                                                                                                                        aCTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTGGACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCCTGc                              >  1:210219/1‑251 (MQ=255)
                                                                                                                                                                                                   cTGGTCGTTGCTGCCGCCCGTGGCTTTCGCGGTGGGCGACTCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCCGTCGTTGTACCCAATTCCACCTTGCTTCCGCTCCTCTGTagtgtcttggttcccgcctg                    >  2:287573/1‑232 (MQ=255)
                                                                                                                                                                                                          ttGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGGCGTTGTACCCAATACCCCATTGTTACCCCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGttt             >  2:601102/1‑250 (MQ=255)
                                                                                                                                                                                                                     ggtggcttggggggggggCGAATCTTCCGGGCGGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTaa  <  1:25036/236‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                              
CACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCCGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAA  >  NC_000913/2510371‑2510831

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