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,0590CT94.3% 86.4 / NA 37P332S (CCG→TCG) dhaRDNA‑binding transcriptional dual regulator DhaR
Reads supporting (aligned to +/- strand):  ref base C (0/2);  major base T (21/12);  minor base G (0/2);  total (21/16)
Fisher's exact test for biased strand distribution p-value = 1.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGT‑TGCTCCATCCTGTGGAACAGATGC‑GGCAGTTGATGACCA‑GTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGAT‑CCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGG  >  NC_000913/1251820‑1252302
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
taccgtcccccGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGT‑TGCTCCATCCTGTGGAACAGATGC‑GGCAGTTGATGACCA‑GTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGAT‑TCGCAAACccg                                                                                                                                                                                                                                           <  1:480653/242‑1 (MQ=255)
     ttACCCGACGTATTGCATCAAGCAGTAAAACAGGCACATCAGCTCCAACACGTAGAAGCAACCTTTGCAAGCTAGCACCCGTTTATTGATGCGGTGATAACACTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGT‑TGCTCCATACTGTGGAACAGATGC‑GGCAGTTGATGAACA‑GTCACTTAGGAAAAGGCAGCCATACCTTCGCTCATATGCCACAGGACGAT‑TCGCaaaacacccgatt                                                                                                                                                                                                                                     >  2:318898/1‑241 (MQ=255)
                                       cacaTCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGT‑TGCTCCATCCTGTGGAACAGATGC‑GGCAGTTGATGACCA‑GTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGAT‑TCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGTCGGCGCGCAGTAGc                                                                                                                                                                                                   >  2:19394/1‑251 (MQ=255)
                                                tCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGT‑TGCTCCATCCTGTGGAACAGATGC‑GGCAGTTGATGACCA‑GTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGAT‑TCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTtcctgtc                                                                                                                                                                                          >  1:501811/1‑251 (MQ=255)
                                                  aaaCACGTAGAAGAAACCTTTAAAAGCCAGCACCATTTTATTGAGGCGGGGATACCCCTTAAACCGATAATAGAAACGCAGGGAACCAGTTTTATTTTGT‑GGCTCCATCCTGTGGAACAGATGC‑GGCGGTTAATGACCA‑GTCAATTAGGAAAAGGCAGCCATACCTTCGCTCATAGGCCACAGGCGGAT‑CCGCAACCCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTtcctgtcct                                                                                                                                                                                        <  2:95980/251‑1 (MQ=255)
                                                        gTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATCACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGT‑TGCTCCATCCTGTGGAACAGATGC‑GGCAGTTGATGACCA‑GTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGAT‑TCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTtg                                                                                                                                                                                  >  1:62175/1‑251 (MQ=255)
                                                             ggCAACTTTGGAAAGCCAGCTCCGGTTTATTGAGTCGGTGATAACCATTAAACCAAAAAAAAAAACGCAGGGACCCAGCTTGATTTTTT‑GGCTCCATCTTGGGGAACAG‑TGCCGGCATTTGATTCCCA‑GGCAATTGGAAAAAGTCACCCACCCGTTCGCTCATATCCCCCAGGTCTATTCCG‑AATCCCGCCGCTCGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGgag                                                                                                                                                                             <  1:102691/250‑1 (MQ=255)
                                                                            ccAGCCCCATTTTATGGATGCGGTGAGACCCCTTAAACCAATAATAAAAACGCGGGGAACAATCTTTTTTGTGTGTCCTC‑ATCCTGGGGAGCATATGC‑GGAAGTTG‑TGGCCAAGTTAATAAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGTGGTT‑TCGCAACCCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCCGTCCTGCTTTGTGGAGAAGAGGTCGTGGGCa                                                                                                                                                              <  2:378344/251‑1 (MQ=255)
                                                                              aGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGT‑TGCTCCATCCTGTGGAACAGATGC‑GGCAGTTGATGACCA‑GTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGAT‑TCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCa                                                                                                                                                              >  1:179985/1‑249 (MQ=255)
                                                                                    agattATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGCAACCAGCTTTATTTTGT‑TGCTCCATACTGTGGAACAGCTGC‑GGCCGTTCATGCCCC‑GTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGAT‑TCGCACACCCGCCGCTTGATTCATTTTGGTCGCCCGGCGGAGAGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTg                                                                                                                                                      >  2:599849/4‑251 (MQ=255)
                                                                                                   tGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGT‑TGCTCCATCCTGTGGAACAGATGC‑GGCAGTTGATGACCA‑GTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCCCAGGACGAT‑TCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTaa                                                                                                                                                  <  1:378760/240‑1 (MQ=255)
                                                                                                   tGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGT‑TGCTCCATCCTGTGGAACAGATGC‑GGCAGTTGATGACCA‑GTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGAT‑TCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTaa                                                                                                                                                  >  2:378760/1‑240 (MQ=255)
                                                                                                    gATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGT‑TGCTCCATCCTGTGGAACAGATGC‑GGCAGTTGATGACCA‑GTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGAT‑TCGCAAACCCGCCGCTTGATTCCTTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGCGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAAttt                                                                                                                                      >  2:435071/1‑250 (MQ=255)
                                                                                                    gATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGT‑TGCTCCATCCTGTGGAACAGATGC‑GGCAGTTGATGACCA‑GTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGAT‑TCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTc                                                                                                                                      >  2:161085/1‑251 (MQ=255)
                                                                                                                 ccGATAATAGAAACGCAGGGAACCAGCTTTATTTTGT‑TGCTCCATCCTGTGGAACAGATGC‑GGCAGTTGATGACCA‑GTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGAT‑TCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAagcg                                                                                                                          >  1:474656/1‑250 (MQ=255)
                                                                                                                  cGATAATAGAAACGCAGGGAACCAGCTTTATTTGGT‑TCCTCCAACCTGTGGACCAGATGC‑GGCAGTTGATGACCA‑GTAAATTAGAAAAAGTCAGCCATACCTCCCCTCAAATGCCACAGGACGAT‑TCGCAACCCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAagcgag                                                                                                                        <  2:616789/251‑1 (MQ=255)
                                                                                                                    ataataGAAACGCAGGGAACCAGCTTTATTTTGT‑TGCTCCATCCTGTGGAACAGATGC‑GGCAGTTGATGACCA‑GTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGAT‑TCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAagcgagcg                                                                                                                      >  1:21290/1‑251 (MQ=255)
                                                                                                                    ataataGAAACGCAGGGAACCAGCTTTATTTTGT‑TGCTCCATCCTGTGGAACAGATGC‑GGCAGTTGATGACCA‑GTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGAT‑TCGCAAACCCGCCGCTTGACTCATTTTGGTCGCCAGGCGGCGCACAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCCCTGCTAAGTCAGGCAATTCATAATGAAagcgagcg                                                                                                                      >  1:74455/1‑251 (MQ=255)
                                                                                                                            aaCGAAGGAAACAAGCTTTATTTTGT‑GGCTCCATCCTGTGAAACAGAGGC‑GGAAGTTAATGACCA‑GTAAATTAGGAAAAGTCACCCATCCCTTCGCTCATAGGCCACAGGACGAT‑TCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAg                                                                                                              <  2:63620/251‑1 (MQ=255)
                                                                                                                              cGCAGGGAACCAGCTTTATTTTGT‑TGCTCCATCCTGGGGAACAGATGC‑GGCAGTTGATGACCA‑GTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGAT‑TCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGt                                                                                                            <  1:538295/251‑1 (MQ=255)
                                                                                                                                  gggAACCAGCTTTATTTTGT‑TGCTCCATCCTGTGGAACAGATGC‑GGCAGTTGATGACCA‑GTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGAT‑TCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCtt                                                                                                        >  1:421067/1‑251 (MQ=255)
                                                                                                                                                   tgt‑tgCTCCATCCTGTGGAACAGATGC‑GGCAGTTGATGACCA‑GTCAATTAGGAAAAGTCAGCCATCCCTTCGCTCATATGCCACAGGACGAT‑TCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATtg                                                                                        <  1:581412/250‑1 (MQ=255)
                                                                                                                                                     t‑tgCTTCATACTGTGGACCAGATGC‑GGAAGTTGATGACCA‑GTCCATTAGGCAACGTCAGCCATACCTTCGCTCATATGACACCGGACGAT‑TCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCACTTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATtgtg                                                                                      >  2:416405/1‑250 (MQ=255)
                                                                                                                                                          tCCATCCGGGGGAACGGATGC‑GGCATTTGAGGACGA‑TTAAATTAGGAAAAGTCAGCCATACCTGCCCTCATATTCCGCAGGACGAT‑GGGTAACCCCGCCGCGTGATTCATTGTGGTCGCCAGGCGGCGCGCAGTAGCGTTCCGGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGt                                                                                   <  1:417907/249‑1 (MQ=255)
                                                                                                                                                                       aaaGATGC‑GGCGTTTGATGGCCA‑TTCAATTAGAAAAAGGCAGCCATCCCTTCGCTCATATGCCGCGGGACGAT‑GCGCAACCCCGCCGCTTGGTTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGGCCTTATATCGCCGTCAATTGTGAGTTATATGGTGAtgctg                                                                    <  2:358338/249‑1 (MQ=255)
                                                                                                                                                                                       gatgaCCA‑GTCAATTAGGAAAAGTCAGACATACCTTCGCTCATATGCCACAGGCCGAT‑TCGAAACCCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGgaagaa                                                      <  1:506810/250‑1 (MQ=255)
                                                                                                                                                                                              caggcAATTAGTAAAAGTCAGCCATCCCTTCGCGCATATGCAACAGGACGGT‑TCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGTAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATtg                                               <  1:387869/247‑1 (MQ=255)
                                                                                                                                                                                                      tAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGAT‑TCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCACGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGGGCTGCAGGTCCTTATATCTCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGACGGAAGAATTTATTGGTGGCggt                                       >  2:532437/1‑249 (MQ=255)
                                                                                                                                                                                                        ggAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGAT‑TCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTACGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCg                                     >  2:227480/1‑251 (MQ=255)
                                                                                                                                                                                                               tCAGCCATACCTTCGCTCATATGCCACAGGACGAT‑TCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGCGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGAc                              >  1:391069/1‑251 (MQ=255)
                                                                                                                                                                                                               ccAGCCATCCCTGCGCTCATAGGCCACAGGACGAT‑TCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGAc                              <  1:502715/250‑1 (MQ=255)
                                                                                                                                                                                                                 aGCCATACCTTCGCTCATATGCCACAGGACGAT‑TCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACggtcc                             >  2:202286/1‑247 (MQ=255)
                                                                                                                                                                                                                    cATACCTTCGCTCATATGCCACAGGACGAT‑TCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATg                          >  1:236864/1‑250 (MQ=255)
                                                                                                                                                                                                                                   atGCCACAGGACGAT‑TCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTg            >  2:189262/1‑249 (MQ=255)
                                                                                                                                                                                                                                   atGCCACAGGACGAT‑TCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTg            >  1:425356/1‑249 (MQ=255)
                                                                                                                                                                                                                                          aGGACGAT‑TCGCAATCCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATTCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTg   <  1:401351/251‑1 (MQ=255)
                                                                                                                                                                                                                                           ggACGAT‑TCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCTCGCAGTAGTTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATCCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTgg  <  1:599849/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
TAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGT‑TGCTCCATCCTGTGGAACAGATGC‑GGCAGTTGATGACCA‑GTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGAT‑CCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGG  >  NC_000913/1251820‑1252302

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