Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,464,290 T→C 7.6% L251P (CTT→CCT)  paaX → transcriptional repressor of phenylacetic acid degradation paa operon, phenylacetyl‑CoA inducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,464,2900TC 7.6% 154.6 / 2.9 58L251P (CTT→CCT) paaXtranscriptional repressor of phenylacetic acid degradation paa operon, phenylacetyl‑CoA inducer
Reads supporting (aligned to +/- strand):  ref base T (23/30);  new base C (2/3);  total (25/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGTCATTCCGCCCGCTGGTGCCGCTTTTAAAAGAGGCGGCAGACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGCCTGTATTTTCAACGT  >  NC_000913/1464147‑1464438
                                                                                                                                               |                                                                                                                                                    
tCAGTCATTCCGCCCGCTGGTGCCGCTTTTAAAAGAGGCGGCAGACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCgg                                                                                                                                                 <  2:105041/149‑1 (MQ=255)
tCAGTCATTCCGCCCGCTGGTGCCGCTTTTAAAAGAGGCGGCAGACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCgg                                                                                                                                                 <  1:662599/149‑1 (MQ=255)
   gTCATTCCGCCCGCTGGTGCCGCTTTTAAAAGAGGCGGCAGACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGcac                                                                                                                                              >  2:690657/1‑149 (MQ=255)
       ttCCGCCCGCTGGTGCCGCTTTTAAAAGAGGCGGCAGACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTg                                                                                                                                          <  1:140763/149‑1 (MQ=255)
       ttCCGCCCGCTGGTGCCGCTTTTAAAAGAGGCGGCAGACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTg                                                                                                                                          <  1:621808/149‑1 (MQ=255)
                    tGCCGCTTTTAAAAGAGGCGGCAGACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACg                                                                                                                             <  1:248341/149‑1 (MQ=255)
                      ccGCTTTTAAAAGAGGCGGCAGACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGgc                                                                                                                           >  1:709637/1‑149 (MQ=255)
                      ccGCTTTTAAAAGAGGCGGCAGACGAGTTAACACAGGAGTGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGgc                                                                                                                           >  2:718157/1‑149 (MQ=255)
                        gCTTTTAAAAGAGGCGGCAGACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGCGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGgcgc                                                                                                                         >  1:485494/1‑149 (MQ=255)
                        gCTTTTAAAAGAGGCGGCAGACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGgcgc                                                                                                                         <  1:540931/149‑1 (MQ=255)
                        gCTTTTAAAAGAGGCGGCAGACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGgcgc                                                                                                                         >  1:485501/1‑149 (MQ=255)
                        gCTTTTAAAAGAGGCGGCAGACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCCTCCGGCACACTGGGCAGGGCATACGgcgc                                                                                                                         >  1:172364/1‑149 (MQ=255)
                        gCTTTTAAAAGAGGCGGCAGACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCCTGTCGACCTTAAAGACCCATTGTTGCCGGAGGAGTTGCCTCCGGCACACTGGGCAGGGCATACGgcgc                                                                                                                         >  1:376495/1‑149 (MQ=255)
                         cTTTTAAAAGAGGCGGCAGACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGcgcg                                                                                                                        >  1:152241/1‑149 (MQ=255)
                                  gagGCGGCAGACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCtgtg                                                                                                               <  2:320527/149‑1 (MQ=255)
                                   agGCGGCAGACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTgtgt                                                                                                              <  2:124940/149‑1 (MQ=255)
                                   agGCGGCAGACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTgtgt                                                                                                              <  1:443974/149‑1 (MQ=255)
                                     gcggcAGACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTCCTTCAGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTAt                                                                                                            <  1:188587/149‑1 (MQ=255)
                                      cggcAGACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTGTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGCCAGGGCATACGGCGCGTCAGCTGTGTATc                                                                                                           <  1:494961/149‑1 (MQ=255)
                                           gACGAGTTGACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACAt                                                                                                      <  1:225357/149‑1 (MQ=255)
                                           gACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACAt                                                                                                      <  1:630526/149‑1 (MQ=255)
                                           gACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACAt                                                                                                      <  1:243550/149‑1 (MQ=255)
                                             cGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCGGTGTATCAACAttt                                                                                                    >  1:456437/1‑149 (MQ=255)
                                                  tacccccGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAg                                                                                               <  1:373436/146‑1 (MQ=255)
                                                     ccccGGAGCGGGCATTTCATATTCAGCTTTTACAGATCCATTTTTATAGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAgcgc                                                                                            >  1:504314/1‑149 (MQ=255)
                                                          gggCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGGCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAgc                                                                                       <  2:666678/147‑1 (MQ=255)
                                                                  aTTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCctgctg                                                                               <  2:490378/149‑1 (MQ=255)
                                                                   tttCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCTTTACAGACCCATTGTTGCCGGAGTAGTTTCTTCCTGCACTCTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCCtgctgc                                                                              <  2:210467/149‑1 (MQ=255)
                                                                   tttCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCCtgctgc                                                                              <  2:510815/149‑1 (MQ=255)
                                                                            cAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGcgtt                                                                     <  1:718157/149‑1 (MQ=255)
                                                                                        aTCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGaaaa                                                         >  1:302320/1‑149 (MQ=255)
                                                                                         accATTTTTAGCGCCGTGCCGGCCCTAAAGACCCATTGTTGCCGGAGGAGTTTCCTCCGGCACACTGGGCAGGGCATACGGCGCGTCATCTGTTTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGaaaaa                                                        <  1:713927/148‑1 (MQ=255)
                                                                                          ccATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACCCTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAg                                                       >  1:338615/1‑149 (MQ=255)
                                                                                          ccATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGCTCGTTAGTGAAAAAg                                                       >  1:367179/1‑149 (MQ=255)
                                                                                          ccATTTTTATCGCCGTGTCGTCCTTAAAGACAAATACCTGACGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAg                                                       >  1:775361/1‑149 (MQ=255)
                                                                                           cATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAgg                                                      >  2:729417/1‑149 (MQ=255)
                                                                                           cATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAgg                                                      <  1:48310/149‑1 (MQ=255)
                                                                                                ttATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGaaa                                                 <  2:257107/149‑1 (MQ=255)
                                                                                                 tgtCGCCGTTTCGCCCTTAAAGGCCCATTGTTGCCGGAGGAGTTGCCTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAAc                                                <  2:172364/147‑1 (MQ=255)
                                                                                                 tATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCCTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAAc                                                <  2:376495/149‑1 (MQ=255)
                                                                                                     gCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCtcg                                            <  1:690657/149‑1 (MQ=255)
                                                                                                                 ttAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAActg                                >  1:498407/1‑149 (MQ=255)
                                                                                                                         ccATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGTTGAACTGCCTGCGcc                        >  1:253083/1‑149 (MQ=255)
                                                                                                                          cATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGAGCGTCAGCTGTGTATCAATATTTATCAGCGCGTAGCGCCTGCTGCTTTATCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCg                       >  1:473505/1‑149 (MQ=255)
                                                                                                                          cATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGAGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGGGCCg                       >  1:268249/1‑149 (MQ=255)
                                                                                                                           aTTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGATGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTGCGCTAGGGAAAAAGGTGAAAACTCGGTCGGTGAACTGCCTGCGCCgg                      >  1:449123/1‑149 (MQ=255)
                                                                                                                               ttgCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTGAGCAGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCCTTAGCGTTCGTTAGTGAAAAAGGTGAAACCGCGGTCGCTTGACAGCCTGCGCCGGGAAg                  >  2:302722/1‑149 (MQ=255)
                                                                                                                               ttgCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGGAg                  >  1:752577/1‑149 (MQ=255)
                                                                                                                               ttgCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAg                  >  1:637989/1‑149 (MQ=255)
                                                                                                                                    ggaggcGTTGCTTCAGGCACACTGGGGAGGGCATAAGGGGCATTAGCGGTTTACCAACATCTGTCAACGCGGGGCGGCTTCTGCTTTGGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGCCTGt             <  2:456437/149‑1 (MQ=255)
                                                                                                                                    ggaggaGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGCCTGt             >  1:213884/1‑149 (MQ=255)
                                                                                                                                    ggaggaGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCCCCGGGAAGCCTGt             >  2:279317/1‑149 (MQ=255)
                                                                                                                                     gaggagTTGCTGCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGCCTGTa            <  1:279939/149‑1 (MQ=255)
                                                                                                                                        ggggtGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGGGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTGTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGCCTGTAttt         <  2:341998/145‑1 (MQ=255)
                                                                                                                                        gatgtGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCATCTGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGCCTGTAttt         <  1:364357/145‑1 (MQ=255)
                                                                                                                                        gagTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGCCTGTAttt         <  1:220243/149‑1 (MQ=255)
                                                                                                                                        gagTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGCCTGTAttt         <  1:732650/149‑1 (MQ=255)
                                                                                                                                        gagTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGCCTGAAttt         <  2:523024/149‑1 (MQ=255)
                                                                                                                                               ttCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGCCGGTGAACTGCCTGCGCCGGGAAGCCTGTATTTTCAACGt  >  1:726616/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TCAGTCATTCCGCCCGCTGGTGCCGCTTTTAAAAGAGGCGGCAGACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGCCTGTATTTTCAACGT  >  NC_000913/1464147‑1464438

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