Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,490,095 G→A intergenic (+424/‑807) aldA → / → cybB aldehyde dehydrogenase A/superoxide oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,490,0950GA97.4% 107.0 / NA 38intergenic (+424/‑807)aldA/cybBaldehyde dehydrogenase A/superoxide oxidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (20/17);  minor base T (0/1);  total (20/18)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

GAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTGGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAGTGGTGCATGACGCCACGGAAACAATGGTGTCGTTGCCATCCAGAGTGTCGTCATTGACGTTATAAACGATAGTTTTCATTTCA  >  NC_000913/1489849‑1490342
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
gaaccgaaATGGCAGCCAATGATATCGGAAGTGACTATTTCTTCATCGGTATAACCAAATGACTCGTAATTGGTGGTTGCTTGTTTTAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGTGAATCGATACCAGTTCTGTGACCGAACCTGTTTTCACCGGCGCGCGTTGCGCTTGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGGCGCTCCCGTAGTGTaggg                                                                                                                                                                                                                                                      >  2:332784/1‑250 (MQ=255)
  accgaaATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTGTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTGTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAAt                                                                                                                                                                                                                                                   <  2:296192/251‑1 (MQ=255)
     gaaATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGat                                                                                                                                                                                                                                                >  1:455695/1‑251 (MQ=255)
                cAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCg                                                                                                                                                                                                                                     >  1:416506/1‑251 (MQ=255)
                    gATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTgc                                                                                                                                                                                                                                 >  1:55110/1‑251 (MQ=255)
                    gATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTgc                                                                                                                                                                                                                                 >  1:508544/1‑251 (MQ=255)
                             agagACTATTTCTTCATCGGTCTAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGacg                                                                                                                                                                                                                         >  2:603399/1‑248 (MQ=255)
                                          tcATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCtt                                                                                                                                                                                                             <  1:339643/249‑1 (MQ=255)
                                              cGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCAc                                                                                                                                                                                                       >  1:240326/1‑251 (MQ=255)
                                                    aaCCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCCCGTAAATCTTTACCACGCGGGc                                                                                                                                                                                                 >  1:295939/1‑251 (MQ=255)
                                                         aaTGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATAGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATcc                                                                                                                                                                                            >  1:352511/1‑251 (MQ=255)
                                                            gccTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCATTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCCGTTTGCCGCTCAGTTCCGGGATCACCAGACCACTGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGTTATTTTCTGCCGCTGCGCGTGACGCACGTAAATCTTTACCACGCGGTCCATccacc                                                                                                                                                                                         >  2:112869/3‑251 (MQ=255)
                                                            gACTCGTTATTGGTGGTTGCTTGTTTACGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATccacc                                                                                                                                                                                         <  1:452438/251‑1 (MQ=255)
                                                                cGTTATTGGTGGTTGCTTGCTTAAGTGCGTTATTCACCTCGTCGGCAGTCACGTTTTTCCGAGAAACGATACCAGTTCAGTGACCTAACCTGTTTTCACCGGCACGCGGTGCGCATGACCTTTCAGTTGGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAAGGATATTTACTGCCTCTGCGCGTGAATCCCGTAAATCTTTACCACGCGGGCCAGCCACCAGTg                                                                                                                                                                                     <  2:109301/251‑1 (MQ=255)
                                                                 gTTATTGGCGGTTGCTTGTTTACGTGCGTTCTTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGACCCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTg                                                                                                                                                                                     >  2:297216/1‑250 (MQ=255)
                                                                                 tGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGtata                                                                                                                                                                     <  2:558428/250‑1 (MQ=255)
                                                                                   tcagaaaatatcaTTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTTGGGAATGATATTTTCTGc                                                                                                                                                                                                                                       <  1:355494/170‑1 (MQ=255)
                                                                                     tAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCTCGCGGGCCAGCCACCAGTGACTGGGTGCCAGTATAGGCa                                                                                                                                                                 <  1:138965/250‑1 (MQ=255)
                                                                                       aGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGCCTGGGTGCCAGTATAGGCATg                                                                                                                                                               >  2:171875/1‑250 (MQ=255)
                                                                                                   aCCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATc                                                                                                                                                  >  2:246690/1‑251 (MQ=255)
                                                                                                     ctctTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCCTGACCTTTCAGTTTGACGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGCGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTACATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGt                                                                                                                                                >  1:288060/1‑251 (MQ=255)
                                                                                                                              gaATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGCCTTCTATCCCGAAACTGTc                                                                                                                        >  1:273065/1‑250 (MQ=255)
                                                                                                                              gaATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTGTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCa                                                                                                                       <  1:315990/251‑1 (MQ=255)
                                                                                                                               aaTCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTTCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCCATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGAGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCAt                                                                                                                      >  2:505097/1‑251 (MQ=255)
                                                                                                                                   gATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCa                                                                                                                   >  1:260173/1‑250 (MQ=255)
                                                                                                                                                                      cACGCGTTGCGCATGACCTTTAAGTTTGCCGCTCCGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGAATGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGTCTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAGTGGTGc                                                                               >  2:279985/1‑251 (MQ=255)
                                                                                                                                                                                 cATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAGACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAGTGGTGCATGACGCCACg                                                                    <  1:235711/251‑1 (MQ=255)
                                                                                                                                                                                          tCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAGTGGTGCATGACGCCACGGAAACAATg                                                           >  2:108644/1‑251 (MQ=255)
                                                                                                                                                                                                 gCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCTAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAGTGGGGCAGGACGCCACGGAAACAATGGTGTCGt                                                    <  2:295939/251‑1 (MQ=255)
                                                                                                                                                                                                            ccGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTGCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCGTTGGCCATCGGCGCAAGACAGTTAGTGGTGCATGACGCCACGGAAACAATGGTGGCGTTGCCATCCaga                                         <  1:422958/251‑1 (MQ=255)
                                                                                                                                                                                                                 aTCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAGTGGTGCATGACGCCACGGAAACAATGGTGTCGTTGCCATCCAGAGTgtc                                    <  2:541210/251‑1 (MQ=255)
                                                                                                                                                                                                                   cACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAGTGGTGCATGACGCCACGGAAACAATGGTGTCGTTGCCATCCAGAGTgtcg                                   <  2:290063/250‑1 (MQ=255)
                                                                                                                                                                                                                       aGACCAATGGCTTTTGCCGCCCCCGAAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAGTGGTGCATGACGCCACGGAAACAATGGTGGCGTTGCCATCCAGAGTGTCGTCAtt                              >  2:612436/1‑251 (MQ=255)
                                                                                                                                                                                                                        ggccAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCGCCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCAGCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCAGCGCCGCAAGACAGGTAGTGGTGCATGACGCCACGGAAACAATGGTGTCGTTGCCACCCACAGTGTCATCATTg                             <  2:82026/249‑1 (MQ=255)
                                                                                                                                                                                                                           cAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAGTGGTGCATGACGCCACGGAAACAATGGTGTCGTTGCCATCCAGAGTGTCGTCATTGACg                          <  1:11729/251‑1 (MQ=255)
                                                                                                                                                                                                                             agggCTTTTGCCGCCCCCGTAGTGTAGGGAATGATAGTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCTCAAGACAGTTAGTGGTGCAGGACGCCACGGAAACAATGGTGTCGTTGCCATCCAGAGTGTCGTCATTGACGtt                        <  2:273065/249‑1 (MQ=255)
                                                                                                                                                                                                                             aTGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAGTGGTGCATGACGCCACGGAAACAATGGTGTCGTTGCCATCCAGAGTGTCGTCATTGACGtt                        >  1:440997/1‑251 (MQ=255)
                                                                                                                                                                                                                               ggCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAGTGGTGCATGACGCCACGGAAACAATGGTGTCGTTGCCATCCAGAGTGTCGTCATTGACGTta                       <  2:455695/250‑1 (MQ=255)
                                                                                                                                                                                                                                    ttGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAGTGGTGCATGACGCCACGGAAACAATGGTGTCGTTGCCATCCAGAGTGTCGTCATTGACGTTATAAAc                  <  1:556452/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   ggtgGGGAATGATATTTTCCGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCAGCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGCCTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAGTGGTGCATGACGCCACGGAAACAATGGTGTCGTTGCCATCCAGAGTGTCGTCATTGACGTTATAAACGATAGTtttcatttca  <  2:439422/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
GAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTGGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAGTGGTGCATGACGCCACGGAAACAATGGTGTCGTTGCCATCCAGAGTGTCGTCATTGACGTTATAAACGATAGTTTTCATTTCA  >  NC_000913/1489849‑1490342

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