Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,102,671 G→C 6.8% *263S (TGA→TCA) ‡ phhA → phenylalanine‑4‑hydroxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,102,6710GC 6.8% 284.3 / 7.9 88*263S (TGA→TCA) ‡phhAphenylalanine‑4‑hydroxylase
Reads supporting (aligned to +/- strand):  ref base G (43/39);  new base C (2/4);  total (45/43)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.82e-01

ATCGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGA  >  NC_002947/5102526‑5102814
                                                                                                                                                 |                                                                                                                                                
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                                                                                                                             >  4:377154/1‑138 (MQ=255)
                             aCTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATaacaac                                                                                                                             <  3:377154/138‑1 (MQ=255)
                               tGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTCAGGCAGCGACCTTgg                                                                                                                                       >  8:196678/1‑125 (MQ=255)
                               tGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTCAGGCAGCGACCTTgg                                                                                                                                       <  7:196678/126‑2 (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                                                                                                        <  7:192421/137‑1 (MQ=255)
                                                   ttCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCAt                                                                                                        >  8:192421/1‑137 (MQ=255)
                                                   ttCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGa                                                                                             <  7:32126/148‑1 (MQ=255)
                                                             cccACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAAt                                                                                                    >  6:190234/1‑131 (MQ=255)
                                                             cccACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAAt                                                                                                    <  5:190234/131‑1 (MQ=255)
                                                              ccACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCg                                                                                                                                              >  3:220208/1‑88 (MQ=255)
                                                              ccACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCg                                                                                                                                              <  4:220208/88‑1 (MQ=255)
                                                                  gAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAATTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGa                                                                             <  3:264964/149‑1 (MQ=255)
                                                                  gAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTg                                                                                                                                        >  5:226068/1‑90 (MQ=255)
                                                                  gAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTg                                                                                                                                        <  6:226068/90‑1 (MQ=255)
                                                                   aGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGaa                                                                            <  5:13019/149‑1 (MQ=255)
                                                                    ggACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAg                                                                           >  2:264508/1‑149 (MQ=255)
                                                                          catGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAAcc                                                                                          >  5:175037/1‑128 (MQ=255)
                                                                          catGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAAcc                                                                                          <  6:175037/128‑1 (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                                                    <  5:368610/149‑1 (MQ=255)
                                                                                           aaGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGaaa                                                    <  6:382674/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              >  2:257261/1‑149 (MQ=255)
                                                                                                                                aCCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGaaa              >  3:361782/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)
                                                                                                                                              ccTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACcga                                           <  6:155808/106‑1 (MQ=255)
                                                                                                                                              ccTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACcga                                           >  5:155808/1‑106 (MQ=255)
                                                                                                                                              ccTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGa  <  1:257261/147‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
ATCGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGA  >  NC_002947/5102526‑5102814

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: