Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,037,957 C→G 6.8% D480H (GAC→CAC)  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,9570CG 6.8% 237.1 / 7.2 74D480H (GAC→CAC) PP_0897class 1 fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (40/29);  new base G (2/3);  total (42/32)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

GCTTGCGCTGTACGTGTGGGAGCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAGGTACACGGCC  >  NC_002947/1037810‑1038102
                                                                                                                                                   |                                                                                                                                                  
gCTTGGGCTGTAGGTGTGGGCGGGGGCTTGCCCGCGAAGCAGCCGGCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTc                                                                                                                                                    <  5:64637/148‑1 (MQ=255)
  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                                                                                                                                         <  5:128738/121‑1 (MQ=255)
                                      gcagcCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTCGCCGTGGACAGCAACgg                                                                                                                                         >  6:128738/1‑121 (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)
                                                                      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                                                                           >  5:185570/1‑142 (MQ=255)
                                                                              ccAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑Gc                                                                           <  6:185570/142‑1 (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           >  3:65428/1‑148 (MQ=255)
                                                                                                                                        aCCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAgg           >  6:237264/1‑148 (MQ=255)
                                                                                                                                                tGTCAACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAGGTACACGGcc  >  2:85070/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCTTGCGCTGTACGTGTGGGAGCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAG‑GCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAGGTACACGGCC  >  NC_002947/1037810‑1038102

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