Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,522,243 C→A 6.2% T31N (ACC→AAC)  idh → isocitrate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,522,2430CA 6.2% 303.3 / 7.1 96T31N (ACC→AAC) idhisocitrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (42/46);  new base A (3/3);  total (45/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GCCTACCCTGCGCACCACGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTACCACCCCTGAAGCCAAC  >  NC_002947/4522095‑4522391
                                                                                                                                                    |                                                                                                                                                    
gcCTACCCTGCGCACCACGAGACGCGATCACGCTCAGCACCCAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGGGTTCAc                                                                                                                                                      <  4:209152/149‑1 (MQ=255)
          cgcACCACGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGc                                                                                                                                             <  1:169531/148‑1 (MQ=255)
          cgcACCACGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCt                                                                                                                                            <  3:293966/149‑1 (MQ=255)
          cgcACCACGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGGTGAACGCTTCga                                                                                                                                              <  6:201609/147‑2 (MQ=255)
           gcACCACGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTg                                                                                                                                           >  1:251822/1‑149 (MQ=255)
            caccacGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGa                                                                                                                                          <  4:283881/149‑1 (MQ=255)
               cacGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACAt                                                                                                                                       <  3:143785/149‑1 (MQ=255)
               cacGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACAt                                                                                                                                       >  3:324810/1‑149 (MQ=255)
                     aCTCGAGCACTCGCAGAACACAGAGAGTTAACCCTCCTGCCCACACGTGCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGt                                                                                                                                 <  4:221149/149‑1 (MQ=255)
                     aCTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCCTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGt                                                                                                                                 <  4:341712/149‑1 (MQ=255)
                     aCTCGAGCACGCTCAGCACACAGAGAGTTAACCCCCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGt                                                                                                                                 <  2:368720/149‑1 (MQ=255)
                       tCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTc                                                                                                                                >  7:195728/1‑148 (MQ=255)
                         gAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCCCTGTCGa                                                                                                                              >  3:140532/1‑148 (MQ=255)
                          aGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGaaa                                                                                                                            <  6:394481/149‑1 (MQ=255)
                               gCTCAGCACACAGAGAGGTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTCCTCGCTGCTGCCGATCATCGACGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTcgc                                                                                                                       <  6:137402/149‑1 (MQ=255)
                                 tCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGa                                                                                                                     >  7:203460/1‑149 (MQ=255)
                                   aGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACa                                                                                                                   >  3:62460/1‑149 (MQ=255)
                                      acacaGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACAttt                                                                                                                <  4:324810/149‑1 (MQ=255)
                                      acacaGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACAttt                                                                                                                <  5:303111/149‑1 (MQ=255)
                                        acaGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTcc                                                                                                              <  4:62460/149‑1 (MQ=255)
                                           gagagTTAACCCGCATGCCCACCCGTTCAAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCctg                                                                                                           >  1:23481/1‑149 (MQ=255)
                                                ttAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCctggctgg                                                                                                      >  1:366619/1‑149 (MQ=255)
                                                ttAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCctggctgg                                                                                                      >  5:120376/1‑149 (MQ=255)
                                                 tAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACa                                                                                                                   <  8:3501/135‑1 (MQ=255)
                                                 tAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACa                                                                                                                   >  7:3501/1‑135 (MQ=255)
                                                   aCCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGcc                                                                                                    <  2:251822/148‑1 (MQ=255)
                                                    cccGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCg                                                                                                   >  3:373130/1‑148 (MQ=255)
                                                       gCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATc                                                                                               >  5:113305/1‑149 (MQ=255)
                                                           gCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGcc                                                                                                    >  5:213376/1‑140 (MQ=255)
                                                           gCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGcc                                                                                                    <  6:213376/140‑1 (MQ=255)
                                                           gCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTg                                                                                           >  8:233648/1‑149 (MQ=255)
                                                             ccACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTgcc                                                                                         <  7:62594/149‑1 (MQ=255)
                                                             ccACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGGTGAACGCTTCga                                                                                                                                              >  2:172812/1‑95 (MQ=255)
                                                             ccACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGGTGAACGCTTCga                                                                                                                                              <  1:172812/96‑2 (MQ=255)
                                                              gctgtgTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTgccg                                                                                        <  4:336475/144‑1 (MQ=255)
                                                                 ccGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCct                                                                                                            >  6:117549/1‑126 (MQ=255)
                                                                 ccGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCt                                                                                     <  7:233648/149‑1 (MQ=255)
                                                                 ccGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCt                                                                                     >  1:113675/1‑149 (MQ=255)
                                                                        aaGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGa                                                                              <  2:366619/149‑1 (MQ=255)
                                                                        aaGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGa                                                                              <  3:262337/149‑1 (MQ=255)
                                                                         atatcatgTATAGCTTGACCGACGAAGGCCCGGCCGTCGGCACTTACTCGCTGCTGGCGATCATCGAAGCGTTCTGCGCTTCGGCTGACTTCGCTGTCGAAACTCGCGACATTTGCCTGGCTGGCCGTATCCTTGCCGGCTTCGCGGAg                                                                             <  4:115847/147‑1 (MQ=255)
                                                                              atcTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACt                                                                        <  5:326642/149‑1 (MQ=255)
                                                                                 tataCCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTggg                                                                     >  6:132048/1‑149 (MQ=255)
                                                                                  ataCCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTc                                                                                   >  1:249825/1‑134 (MQ=255)
                                                                                  ataCCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTc                                                                                   <  2:249825/134‑1 (MQ=255)
                                                                                  ataCCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGgc                                                                    >  4:91454/1‑149 (MQ=255)
                                                                                   taCCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTgg                                                                      <  6:120376/146‑1 (MQ=255)
                                                                                   taCCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTgg                                                                      <  4:140532/146‑1 (MQ=255)
                                                                                      cTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCg                                                                                <  6:137886/133‑1 (MQ=255)
                                                                                      cTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGgcg                                                                   >  3:71891/1‑146 (MQ=255)
                                                                                      cTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGgcg                                                                   >  5:306136/1‑146 (MQ=255)
                                                                                      cTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCa                                                                 <  4:71891/148‑1 (MQ=255)
                                                                                      cTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTACCGCCTTCCCGGAGCAACTGGgcg                                                                   >  1:92289/1‑146 (MQ=255)
                                                                                          aCCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCagag                                                              >  5:34813/1‑147 (MQ=255)
                                                                                           ccGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAg                                                           >  2:21613/1‑149 (MQ=255)
                                                                                             gacgaAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATc                                                                                               >  6:185969/1‑111 (MQ=255)
                                                                                             gacgaAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATc                                                                                               <  5:185969/111‑1 (MQ=255)
                                                                                              tcgaAGCCCACGCCGTCGCCAGTTACTCGCTGCTGCCTATCATCGAAGCGTTCACCGCTTCGGCTTACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTTTGGTTGCCGCCTTCCCGGAGCAACTGGGCGGAGAGAAGc                                                          <  5:190188/146‑1 (MQ=255)
                                                                                              acgaAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTACCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCaa                                                        <  2:92289/149‑1 (MQ=255)
                                                                                                gaAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGt                                                      <  6:113305/149‑1 (MQ=255)
                                                                                                gaAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGt                                                      <  3:91454/149‑1 (MQ=255)
                                                                                                      cccGCCCTCGCCACTTACTCGCTGCTTGCGATCATGGAAGCGTTCAGCGCTTCGGCTGACTTCGCTGTGGAAACTCGCGAGATTTCGCTGTCTGGCAGTTTGCTTGCCGCCTTACCGGAGCAACTTGGCGCAGAGAAGCATGTTGGCGa                                                <  8:203460/149‑1 (MQ=255)
                                                                                                      cccGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGa                                                >  7:248051/1‑149 (MQ=255)
                                                                                                      cccGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCGTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCGTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGa                                                >  1:300407/1‑149 (MQ=255)
                                                                                                           ccTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCAcc                                           >  5:78887/1‑149 (MQ=255)
                                                                                                            cTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTCTCGAAACTCGCGACATTTCACTGGCTGGCCGTATACTTGCCGCCTTACCGGAGAAAATGGGCGCAGAGAAGCAAGTAGGCGATCAACt                                          >  1:111445/1‑149 (MQ=255)
                                                                                                             tCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTg                                         >  5:301740/1‑149 (MQ=255)
                                                                                                               gCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTTGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGc                                       <  3:6596/149‑1 (MQ=255)
                                                                                                               gCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGc                                       <  1:239694/149‑1 (MQ=255)
                                                                                                                 cACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGa                                                <  2:96963/138‑1 (MQ=255)
                                                                                                                 cACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGa                                                >  1:96963/1‑138 (MQ=255)
                                                                                                                    ttACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAAc                                  >  2:318208/1‑149 (MQ=255)
                                                                                                                       cTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTgg                               <  1:318208/149‑1 (MQ=255)
                                                                                                                          gctgctGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCg                                                 >  4:4272/1‑128 (MQ=255)
                                                                                                                          gctgctGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCg                                                 <  3:4272/128‑1 (MQ=255)
                                                                                                                          gctgctGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGa                                                <  4:274210/129‑1 (MQ=255)
                                                                                                                          gctgctGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGa                                                >  3:274210/1‑129 (MQ=255)
                                                                                                                          gctgctGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTggg                              <  4:373130/147‑1 (MQ=255)
                                                                                                                                gcCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTgg                      >  8:216181/1‑149 (MQ=255)
                                                                                                                                 ccGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTgg                      <  6:34813/148‑1 (MQ=255)
                                                                                                                                 ccGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGc                     <  5:132048/149‑1 (MQ=255)
                                                                                                                                 ccGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGc                     >  3:235600/1‑149 (MQ=255)
                                                                                                                                 ccGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGc                     >  7:367257/1‑149 (MQ=255)
                                                                                                                                    tcagccGAAGCGGTGAACGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCag                                                                <  4:8507/98‑1 (MQ=255)
                                                                                                                                    tcagccGAAGCGGTGAACGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCag                                                                >  3:8507/6‑103 (MQ=255)
                                                                                                                                        ccGAAGCGGTGAACGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTCGGCCAGCTGGCTaccgcc              >  5:348316/2‑149 (MQ=255)
                                                                                                                                          gAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTACCAccc             <  7:142888/148‑1 (MQ=255)
                                                                                                                                          gAAGCGTGCACCGCTTCGGCTGACATCGCTGGCGCAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTACCAccc             <  2:300407/148‑1 (MQ=255)
                                                                                                                                           aaGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTACCACCCCt           <  1:21613/149‑1 (MQ=255)
                                                                                                                                           aaGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACAGGGCGCAGAGAAGCAAGTAGGCGATCACCTGCCGGAACTCGGCCAGCTGGCTACCACACCt           >  4:64494/1‑149 (MQ=255)
                                                                                                                                            aGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTgg                      >  5:317082/1‑137 (MQ=255)
                                                                                                                                            aGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTgg                      <  6:317082/137‑1 (MQ=255)
                                                                                                                                              cGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTACCACCCCTGaa        <  6:78887/149‑1 (MQ=255)
                                                                                                                                              cGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCCGGCGGAACTGGGCCAGCTGGCTACCACCCCTGaa        <  6:306136/149‑1 (MQ=255)
                                                                                                                                               gTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGAGGAACTGGGCCAGCTGGCTACCACCCCTGAAg       >  5:106360/1‑149 (MQ=255)
                                                                                                                                                ttCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTACCACCCCTGAAGc      >  4:376113/1‑149 (MQ=255)
                                                                                                                                                   aCCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTACCACCCCTGAAGCCaa   >  2:289193/1‑149 (MQ=255)
                                                                                                                                                    ccGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTACCACCCCTGAAGCCAAc  <  2:113675/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCCTACCCTGCGCACCACGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTACCACCCCTGAAGCCAAC  >  NC_002947/4522095‑4522391

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