Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 789,464 G→A I172I (ATC→ATT galK ← galactokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913789,4640GA100.0% 240.3 / NA 73I172I (ATC→ATTgalKgalactokinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (38/35);  total (38/35)

CGTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACT  >  NC_000913/789333‑789598
                                                                                                                                   |                                                                                                                                      
cGTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGc                                                                                                                                 >  1:849051/1‑139 (MQ=255)
cGTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGc                                                                                                                                 >  2:587060/1‑139 (MQ=255)
      aaGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTAc                                                                                                                           <  2:827554/139‑1 (MQ=255)
      aaGTTACTGTGGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCAGGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTAc                                                                                                                           <  2:1137298/139‑1 (MQ=255)
       aGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACa                                                                                                                          <  1:831941/139‑1 (MQ=255)
            cTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTa                                                                                                                     <  2:777064/139‑1 (MQ=255)
                 gatgatgaCGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAAc                                                                                                                >  1:74077/1‑139 (MQ=255)
                      tgaCGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGtt                                                                                                           >  2:472140/1‑139 (MQ=255)
                       gacgacAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGttt                                                                                                          >  1:336454/1‑139 (MQ=255)
                         cgacAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTc                                                                                                        <  2:849051/139‑1 (MQ=255)
                         cgacAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTc                                                                                                        <  2:390763/139‑1 (MQ=255)
                         cgacAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTc                                                                                                        >  2:73081/1‑139 (MQ=255)
                         cgacAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTc                                                                                                        >  1:448528/1‑139 (MQ=255)
                              gCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCtt                                                                                                   <  2:461863/139‑1 (MQ=255)
                              gCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCACGGAAATTAGCTGATCCATAATCCCGCAGTTAc                                                                                                                           <  1:468632/115‑1 (MQ=255)
                              gCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCACGGAAATTAGCTGATCCATAATCCCGCAGTTAc                                                                                                                           >  2:468632/1‑115 (MQ=255)
                                         gggCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCg                                                                                                                                  <  1:425015/97‑1 (MQ=255)
                                         gggCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCg                                                                                                                                  >  2:425015/1‑97 (MQ=255)
                                         gggCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTa                                                                                        >  1:145111/1‑139 (MQ=255)
                                          ggCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTaa                                                                                       >  2:138295/1‑139 (MQ=255)
                                               ggAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATccc                                                                                                                                   >  2:581133/1‑90 (MQ=255)
                                               ggAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATccc                                                                                                                                   <  1:581133/90‑1 (MQ=255)
                                               ggAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGcgcg                                                                                  >  1:568239/1‑139 (MQ=255)
                                                    cTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCg                                                                                                                                  >  1:980019/1‑86 (MQ=255)
                                                    cTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCg                                                                                                                                  <  2:980019/86‑1 (MQ=255)
                                                     tGCTTTGGTCCCCAGTGAGCGGCAATCGATCAACAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTtgt                                                                            >  1:723217/1‑139 (MQ=255)
                                                        tttGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTgcg                                                                         <  2:448528/139‑1 (MQ=255)
                                                         ttGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTgcgc                                                                        >  2:644911/1‑139 (MQ=255)
                                                           ggTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCg                                                                      >  2:727051/1‑139 (MQ=255)
                                                            gTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTgcgc                                                                        >  2:1175302/1‑136 (MQ=255)
                                                            gTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTgcgc                                                                        <  1:1175302/136‑1 (MQ=255)
                                                            gTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGt                                                                     <  1:693998/139‑1 (MQ=255)
                                                                ccAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAg                                                                 <  2:568239/139‑1 (MQ=255)
                                                                      aGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAg                                                           >  2:394040/1‑139 (MQ=255)
                                                                         ggCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAg                                                                 >  1:482118/1‑130 (MQ=255)
                                                                         ggCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAg                                                                 <  2:482118/130‑1 (MQ=255)
                                                                           cAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAgatgat                                                      <  1:644911/139‑1 (MQ=255)
                                                                            aaTCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTAc                                                                                                                           <  1:467826/69‑1 (MQ=255)
                                                                            aaTCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTAc                                                                                                                           >  2:467826/1‑69 (MQ=255)
                                                                            aaTCGATCAGCAAGGCATAATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATa                                                     >  1:677456/1‑139 (MQ=255)
                                                                                  tCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAgctg                                               <  2:388113/139‑1 (MQ=255)
                                                                                   cAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAgctgc                                              <  1:215768/139‑1 (MQ=255)
                                                                                         ggCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATccc                                                                                                                                   <  1:1133584/48‑1 (MQ=255)
                                                                                         ggCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATccc                                                                                                                                   >  2:1133584/1‑48 (MQ=255)
                                                                                         ggCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAgc                                                 <  2:175838/130‑1 (MQ=255)
                                                                                         ggCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAgc                                                 >  1:175838/1‑130 (MQ=255)
                                                                                          gCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATa                                       <  1:702595/139‑1 (MQ=255)
                                                                                           cATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATAc                                      >  2:665637/1‑139 (MQ=255)
                                                                                                tctttcttGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGctgctg                                            <  1:954631/128‑1 (MQ=255)
                                                                                                tctttcttGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGctgctg                                            >  2:954631/1‑128 (MQ=255)
                                                                                                    tcttGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTgcg                                                                         >  2:493494/1‑95 (MQ=255)
                                                                                                    tcttGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTgcg                                                                         <  1:493494/95‑1 (MQ=255)
                                                                                                    tcttGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGAc                             >  2:123518/1‑139 (MQ=255)
                                                                                                    tcttGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGAc                             >  2:861998/1‑139 (MQ=255)
                                                                                                        gCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACcgcg                         <  2:723217/139‑1 (MQ=255)
                                                                                                        gCCCAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTCGTTTTCTGCTTCCTGACCGTTAAGCACGACTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACcgcg                         <  1:232301/139‑1 (MQ=255)
                                                                                                         ccGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGa                        >  2:112919/1‑139 (MQ=255)
                                                                                                          cGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAg                                                  >  2:511844/1‑112 (MQ=255)
                                                                                                          cGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAg                                                  <  1:511844/112‑1 (MQ=255)
                                                                                                           gAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACt                      >  1:566023/1‑139 (MQ=255)
                                                                                                           gAGCGCGGAAATTAGCTGATCCATAATCCAGCAGTTACAGCCTACAAACTGGTTTTc                                                                                                        <  2:651277/57‑1 (MQ=255)
                                                                                                           gAGCGCGGAAATTAGCTGATCCATAATCCAGCAGTTACAGCCTACAAACTGGTTTTc                                                                                                        >  1:651277/1‑57 (MQ=255)
                                                                                                            aGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACtt                     >  2:669059/1‑139 (MQ=255)
                                                                                                              cgcgGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTcc                   >  2:513349/1‑139 (MQ=255)
                                                                                                              cgcgGAAATTAGCTGATACATAATCCCGCAGTTACAGCCTACAAACTGGTTTTc                                                                                                        >  1:852231/1‑54 (MQ=255)
                                                                                                              cgcgGAAATTAGCTGATACATAATCCCGCAGTTACAGCCTACAAACTGGTTTTc                                                                                                        <  2:852231/54‑1 (MQ=255)
                                                                                                                cgGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAg                 <  1:112919/139‑1 (MQ=255)
                                                                                                                         gCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGga        <  2:336454/139‑1 (MQ=255)
                                                                                                                             aTCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGAt                                                                                >  1:216700/1‑63 (MQ=255)
                                                                                                                             aTCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGAt                                                                                <  2:216700/63‑1 (MQ=255)
                                                                                                                             aTCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGgaagaa    >  2:1079285/1‑139 (MQ=255)
                                                                                                                               ccATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACt  <  1:73081/139‑1 (MQ=255)
                                                                                                                               ccATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACt  <  2:677456/139‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                      
CGTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACT  >  NC_000913/789333‑789598

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