Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,252,059 C→T P332S (CCG→TCG)  dhaR → DNA‑binding transcriptional dual regulator DhaR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,252,0590CT92.3% 91.2 / ‑3.8 39P332S (CCG→TCG) dhaRDNA‑binding transcriptional dual regulator DhaR
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (21/15);  total (22/17)
Fisher's exact test for biased strand distribution p-value = 5.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACC‑GATAATA‑GAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAG‑CCATACCTTCGCTCATATGCCACAGGACG‑ATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCT  >  NC_000913/1251811‑1252300
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                 
acaaactcttcgttttccccGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACC‑GATAATA‑GAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAG‑CCATACCTTCGCTCATATGCCACAGGACG‑ATc                                                                                                                                                                                                                                                   <  1:245152/232‑1 (MQ=255)
                 cccGCCGTAATAAAACAGGCGATAAAACCGGAACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAACCC‑GATAATA‑GAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGCACAAATGCGGAAGTTGATGACCCGTCAATTAGGAAAATTAAG‑CCATACCTTCGCTCATATGCCAAAGGACG‑ATTCGCAAAACCGCCACttact                                                                                                                                                                                                                                >  2:672765/1‑248 (MQ=255)
                   cgccgTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACC‑GATAATA‑GAAACGCAGGGTACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAG‑CCATACCTTCGCTCATATGCCACAGGACG‑ATCCGCAAACCCGCCGCTTGATTc                                                                                                                                                                                                                              >  1:274385/1‑251 (MQ=255)
                            gCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACC‑GATAATA‑GAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAG‑CCATACCTTCGCTCATATGCAACAGGACG‑ATTCGCAAACCCGCCGCTTGATTCATTTTGGTc                                                                                                                                                                                                                     >  1:598177/1‑251 (MQ=255)
                               ccaaGCCATAAAACAGGCACATCCGCTCAAACACGTAGAAGCACCCTTTGACAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTCAACC‑GATAATA‑GAAACGCAGGGAACCAGCTTTCTTTTGTTGATCCATCCTGTGGAACAGCTGCGGCAGATGATGAACAGTCATTTAGGAAAAGTAGG‑CCAAACCTTCGCTCATATGCCACAGAACG‑ATTGGCAAACCCGCCGCTTGATTCCTTTTgtgcgcc                                                                                                                                                                                                                  >  2:16158/2‑245 (MQ=255)
                                      aTAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACC‑GATAATA‑GAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAG‑CCATACCTTCGCTCATATGCCACAGGACG‑ATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAggcggg                                                                                                                                                                                                           >  1:597955/1‑250 (MQ=255)
                                                 ccacccGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACC‑GATAATA‑GAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAG‑CCATACCTTCGCTCATATGCCACAGGACG‑ATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCt                                                                                                                                                                                                <  1:12406/247‑1 (MQ=255)
                                                     ccGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACC‑GATAATA‑GAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAG‑CCATACCTTCGCTCATATGCCACAGGACG‑ATTCGCAAACCCGCCGCTTGATTCATTTTGGGCGCCAGGCGGCGCGCAGTAGCTTtcc                                                                                                                                                                                            >  2:310596/1‑251 (MQ=255)
                                                            aaCACGTAGAAGCACCCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACC‑GATAATA‑GAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAG‑CCATACCTTCGCTCATATGCCACAGGACG‑ATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTtcctgtcctg                                                                                                                                                                                     <  1:310596/251‑1 (MQ=255)
                                                            aaCACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAACCC‑GATAATA‑GAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAG‑CCATCCCTTCGCTCATATGCCACAGGACG‑ATTCGCAACCCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTtcctgtcctg                                                                                                                                                                                     <  1:660349/251‑1 (MQ=255)
                                                              cacGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACC‑GATAATA‑GAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAG‑CCATACCTTCGCTCATATGCCACAGGACG‑ATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCt                                                                                                                                                                                   >  1:182470/1‑251 (MQ=255)
                                                                  tAGAAGCACCCTTTGAACGCCAGCACCAGTTTATTGATGGGTTAATAACCCTTAA‑CC‑GGAAATAAGAACCGCAGGGAACCATCTTTTTTTTGTTGCTCCACCCTGTGGAACAGAGGCGGCAGTTGATGACCATTAAATTAGGAAAATTCAG‑CCATCCCTCCGCTCATATGCCACAGGGCG‑TTTCGAAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTtg                                                                                                                                                                                <  1:55843/250‑1 (MQ=255)
                                                                   aGAAGCAACTTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACC‑GATAATA‑GAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACAAGTCAATTAGGAAAAGTCAG‑CCATACCTTCGCTCATATGCCACAGGACG‑ATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTtgt                                                                                                                                                                               <  1:668862/250‑1 (MQ=255)
                                                                            cTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACC‑GATAATA‑GAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAG‑CCATACCTTCGCTCATATGCCACAGGACG‑CTTCGCAAACACGCCGCTTGATTCATTTTGGTCGCCcggcgg                                                                                                                                                                                                            >  2:540362/1‑212 (MQ=255)
                                                                               tGAAAGCAAGCCCCAGTTTATTGTGGGGGGGAAAACCCTAAAACC‑GAAAATA‑GAAACGCAGGGAACCAGCTTAATTTTTTTGCTCCACCCGGTGGAACAAGGGCGGCATTTGATGACCAGTCAAGTAGGAAAAGACAG‑CCTTACTTCCCCTCATATGCCAAAGGACG‑ATTCGCAAACCCGCCGCTTGATTCATTTTTGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGc                                                                                                                                                                   <  1:471904/250‑1 (MQ=255)
                                                                                        cacccgggttatggtgtcGGTGATA‑CCCTTAAACCCGATAAAA‑GAAACGGAGGGAACAAGCTTATTTTTGTAGCTCCTTCTTGTGCAACAGTTGCGGCATGTGACGCCCAGAAAGTTAGAAAAAGGCAG‑CTATCCCTGCGCTCAGATGCCCCGGGCTG‑ATTCGAAAACCCGCCGCGTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAgg                                                                                                                                                         <  2:438343/234‑1 (MQ=255)
                                                                                                             gATAACCCTTAAACC‑GATAATA‑GAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAG‑CCATACCTTCGCTCATATGCCACAGGACG‑ATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTc                                                                                                                                    >  1:573996/1‑251 (MQ=255)
                                                                                                                accccTTAAACC‑AAAAAAA‑AAAACGTAGGGAACTTGTTTTTTTTTGTTGCTCCATCCTGGTGAACAGATGCGGCAGTTGTGGACAATACAATAAGAAAACGTGGG‑CCATCCCCTCCCATATATGCCAC‑GGGCGTATTCGCAAAGCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCGGTTCTGTCCTGCTTTGTGGAGAAGAGGGCGTGTGCAAGGCACTGCTAAGTCAGGCAATTCATa                                                                                                                                 <  1:61074/249‑1 (MQ=255)
                                                                                                                    cTTAAACC‑GATAATA‑GAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAG‑CCATACCTTCGCTCATATGCCACAGGACG‑ATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATg                                                                                                                              >  1:167928/1‑250 (MQ=255)
                                                                                                                               aaaaaagaaccggaGGGAACCAGCTTTATTTTTTTGCTCCATCCTGGGAAAAAGGGGCGACAGTTGATGACAAGTCAATAAGGAAAAGTCAG‑CTTTACCTTCCCTCATATCCCCCAGGTCG‑ATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAagcgag                                                                                                                      <  1:510497/236‑1 (MQ=255)
                                                                                                                                   a‑GAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAG‑CCATACCTTCGCTCATATGCCACAGGACG‑ATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGtgctg                                                                                                               >  1:304804/1‑251 (MQ=255)
                                                                                                                                     gAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAG‑CCATACCTTCGCTCATATGCCACAGGACG‑ATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAAGGAAAGCGGACGtgctgc                                                                                                              >  2:202825/1‑251 (MQ=255)
                                                                                                                                           gggggAACCAACTTTTGTTTGTTGCTCCACCGTGTGGAACAGAGGCGGAAGTTGAGGACCTGTCAATTAGAAAAAGTCCG‑CCATCCCTCCCCTTATATGCCCCAGGACG‑TTTCGAAACCCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTcc                                                                                                        <  2:182470/249‑1 (MQ=255)
                                                                                                                                            aggaaaCCAGCTTTATTTTGTTGCTCCACCCTGTGGAACAAATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAG‑CCATCCCTTCGCTCATATGCCACAGGACG‑ATTCGCAACCCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCt                                                                                                       <  1:418529/247‑1 (MQ=255)
                                                                                                                                             gggAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAG‑CCATACCTTCGCTCATATGCCACAGGACG‑ATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAgg                                                                                                                                                         >  2:351544/1‑200 (MQ=255)
                                                                                                                                             gggAACCAGCTTTATTTTGTTGCCCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAG‑CCATCCCTTCGCTCATATGCCACAGGACG‑ATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAgg                                                                                                                                                         <  1:351544/200‑1 (MQ=255)
                                                                                                                                                    aGCTTTAGTTTGTTGATACATCGTGTGGAATAGATGCGGCAGCTGATGACCAGTCATATAGGAAAAGTCTGACCA‑ACCTTCGATCATATGCCACAGAACG‑ATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGAGCGCAGTAGCTTTCCTGTCATGCTTTGTGACGAAGAGGACGTGAGCAAGGCACTGCTAAGTCAGGCAATTCATAATGCAAGCGAGCGGTCTACAGGTCATTAGATCGc                                                                                               >  2:122267/1‑251 (MQ=255)
                                                                                                                                                                         ccTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAG‑CCATACCTTCGCTCATATGCCACAGGACG‑ATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTTAATTGTGAGTTATATgg                                                                          >  1:391832/1‑251 (MQ=255)
                                                                                                                                                                           tgtgGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAG‑CCATACCTTCGCTCATATGCCACAGGACG‑ATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTTTGAGTTATATgg                                                                          >  2:537554/1‑249 (MQ=255)
                                                                                                                                                                           ggtgGAACAGATGCGGCAGTTGATGCCCAGTCAATTAGGAAAAGTCAG‑CCATCCCTTCGCTCATATGCCACAGGCCG‑TTTCGCAACCCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTg                                                                        <  1:487374/250‑1 (MQ=255)
                                                                                                                                                                            gtgGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAG‑CCATACCTTCGCTCATATGCCACAGGACG‑ATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGGCCTTATATCGCCGTCACTTGTGAGTTATATGGTgg                                                                       >  1:618918/1‑250 (MQ=255)
                                                                                                                                                                                     atggggaagtgtaggaCCAGTAAATTAGGAAAAGTAAG‑CCATCCCTCCGCTCATATGCCCCAGGCCG‑TTTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCt                                                              <  1:480920/237‑1 (MQ=255)
                                                                                                                                                                                        cGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAG‑CCATACCTTCGCTCATATGCCACAGGACG‑ATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGc                                                           >  1:534404/1‑251 (MQ=255)
                                                                                                                                                                                          gCAGTTGATGACCAGTCAATTAGGAAAAGTCAG‑CCATACCTTCGCTCATATGCCACAGGACG‑ATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCgg                                                         >  1:443395/1‑251 (MQ=255)
                                                                                                                                                                                                    aCCAGTCAATTAGGAAAAGTCAG‑CCATACCTTCGCTCATATGCCACAGGACG‑ATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTAt                                               <  1:537554/251‑1 (MQ=255)
                                                                                                                                                                                                                ggAAAAGTCAA‑CCATACCTTCGCTCATATGCCACAGGACG‑ATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCg                                   >  1:50798/1‑251 (MQ=255)
                                                                                                                                                                                                                        cAG‑CCATACCTTCGCTCATATGCCACAGGACG‑ATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCACTTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGATGCGCTGGCGGAAGAATTTTTTGGTGGCGATCGCGCgggc                            >  2:66533/1‑248 (MQ=255)
                                                                                                                                                                                                                             catcccTTCGCTCATAGGCCACAGGACG‑AGTCGCAACCCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGa                       <  1:131734/247‑1 (MQ=255)
                                                                                                                                                                                                                                                 acaGGACG‑ATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGc   >  1:524738/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   gagtcgtttCCGAAAACCCGCCGCTTGAATCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTGCCTGGCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAAGTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTAGATCGCCGTTAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCt  <  2:562004/243‑1 (MQ=255)
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                 
CTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACC‑GATAATA‑GAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAG‑CCATACCTTCGCTCATATGCCACAGGACG‑ATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCT  >  NC_000913/1251811‑1252300

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