Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,003,502 G→T P226Q (CCG→CAG)  rarD ← putative transporter RarD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,003,5020GT93.3% 120.2 / ‑6.5 46P226Q (CCG→CAG) rarDputative transporter RarD
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base T (22/20);  minor base C (3/0);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 2.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCGGTACGGTAG‑TGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACAT‑GCCGGTTTGCGC‑TTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACC  >  NC_000913/4003258‑4003737
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              
cGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGctgtacgg                                                                                                                                                                                                                                          >  1:279487/1‑251 (MQ=255)
     tatTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACCC‑AACAGCTGTACGGTAG‑Tg                                                                                                                                                                                                                                    >  1:361528/1‑251 (MQ=255)
         ttataagaatcaatttagttTTGAATTGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGACACCCAAATAAAGGCGAAAGGCACAATAGTATCGGCACCCGGTTTTTAAACATAAAACGTAACAGCCAGCAGGAACATCAGCGTCGGGCCAATTTACTGGAAAAAGCCGAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACTC‑AACAGCTGTAAGGTAG‑TTACa                                                                                                                                                                                                                                 <  2:499887/238‑1 (MQ=255)
             cAGAATATATTTATTTTTTATTTGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGTAG‑TGACAATAcc                                                                                                                                                                                                                            >  1:120010/1‑251 (MQ=255)
                         attttttatttGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGTAG‑TGACAATACCGGCGGCGATca                                                                                                                                                                                                                 <  1:319268/250‑1 (MQ=255)
                          ttttttatttGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGAAAAAAAC‑AACAGCCGTACGGGAG‑TGACaaaaacgggggcggtacac                                                                                                                                                                                                               >  2:440133/1‑234 (MQ=255)
                           tttttatttGGACGTTCTACGCTGCGTATACATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAAAAC‑AACAGCTGTAAGGGAG‑TGACAATACCGGCGGCGATcagc                                                                                                                                                                                                               >  1:162755/1‑250 (MQ=255)
                            ttttatttGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGTAG‑TGACAATACCGGCGGCGATcagca                                                                                                                                                                                                              <  2:361528/250‑1 (MQ=255)
                             tttatttGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGTAG‑TGGAAATACCGGCGGCGacccgc                                                                                                                                                                                                               >  2:243454/1‑243 (MQ=255)
                              ttatttGGACGTTCTACGCTGCGTATAAATCGCATCCATCACACAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGTAG‑TGACAATACCGGCGGCGAACCGCAGt                                                                                                                                                                                                            >  1:55908/1‑250 (MQ=255)
                                   tGGACGTTCTACGCTGCGTATAAAGCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGTAG‑TGACAATACCGGCGGCGATCAGCAGTAAAtt                                                                                                                                                                                                       <  2:263310/250‑1 (MQ=255)
                                    ggACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGTAG‑TGACAATACCGGCGGCGATCAGCAGTAAATTCa                                                                                                                                                                                                     <  1:243454/251‑1 (MQ=255)
                                     gACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGTAG‑TGACAATACCGGCGGCGATCAGCAGTAAATTCAg                                                                                                                                                                                                    <  1:243626/251‑1 (MQ=255)
                                       cGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGTAG‑TGACAATACCGGCGGGGATCAGCAGTAAATTCAGCg                                                                                                                                                                                                  >  1:556410/1‑251 (MQ=255)
                                         ttCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCCGTACGGTAGCGGACA‑CACCGGCGGCGGTCAGCAGGAAATTCAGCGAc                                                                                                                                                                                                >  2:330617/1‑251 (MQ=255)
                                                                     cacaAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCGTATCGGCACCCGTTTGTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGGTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGTAG‑TGACAATACCGTCTGCGATCAGCAGGAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTg                                                                                                                                                                    <  2:500228/251‑1 (MQ=255)
                                                                         aaaaTGGCCAGCGCCACCCAAATAAAGGCGAAAGGCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACAGCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGTAG‑TGACAATACCTGCGGCGAGCAGCAGTAAATTCAGCGAGATCGGGTTGTGCCCTAAATGGCGGGTTGAGCt                                                                                                                                                                <  2:577782/251‑1 (MQ=255)
                                                                              tGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGTAG‑TGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCg                                                                                                                                                            <  2:150934/250‑1 (MQ=255)
                                                                              tGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTGTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGTAG‑TGACAATACCGGCGGCGATCAGCAGGAAATTCAGCGACATCGGGTTGTGCCCCGTATGTCGGGTGGGGCTGTCg                                                                                                                                                            <  2:331928/250‑1 (MQ=255)
                                                                                ccAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGGAG‑TGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGGTTTGCCCCATATGGCTGGTTGAGCTGTCGGc                                                                                                                                                          >  1:284260/1‑250 (MQ=255)
                                                                                      ccACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGTAG‑TGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGGTGAACTGTCGGCAATAGc                                                                                                                                                    >  1:408049/1‑250 (MQ=255)
                                                                                          ccAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATCAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAACACAC‑AACAGCTGTACGGTAG‑GGACCATACCGGCGGCGGCCCGCAGTAAATTCAGCGACATCGGGGTTTGCACCATATGGCTGGTTTAGCttttgggcataaacaaa                                                                                                                                                >  2:234249/1‑234 (MQ=255)
                                                                                                        aaGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGTAG‑TGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAAGTgccgcc                                                                                                                                 <  2:329806/251‑1 (MQ=255)
                                                                                                          gTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGTAG‑TGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCAc                                                                                                                               <  1:472282/251‑1 (MQ=255)
                                                                                                             aCCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGTAG‑TGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTTAGCTTTCGGCAATAGCAAACAGGTAAATTGCCGCCACggg                                                                                                                            >  1:194981/1‑251 (MQ=255)
                                                                                                              ccATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAAAAC‑AAAAGCTGTACGGTAG‑TGACAAAACCGGCGGCGAACAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGATGGTTGAGCTGTAGGCACTAGCAAACAGGTAAATTGCCGCCACgggg                                                                                                                           >  1:551175/1‑250 (MQ=255)
                                                                                                                    tATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGTAG‑TGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGcagcag                                                                                                                      <  1:187555/250‑1 (MQ=255)
                                                                                                                    tATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AAAAGCTGTACGGGAG‑GGACAATACCGGGGGCGATCAGCAGTAAATTCAGCGACATCGGGGTTTGCCCCATATGGCTGGGTGAGCTGTCGGGAATAGCAAACAGGGAAATTGCCGGCAAGGGcagcagg                                                                                                                     >  1:231330/1‑250 (MQ=255)
                                                                                                                       cggccctcGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTTGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGTAG‑TGACTATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAGCAGGTAAATTGCCGCCACGGGCAGCAGCca                                                                                                                   <  1:82731/243‑1 (MQ=255)
                                                                                                                                             aaaCGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGCAAAACAC‑AACCGCTGTAAGGTAG‑TGACAAAACCGGCGGCGATCAACAGTAAATTTAACGACATCGGGGTTTGCCCCATAAGGCTGGTTGAGCTGTCGGCAATAGCAAAAAGGGTAATTGCCGCCCCGGGGAGCAGCCACCTGGGTTTGATTAACAT‑GCCgg                                                                                           >  1:405168/1‑251 (MQ=255)
                                                                                                                                                 gTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGTAG‑TGACCATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGGTGAGCTGGCGGCAATAGCAAACAGGGAAATTGCCGCCACCGGCAGCAGCCACATGGGTTTGATTAACAA‑GCCGGTTTg                                                                                       >  1:392389/1‑251 (MQ=255)
                                                                                                                                                 gTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGTAG‑TGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGCGTTTTGCCCCATATGGCTGGTTTAGCTGTCGGCAATAGCAAACAGGTAAATTGCCCCCACGGGCAGCAGC‑ACATGGTTTCGATTAACAT‑GCCGGTTTgc                                                                                      <  1:577937/251‑1 (MQ=255)
                                                                                                                                                 gTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGTAG‑TGACAATAACGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGGAAATTTGCGCCACGGGCAGCCGCCACATGGTTTCGATTAAAAT‑GCCCGtttt                                                                                       >  1:24188/1‑250 (MQ=255)
                                                                                                                                                   cacaGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGTAG‑TGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACAT‑GCCGGTTTgcg                                                                                     <  1:234252/251‑1 (MQ=255)
                                                                                                                                                   cacaGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGTAG‑TGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACAT‑GCCGGTTTgcg                                                                                     <  1:209566/251‑1 (MQ=255)
                                                                                                                                                     cagccagcAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGTAG‑TGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATATCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACAT‑GCCGGTTTGCGC‑t                                                                                  <  2:325766/251‑1 (MQ=255)
                                                                                                                                                                               ccAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑CACAGCTGTACGGTAG‑TGACAATACCGGCGGCGGTCAGCAGGAAATTCCGCGACATCGGGTTTTGCCCCATATGGCCGGTTTAGCTGTCGGCAAAAGCCAACAGGTAAATTTCCGCCACCGGCAACAGCCACCTTGTTTCGATTAACAT‑GCCGGGTTGCGCCTTCA‑CGGCCATCCTCTTGCGtccccc                                                        >  2:142301/1‑246 (MQ=255)
                                                                                                                                                                                 aaTGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAACTGAGCGGGAG‑TGAAAATACCGGCGGCGATCAGCAGTAAATTCAGCGACCTCGGGGTTTGCCCCATATGGCTGGTTGAGCTGTCGGGACTCGCAAACAGGTAAATTGTCGCCAAGGGCAGCAGCCCCATGGGTTCGGTTAAAAT‑GCCGGTTTGCGC‑TTCAAAGGGAAACATTTTGCGTACCAGAc                                                      >  1:303452/1‑251 (MQ=255)
                                                                                                                                                                                  aTGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGAAGCGGTAAAACAA‑AAAAGATGTACGGTAG‑TGACAATACCGGCGGCGAACAGCAGTAAATTTAGCGACCTCGGGTTTTGCCCCATATGGCTGGGTGAGCTGTCGGCAATAGCAAACAGGTAAATTGGCGCCAAGGGCAGCAGCCACATGGGTTCGATTAACAT‑GCCGGGTTGCGC‑CTCAACGGCAATCTTCTTTCGTACCAGAcc                                                     >  1:244092/1‑251 (MQ=255)
                                                                                                                                                                                   tGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGTAG‑TGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGGTTTGCCCCATATGGCTGGTTGAGCTGTCGGCCATAGCAAACAGGGAAATTTCCGCCCCGGGCAGCAGCCACAAGGTTTCGATTAACAT‑GCCGGGTTGCGC‑TTCAACGGCAATCTTCCTGCCTACCAGAccc                                                    >  2:475639/1‑250 (MQ=255)
                                                                                                                                                                                    gTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGTAG‑TGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACAT‑GCCGGTTTGCGC‑TTCAACGGCAATCTTCTTGCGTACCAGACCg                                                    <  2:55908/250‑1 (MQ=255)
                                                                                                                                                                                                  cTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGTAG‑TGACAATACCGGCGGGGATCAGCAGTAAAATCAGCGGCATCGGGTTTTGCCCCATATGGCCGGTTTAGCTTTCGGCAATAACAAACAGGGAAAATGCCGCCCCGGGCAGCCGCCCCCTGGTTTCGATTAACCT‑GCCGGTTTGCGC‑TTCCACCGCCATATTCTTTCGTAACACACCCTAGAAGGCcaaaat                                      >  2:135115/1‑248 (MQ=255)
                                                                                                                                                                                                  cTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGGAG‑TGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTTAGCTGTCGGCAATAGCAAACAGGTAAATTTCCGCCACGGGCAGCAGCCACATGGGTTCTATTAACAT‑ACCGGGTTTCGC‑TTCAAAGGCAATCTTCTTGCGTACCACACCGTAGAAGGCaaaaat                                      >  2:182635/1‑248 (MQ=255)
                                                                                                                                                                                                          gagagaCGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCTGTACGGTAG‑TGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACAT‑GCCGGTTTGCGC‑TTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTAAATGCCa                              <  1:458178/250‑1 (MQ=255)
                                                                                                                                                                                                             ggaCGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCGGTACGGTAG‑TGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACAT‑GCCGGTTTGCGC‑TTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTAAATGCCAGTcc                          <  2:279487/250‑1 (MQ=255)
                                                                                                                                                                                                             agaCGCAAGCGCGAGGCAGCGGCGGTAAAACCCCAAC‑GCCGTACGGAAG‑GGCCAATACCCGCGGCGATCAGCAGTAAATTCAGCGACATCTGGTTTTGGCCCATATGGCCGGGTGAGCTGCCGGCAATAGCAAACAAGGAAATTGCCGCCCCCGGCCGCCACCACAATGTTTCCACTAACCTCGCCGGTT‑GCCC‑CTCCACGGCCAACTTCTTGTCTACCAAACCGTAAAAAGCAAAACTAAATACCaagac                          >  2:423377/1‑247 (MQ=255)
                                                                                                                                                                                                                                      tAAAACAC‑AACAGCTGTACGGTAG‑TGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACAT‑GCCGGTTTGCGC‑TTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAAcc  >  1:383498/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              
CGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAC‑AACAGCGGTACGGTAG‑TGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACAT‑GCCGGTTTGCGC‑TTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACC  >  NC_000913/4003258‑4003737

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: ‑