Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,037,961 G→A 6.6% A478A (GCC→GCT PP_0897 ← class 1 fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,037,9610GA 6.6% 229.2 / 3.9 76A478A (GCC→GCTPP_0897class 1 fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (39/32);  new base A (2/3);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.34e-01

TGCGCTGTACGTGTGGGAGCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAGGTACACGGCCTTGTTGT  >  NC_002947/1037813‑1038109
                                                                                                                                                    |                                                                                                                                                     
tGCGCTGTACGTGTGGGAGCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACg                                                                                                                                                       >  4:245073/1‑149 (MQ=255)
  cgcTGTACGTGTGGGAGCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACgg                                                                                                                                                      <  4:102011/148‑1 (MQ=255)
    cTGTACGTGTGGGAGCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGa                                                                                                                                                   >  6:218664/1‑149 (MQ=255)
    cTGTACGTGTGGGAGCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGa                                                                                                                                                   >  2:207384/1‑149 (MQ=255)
         cGTGTGGGAGCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTCGCCGTGGACAGCAACGGt                                                                                                                                               <  6:354304/148‑1 (MQ=255)
              gggAGCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTCGCCGTGGACAGCAACGGt                                                                                                                                               >  4:221415/1‑143 (MQ=255)
              gggAGCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTCGCCGTGGACAGCAACGGt                                                                                                                                               <  3:221415/143‑1 (MQ=255)
                       cTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCtt                                                                                                                                >  8:207937/1‑149 (MQ=255)
                       cTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCtt                                                                                                                                <  7:65422/149‑1 (MQ=255)
                         tGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGa                                                                                                                              >  7:49382/1‑149 (MQ=255)
                            ccGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACtt                                                                                                                           <  7:77529/149‑1 (MQ=255)
                             cgcgAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGCCGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCCCGGCGACGGTCCCCGGCATGTCCTGGCCTTc                                                                                                                          >  4:223512/1‑149 (MQ=255)
                             cgcgAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGTCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTc                                                                                                                          >  4:353771/1‑149 (MQ=255)
                              gcgAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCGTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCg                                                                                                                         <  1:324644/149‑1 (MQ=255)
                                   gcagcCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTc                                                                                                                    <  2:79575/149‑1 (MQ=255)
                                   gcagcCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTc                                                                                                                    <  6:319570/149‑1 (MQ=255)
                                   gcagcCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTCGCCGTGGACAGCAACgg                                                                                                                                                >  6:128738/1‑121 (MQ=255)
                                   gcagcCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTCGCCGTGGACAGCAACgg                                                                                                                                                <  5:128738/121‑1 (MQ=255)
                                     agcCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGTCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCg                                                                                                                   >  2:219083/1‑148 (MQ=255)
                                      gccgccGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGt                                                                                                                                               <  8:88606/119‑1 (MQ=255)
                                      gccgccGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGt                                                                                                                                               >  7:88606/1‑119 (MQ=255)
                                            gcAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTCATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGc                                                                                                           >  7:35643/1‑149 (MQ=255)
                                                    aTCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGc                                                                                                   >  5:330851/1‑149 (MQ=255)
                                                         gTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGc                                                                                              >  6:153403/1‑149 (MQ=255)
                                                             cttcacttcGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCgg                                                                                          >  1:203196/1‑149 (MQ=255)
                                                              ttcacttcGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGc                                                                                         >  2:179637/1‑149 (MQ=255)
                                                                cacttcGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAATGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGc                                                                                         <  6:163051/147‑1 (MQ=255)
                                                                 acttcGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGaa                                                                                      <  7:207937/149‑1 (MQ=255)
                                                                  cttcGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAg                                                                                     >  4:170929/1‑149 (MQ=255)
                                                                   ttcGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAgg                                                                                   >  4:353673/1‑149 (MQ=255)
                                                                     cGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGc                                                                                  >  5:199694/1‑148 (MQ=255)
                                                                      gACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGc                                                                                  >  5:14723/1‑147 (MQ=255)
                                                                      gACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCa                                                                                <  6:199694/149‑1 (MQ=255)
                                                                          gCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGATGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTAc                                                                            >  8:21804/1‑149 (MQ=255)
                                                                           ccAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGc                                                                                  <  6:185570/142‑1 (MQ=255)
                                                                           ccAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGc                                                                                  >  5:185570/1‑142 (MQ=255)
                                                                           ccAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTAcc                                                                           >  4:89308/1‑149 (MQ=255)
                                                                             aGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCtt                                                                         >  4:317461/1‑149 (MQ=255)
                                                                               gCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTc                                                                        >  2:88938/1‑148 (MQ=255)
                                                                                cTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGTCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGa                                                                      <  1:219083/149‑1 (MQ=255)
                                                                                 tttCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGTGGCAGAACCTTGGa                                                                      >  1:127333/1‑149 (MQ=255)
                                                                                     ggCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGc                                                                 <  3:245073/149‑1 (MQ=255)
                                                                                       cAATCTAGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGc                                                                                                           >  7:293934/1‑106 (MQ=255)
                                                                                       cAATCTAGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGc                                                                                                           <  8:293934/106‑1 (MQ=255)
                                                                                         aTCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGAt                                                             <  3:353673/149‑1 (MQ=255)
                                                                                              gCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGTAGGCCAGTACCTTCGACTTGCGGATGGc                                                          <  1:207384/147‑1 (MQ=255)
                                                                                              gCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATg                                                            <  1:88938/145‑1 (MQ=255)
                                                                                              gCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGc                                                          <  5:218664/147‑1 (MQ=255)
                                                                                              gCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCt                                                        >  6:318628/1‑149 (MQ=255)
                                                                                                tctGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTgg                                                      >  2:352988/1‑149 (MQ=255)
                                                                                                tctGCCACAGGGCAGGGCCAGTGATGTGTACCGATTAACAGTTGCTGTCCAAGGCGACGGTCACCGGAATGTACATGACTTAGAACTCGAAGAAAGCATCCAAGCCCAGAACGGAGAACGCCAGAAAAGTAGACTTGCGGaagtccac                                                       >  8:348856/1‑141 (MQ=255)
                                                                                                 ctGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTggg                                                     >  7:41233/1‑149 (MQ=255)
                                                                                                  tGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCTGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTgg                                                      <  8:35643/147‑1 (MQ=255)
                                                                                                    ccACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTggg                                                     <  8:41233/146‑1 (MQ=255)
                                                                                                    ccACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGc                                                    <  3:89308/147‑1 (MQ=255)
                                                                                                    ccACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGc                                                    <  1:179637/147‑1 (MQ=255)
                                                                                                     cacaGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGAc                                                 >  7:84497/1‑149 (MQ=255)
                                                                                                     cacaGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGAc                                                 >  3:358048/1‑149 (MQ=255)
                                                                                                      acaGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTACTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGAc                                                 >  5:10410/1‑148 (MQ=255)
                                                                                                       cagggcagggCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGAc                                                 >  8:162516/1‑147 (MQ=255)
                                                                                                        agggcagggcCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAg                                              <  7:21804/149‑1 (MQ=255)
                                                                                                         gggcagggcCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAgg                                             <  8:49382/149‑1 (MQ=255)
                                                                                                            cagggcCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAg                                          <  5:318628/149‑1 (MQ=255)
                                                                                                                    gTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGcgcc                                  <  6:330851/149‑1 (MQ=255)
                                                                                                                      gATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGcgccgc                                >  4:320048/1‑149 (MQ=255)
                                                                                                                           gtACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACg                           >  7:176919/1‑149 (MQ=255)
                                                                                                                              ccGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGcc                        <  2:203196/149‑1 (MQ=255)
                                                                                                                                gATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCAt                      >  5:15603/1‑149 (MQ=255)
                                                                                                                                  ttCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCa                    <  3:353771/149‑1 (MQ=255)
                                                                                                                                    cACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAg                   <  3:320048/148‑1 (MQ=255)
                                                                                                                                     aCCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAgg                  >  6:237264/1‑148 (MQ=255)
                                                                                                                                     aCCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAgg                  >  3:65428/1‑148 (MQ=255)
                                                                                                                                             tGTCAACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAGGTACACGGcc         >  2:85070/1‑149 (MQ=255)
                                                                                                                                                 cACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAGGTACACGGCCttg      <  5:12208/148‑1 (MQ=255)
                                                                                                                                                 cACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAGGTACACGGCCttg      <  8:84497/148‑1 (MQ=255)
                                                                                                                                                  aCGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAGGTACACGGCCttgtt    >  8:345075/1‑149 (MQ=255)
                                                                                                                                                  aCGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAGGTACACGGCCttgtt    <  3:223512/149‑1 (MQ=255)
                                                                                                                                                   cGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAGGTACACGGCCttgttg   <  6:10410/149‑1 (MQ=255)
                                                                                                                                                    ggCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAGGTACACGGCCTtgttgt  >  5:85373/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGCGCTGTACGTGTGGGAGCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAGGTACACGGCCTTGTTGT  >  NC_002947/1037813‑1038109

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: