Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,037,958 C→G 6.8% V479V (GTG→GTC)  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,9580CG 6.8% 216.1 / 7.0 74V479V (GTG→GTC) PP_0897class 1 fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (39/29);  new base G (2/3);  total (42/32)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.27e-01

TTGCGCTGTACGTGTGGGAGCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAGGTACACGGCCTTG  >  NC_002947/1037812‑1038105
                                                                                                                                                  |                                                                                                                                                    
ttGCGCTGTACGTGTGGGAGCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCAc                                                                                                                                                    >  1:79575/1‑149 (MQ=255)
 tGCGCTGTACGTGTGGGAGCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACg                                                                                                                                                   >  4:245073/1‑149 (MQ=255)
   cgcTGTACGTGTGGGAGCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACgg                                                                                                                                                  <  4:102011/148‑1 (MQ=255)
     cTGTACGTGTGGGAGCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGa                                                                                                                                               >  2:207384/1‑149 (MQ=255)
     cTGTACGTGTGGGAGCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGa                                                                                                                                               >  6:218664/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                                                                                                                <  6:319570/149‑1 (MQ=255)
                                    gcagcCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTc                                                                                                                <  2:79575/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                                                                                                                                           >  7:88606/1‑119 (MQ=255)
                                       gccgccGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGt                                                                                                                                           <  8:88606/119‑1 (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)
                                                                    ttcGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAg‑g                                                                               >  4:353673/1‑149 (MQ=255)
                                                                      cGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑Gc                                                                              >  5:199694/1‑148 (MQ=255)
                                                                       gACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑Gc                                                                              >  5:14723/1‑147 (MQ=255)
                                                                       gACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCa                                                                            <  6:199694/149‑1 (MQ=255)
                                                                           gCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGATGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTAc                                                                        >  8:21804/1‑149 (MQ=255)
                                                                            ccAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑Gc                                                                              <  6:185570/142‑1 (MQ=255)
                                                                            ccAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑Gc                                                                              >  5:185570/1‑142 (MQ=255)
                                                                            ccAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTAcc                                                                       >  4:89308/1‑149 (MQ=255)
                                                                              aGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCtt                                                                     >  4:317461/1‑149 (MQ=255)
                                                                                gCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTc                                                                    >  2:88938/1‑148 (MQ=255)
                                                                                 cTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGTCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGa                                                                  <  1:219083/149‑1 (MQ=255)
                                                                                  tttCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGTGGCAGAACCTTGGa                                                                  >  1:127333/1‑149 (MQ=255)
                                                                                      ggCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGc                                                             <  3:245073/149‑1 (MQ=255)
                                                                                        cAATCTAGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGc                                                                                                       <  8:293934/106‑1 (MQ=255)
                                                                                        cAATCTAGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGc                                                                                                       >  7:293934/1‑106 (MQ=255)
                                                                                          aTCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGAt                                                         <  3:353673/149‑1 (MQ=255)
                                                                                               gCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGTAG‑GCCAGTACCTTCGACTTGCGGATGGc                                                      <  1:207384/147‑1 (MQ=255)
                                                                                               gCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATg                                                        <  1:88938/145‑1 (MQ=255)
                                                                                               gCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGc                                                      <  5:218664/147‑1 (MQ=255)
                                                                                               gCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCt                                                    >  6:318628/1‑149 (MQ=255)
                                                                                                 tctGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTgg                                                  >  2:352988/1‑149 (MQ=255)
                                                                                                 tctGCCACAGGGCAGGGCCAGTGATGTGTACCGATTAACAGTTGCTGTCCAAGGCGACGGTCACCGGAATGTACATGACTTAGAACTCGAAGAAAGCATCCAAGCCCAGAACGGAGAAC‑GCCAGAAAAGTAGACTTGCGGaagtccac                                                   >  8:348856/1‑141 (MQ=255)
                                                                                                  ctGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTggg                                                 >  7:41233/1‑149 (MQ=255)
                                                                                                   tGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCTGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTgg                                                  <  8:35643/147‑1 (MQ=255)
                                                                                                     ccACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTggg                                                 <  8:41233/146‑1 (MQ=255)
                                                                                                     ccACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTGGGc                                                <  1:179637/147‑1 (MQ=255)
                                                                                                     ccACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTGGGc                                                <  3:89308/147‑1 (MQ=255)
                                                                                                      cacaGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTGGGCGAc                                             >  3:358048/1‑149 (MQ=255)
                                                                                                      cacaGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTGGGCGAc                                             >  7:84497/1‑149 (MQ=255)
                                                                                                       acaGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTACTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTGGGCGAc                                             >  5:10410/1‑148 (MQ=255)
                                                                                                        cagggcagggCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTGGGCGAc                                             >  8:162516/1‑147 (MQ=255)
                                                                                                         agggcagggcCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAg                                          <  7:21804/149‑1 (MQ=255)
                                                                                                          gggcagggcCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAgg                                         <  8:49382/149‑1 (MQ=255)
                                                                                                             cagggcCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAg                                      <  5:318628/149‑1 (MQ=255)
                                                                                                                     gTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGcgcc                              <  6:330851/149‑1 (MQ=255)
                                                                                                                       gATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGcgccgc                            >  4:320048/1‑149 (MQ=255)
                                                                                                                            gtACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACg                       >  7:176919/1‑149 (MQ=255)
                                                                                                                               ccGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGcc                    <  2:203196/149‑1 (MQ=255)
                                                                                                                                 gATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCAt                  >  5:15603/1‑149 (MQ=255)
                                                                                                                                   ttCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCa                <  3:353771/149‑1 (MQ=255)
                                                                                                                                     cACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAg               <  3:320048/148‑1 (MQ=255)
                                                                                                                                      aCCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAgg              >  6:237264/1‑148 (MQ=255)
                                                                                                                                      aCCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAgg              >  3:65428/1‑148 (MQ=255)
                                                                                                                                              tGTCAACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAGGTACACGGcc     >  2:85070/1‑149 (MQ=255)
                                                                                                                                                  cACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAGGTACACGGCCttg  <  5:12208/148‑1 (MQ=255)
                                                                                                                                                  cACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAGGTACACGGCCttg  <  8:84497/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TTGCGCTGTACGTGTGGGAGCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAGGTACACGGCCTTG  >  NC_002947/1037812‑1038105

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: ‑