Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 498,344 C→T A97V (GCG→GTG)  hemH → ferrochelatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913498,3440CT100.0% 261.0 / NA 73A97V (GCG→GTG) hemHferrochelatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (31/42);  total (31/42)

CGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGG  >  NC_000913/498207‑498475
                                                                                                                                         |                                                                                                                                   
cgTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAgtgc                                                                                                                                   >  2:951319/1‑140 (MQ=255)
cgTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAgtgc                                                                                                                                   >  2:1002934/1‑140 (MQ=255)
cgTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAgtgc                                                                                                                                   >  1:856961/1‑140 (MQ=255)
 gTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCt                                                                                                                                  <  2:596290/140‑1 (MQ=255)
 gTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCt                                                                                                                                  <  1:1128998/140‑1 (MQ=255)
      tttGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGaa                                                                                                                             <  1:722229/140‑1 (MQ=255)
          ccGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAg                                                                                                                         >  2:379099/1‑140 (MQ=255)
            gCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCt                                                                                                                       >  2:1184946/1‑140 (MQ=255)
                 gcTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGc                                                                                                                  <  2:609400/140‑1 (MQ=255)
                     gccgcGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGc                                                                                                              <  1:1002934/140‑1 (MQ=255)
                     gccgcGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGc                                                                                                              <  1:379099/140‑1 (MQ=255)
                       cgcgTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCa                                                                                                            <  1:287848/140‑1 (MQ=255)
                          gtgtGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCa                                                                                                         <  2:890858/140‑1 (MQ=255)
                                       gTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCg                                                                                            <  1:951319/140‑1 (MQ=255)
                                       gTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCg                                                                                            >  2:878551/1‑140 (MQ=255)
                                            ccTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGAt                                                                                       <  1:15724/140‑1 (MQ=255)
                                            ccTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGAt                                                                                       <  1:752757/140‑1 (MQ=255)
                                               ctgtctgGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGaa                                                                                    >  1:247437/1‑140 (MQ=255)
                                               ctgtctgGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGaa                                                                                    >  2:193117/1‑140 (MQ=255)
                                                 gtctgGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACt                                                                                  >  1:880566/1‑140 (MQ=255)
                                                  tctgGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAg                                                                                                 >  2:324665/1‑124 (MQ=255)
                                                  tctgGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAg                                                                                                 <  1:324665/124‑1 (MQ=255)
                                                  tctgGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTc                                                                                 >  2:1117762/1‑140 (MQ=255)
                                                      gatggaAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTgg                                                                             <  1:193117/140‑1 (MQ=255)
                                                       atggaAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGc                                                                            >  1:1195042/1‑140 (MQ=255)
                                                       atggaAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGc                                                                            <  2:283649/140‑1 (MQ=255)
                                                       atggaAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGc                                                                            >  1:773883/1‑140 (MQ=255)
                                                        tggaAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCa                                                                           <  1:213203/140‑1 (MQ=255)
                                                        tggaAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCa                                                                           <  2:856961/140‑1 (MQ=255)
                                                         ggaAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCag                                                                          <  1:878551/140‑1 (MQ=255)
                                                         ggaAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCag                                                                          <  1:329166/140‑1 (MQ=255)
                                                          gaAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCaga                                                                         <  1:140786/140‑1 (MQ=255)
                                                          gaAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCaga                                                                         <  2:773883/140‑1 (MQ=255)
                                                            aGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGc                                                                            <  1:657344/135‑1 (MQ=255)
                                                            aGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGc                                                                            >  2:657344/1‑135 (MQ=255)
                                                            aGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCag                                                                          <  2:709052/137‑1 (MQ=255)
                                                            aGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGc                                                                       >  2:911410/1‑140 (MQ=255)
                                                            aGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCATAGATGAACTCCTGGCag                                                                          >  1:709052/1‑137 (MQ=255)
                                                              gTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCAt                                                                     <  1:609317/140‑1 (MQ=255)
                                                                ggCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGt                                                                   <  2:1195042/140‑1 (MQ=255)
                                                                ggCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGt                                                                   >  1:430788/1‑140 (MQ=255)
                                                                ggCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGt                                                                   <  1:1184946/140‑1 (MQ=255)
                                                                 gCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTa                                                                  >  2:920908/1‑140 (MQ=255)
                                                                 gCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTa                                                                  >  1:356132/1‑140 (MQ=255)
                                                                 gCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTa                                                                  >  1:1089389/1‑140 (MQ=255)
                                                                   tCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGa                                                                <  1:1117762/140‑1 (MQ=255)
                                                                   tCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGa                                                                <  1:318972/140‑1 (MQ=255)
                                                                     gccgcTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATc                                                              >  2:736996/1‑140 (MQ=255)
                                                                     gccgcTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATc                                                              <  2:247437/140‑1 (MQ=255)
                                                                      ccgcTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCa                                                             <  1:8366/140‑1 (MQ=255)
                                                                         cTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCa                                                                           <  2:516334/123‑1 (MQ=255)
                                                                         cTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCa                                                                           >  1:516334/1‑123 (MQ=255)
                                                                           gATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATtg                                                        <  2:430788/140‑1 (MQ=255)
                                                                             tGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATtgtg                                                      <  2:356132/140‑1 (MQ=255)
                                                                               gTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTgtggt                                                    <  1:125769/140‑1 (MQ=255)
                                                                                      gccgccAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCg                                             >  1:1181200/1‑140 (MQ=255)
                                                                                       ccgccAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGc                                            >  1:101320/1‑140 (MQ=255)
                                                                                         gccAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCtt                                          >  2:511456/1‑140 (MQ=255)
                                                                                          ccAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCttt                                         <  2:1181200/140‑1 (MQ=255)
                                                                                               aaCAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCg                                    <  2:1089389/140‑1 (MQ=255)
                                                                                                 cagcagGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCa                                  >  1:665743/1‑140 (MQ=255)
                                                                                                 cagcagGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCa                                  >  1:723842/1‑140 (MQ=255)
                                                                                                  agcagGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCaa                                 <  2:101320/140‑1 (MQ=255)
                                                                                                       gcgcTGGCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATtctc                            <  1:911410/140‑1 (MQ=255)
                                                                                                            ggCACAACGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGtt                       >  2:662536/1‑140 (MQ=255)
                                                                                                                   cGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTAcggt                <  2:723842/140‑1 (MQ=255)
                                                                                                                   cGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTAcggt                <  1:920908/140‑1 (MQ=255)
                                                                                                                   cGTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTAcggt                <  2:880566/140‑1 (MQ=255)
                                                                                                                    gTTTACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTAcggtc               <  1:736996/140‑1 (MQ=255)
                                                                                                                       tACCGGAGATGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTAcggtcggt            <  2:412191/140‑1 (MQ=255)
                                                                                                                               aTGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATg    >  2:695760/1‑140 (MQ=255)
                                                                                                                               aTGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATg    >  2:148301/1‑140 (MQ=255)
                                                                                                                               aTGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATg    >  1:23989/1‑140 (MQ=255)
                                                                                                                                tGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATgg   >  1:1057422/1‑140 (MQ=255)
                                                                                                                                tGCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATgg   >  2:857110/1‑140 (MQ=255)
                                                                                                                                 gCCCGTAGTGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATggg  >  2:554773/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                   
CGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGG  >  NC_000913/498207‑498475

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: