Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,102,667 G→C 7.0% A262P (GCC→CCC) ‡ phhA → phenylalanine‑4‑hydroxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,102,6670GC 7.0% 278.2 / 7.2 86A262P (GCC→CCC) ‡phhAphenylalanine‑4‑hydroxylase
Reads supporting (aligned to +/- strand):  ref base G (42/38);  new base C (2/4);  total (44/42)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

CCGCATCGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTG  >  NC_002947/5102522‑5102812
                                                                                                                                                 |                                                                                                                                                  
ccGCATCGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGc                                                                                                                                                   <  3:227957/147‑1 (MQ=255)
    aTCGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGc                                                                                                                                             >  2:31613/1‑149 (MQ=255)
             cTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTCAGGCAGCGACCTTgg                                                                                                                                     <  1:165158/148‑2 (MQ=255)
             cTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTCAGGCAGCGACCTTgg                                                                                                                                     <  7:131939/148‑2 (MQ=255)
              tGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGcc                                                                                                                                    >  2:93209/1‑148 (MQ=255)
               gCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTATGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGa                                                                                                                                  >  7:16228/1‑149 (MQ=255)
                cAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGcc                                                                                                                                    <  3:90249/146‑1 (MQ=255)
                cAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGa                                                                                                                                  <  5:176489/148‑1 (MQ=255)
                      gtgtATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATaacaac                                                                                                                           >  1:133703/1‑149 (MQ=255)
                          aTTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCg                                                                                                                        >  3:188394/1‑148 (MQ=255)
                           tttCGTACTGCCGAACATGAAGCGCCTGATCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGaa                                                                                                                      >  1:23306/1‑149 (MQ=255)
                            ttCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAAc                                                                                                                     >  7:87787/1‑149 (MQ=255)
                              cGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCAAGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCg                                                                                                                   >  6:182336/1‑149 (MQ=255)
                              cGTACGGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGcc                                                                                                                                    >  7:339647/1‑132 (MQ=255)
                              cGTACGGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGcc                                                                                                                                    <  8:339647/132‑1 (MQ=255)
                               gTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCAGAGCGACCTTGCCGATAACAACTCGAACCg                                                                                                                   <  1:31613/148‑1 (MQ=255)
                                 aCTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATaacaac                                                                                                                           <  3:377154/138‑1 (MQ=255)
                                 aCTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATaacaac                                                                                                                           >  4:377154/1‑138 (MQ=255)
                                   tGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTCAGGCAGCGACCTTgg                                                                                                                                     <  7:196678/126‑2 (MQ=255)
                                   tGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTCAGGCAGCGACCTTgg                                                                                                                                     >  8:196678/1‑125 (MQ=255)
                                     ccGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCAATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAaca                                                                                                            >  4:264964/1‑149 (MQ=255)
                                       gAACATGAAGCGCCTGTTAGACCTGGCCCACGAGGACATAATGGGAATGGTCCATAAAGCCATGCAGATGGGTCTGAATGCACCGAAGTTTACACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAacac                                                                                                           >  4:135977/1‑148 (MQ=255)
                                           aTGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCAt                                                                                                      <  2:23306/149‑1 (MQ=255)
                                                 cgcCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGc                                                                                                <  1:346553/149‑1 (MQ=255)
                                                   ccTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCt                                                                                              >  6:13019/1‑149 (MQ=255)
                                                    cTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCtt                                                                                             <  4:188394/149‑1 (MQ=255)
                                                       ttCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCAt                                                                                                      >  8:192421/1‑137 (MQ=255)
                                                       ttCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCAt                                                                                                      <  7:192421/137‑1 (MQ=255)
                                                       ttCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGa                                                                                           <  7:32126/148‑1 (MQ=255)
                                                                 cccACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAAt                                                                                                  <  5:190234/131‑1 (MQ=255)
                                                                 cccACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAAt                                                                                                  >  6:190234/1‑131 (MQ=255)
                                                                  ccACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCg                                                                                                                                            <  4:220208/88‑1 (MQ=255)
                                                                  ccACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCg                                                                                                                                            >  3:220208/1‑88 (MQ=255)
                                                                      gAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAATTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGa                                                                           <  3:264964/149‑1 (MQ=255)
                                                                      gAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTg                                                                                                                                      <  6:226068/90‑1 (MQ=255)
                                                                      gAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTg                                                                                                                                      >  5:226068/1‑90 (MQ=255)
                                                                       aGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGaa                                                                          <  5:13019/149‑1 (MQ=255)
                                                                        ggACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAg                                                                         >  2:264508/1‑149 (MQ=255)
                                                                              catGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAAcc                                                                                        <  6:175037/128‑1 (MQ=255)
                                                                              catGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAAcc                                                                                        >  5:175037/1‑128 (MQ=255)
                                                                                  ggCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCgcgc                                                               <  8:87787/149‑1 (MQ=255)
                                                                                    cATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTc                                                                                                        <  2:324549/106‑1 (MQ=255)
                                                                                    cATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTc                                                                                                        >  1:324549/1‑106 (MQ=255)
                                                                                        gTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCTGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCTATgc                                                         <  3:200078/149‑1 (MQ=255)
                                                                                         tCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATgcg                                                        >  7:43900/1‑149 (MQ=255)
                                                                                          ccATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATgcgc                                                       >  8:176784/1‑149 (MQ=255)
                                                                                           cATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACAGCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGcc                                                      >  5:382674/1‑149 (MQ=255)
                                                                                             tAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCt                                                                      <  1:219826/131‑1 (MQ=255)
                                                                                             tAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCt                                                                      >  2:219826/1‑131 (MQ=255)
                                                                                             tAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCg                                                     >  6:368610/1‑148 (MQ=255)
                                                                                               aaGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGaaa                                                  <  6:382674/149‑1 (MQ=255)
                                                                                               aaGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGaaa                                                  <  5:368610/149‑1 (MQ=255)
                                                                                                  ccATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGt                                               >  8:110166/1‑148 (MQ=255)
                                                                                                  ccATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTc                                              <  1:93209/149‑1 (MQ=255)
                                                                                                      gCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTCAGGCAGCGACCTTg                                                                                                                                      <  7:236260/58‑1 (MQ=21)
                                                                                                      gCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTCAGGCAGCGACCTTg                                                                                                                                      >  8:236260/1‑58 (MQ=21)
                                                                                                         gCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACcgac                                        >  7:354723/1‑148 (MQ=255)
                                                                                                            gggTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGCACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCAccgccg                                       >  7:16912/1‑146 (MQ=255)
                                                                                                                   catgcaCCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATgc                                                         <  2:146396/122‑1 (MQ=255)
                                                                                                                   catgcaCCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATgc                                                         >  1:146396/1‑122 (MQ=255)
                                                                                                                    atgcaCCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGcc                            >  4:173043/1‑149 (MQ=255)
                                                                                                                      gcaCCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGTAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGc                                                                        <  5:141904/104‑1 (MQ=255)
                                                                                                                      gcaCCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGc                                                                        >  6:141904/1‑104 (MQ=255)
                                                                                                                        aCCGAAGTTTCCACCCAACGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGACGAGCTGGCCGAGc                        <  8:233011/149‑1 (MQ=255)
                                                                                                                          cGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCt                       >  4:165012/1‑148 (MQ=255)
                                                                                                                           gAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGc                        >  4:310000/1‑146 (MQ=255)
                                                                                                                           gAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGc                        <  3:310000/146‑1 (MQ=255)
                                                                                                                                ttCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCcgc                >  6:352429/1‑149 (MQ=255)
                                                                                                                                  ccACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATGcgcg               >  7:31846/1‑148 (MQ=255)
                                                                                                                                   cACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGTAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGa              <  5:182336/148‑1 (MQ=255)
                                                                                                                                    aCCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGaaa            >  3:361782/1‑149 (MQ=255)
                                                                                                                                    aCCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGaaa            >  2:257261/1‑149 (MQ=255)
                                                                                                                                      ccAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGAt                                                           <  7:208616/101‑1 (MQ=255)
                                                                                                                                      ccAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGAt                                                           >  8:208616/1‑101 (MQ=255)
                                                                                                                                      ccAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGACTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATc          >  6:315891/1‑149 (MQ=255)
                                                                                                                                       cAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATcc         <  7:110166/149‑1 (MQ=255)
                                                                                                                                       cAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATcc         >  3:82310/1‑149 (MQ=255)
                                                                                                                                        aaGGTCTCTGCCTGAGCGACCTTGCCGATAACAAATCGAACCTGAAAACACCTCATGAATGCCTTGAACGAAGCCAATTGCGAAGCCTTCCTCGCCGATGCGCCGAAAGTCACGGACGAAGAGCTGGCCGAGCTGATCAGCGAAATccc        <  3:135977/149‑1 (MQ=255)
                                                                                                                                        aaGGTCGCTGCCTGAGCGACCTTGCCGATATCTACTCGAACCGGAAAACTCGTCTTGAATGGCTTGAACCAAGCGCATTGCGAAGCCTGCCGGGCCGATGCGCCGAAAGTCACCGACGAAGAGCTTGCCGAGCTGATCCGCGAAATccc        <  8:16912/149‑1 (MQ=255)
                                                                                                                                        aaGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATccc        >  5:207536/1‑149 (MQ=255)
                                                                                                                                         aGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCg       <  8:43900/149‑1 (MQ=255)
                                                                                                                                           gTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACTAAGAGCTGGCCGAGCTGATCCGCGAAATCCCg       <  1:237929/147‑1 (MQ=255)
                                                                                                                                           gTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCg       >  2:237929/1‑147 (MQ=255)
                                                                                                                                           gTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGa     <  8:16228/149‑1 (MQ=255)
                                                                                                                                            tCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGa     <  1:264508/149‑1 (MQ=255)
                                                                                                                                              gCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTg  <  3:173043/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CCGCATCGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTG  >  NC_002947/5102522‑5102812

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