Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 789,580 C→T E134K (GAA→AAA)  galK ← galactokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913789,5800CT98.9% 323.0 / NA 91E134K (GAA→AAA) galKgalactokinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (50/40);  minor base A (0/1);  total (50/41)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTT  >  NC_000913/789453‑789715
                                                                                                                               |                                                                                                                                       
aGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGga                                                                                                                             >  1:1264668/1‑140 (MQ=255)
      tCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACt                                                                                                                       <  2:928604/140‑1 (MQ=255)
      tCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGAAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACt                                                                                                                       <  1:649865/140‑1 (MQ=255)
          tGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAAc                                                                                                                   >  2:774856/1‑140 (MQ=255)
               ccGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGc                                                                                                              >  1:272392/1‑140 (MQ=255)
               ccGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGc                                                                                                              >  2:1425434/1‑140 (MQ=255)
               ccGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGc                                                                                                              >  2:1253178/1‑140 (MQ=255)
                cGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCg                                                                                                                >  1:84633/1‑137 (MQ=255)
                cGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCg                                                                                                                <  2:84633/137‑1 (MQ=255)
                cGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCa                                                                                                             >  2:1353357/1‑140 (MQ=255)
                cGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCa                                                                                                             <  1:998094/140‑1 (MQ=255)
                cGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCa                                                                                                             >  1:847278/1‑140 (MQ=255)
                cGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCa                                                                                                             <  1:1253178/140‑1 (MQ=255)
                cGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCa                                                                                                             >  2:787226/1‑140 (MQ=255)
                 gCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCAc                                                                                                            >  1:52865/1‑140 (MQ=255)
                 gCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCAc                                                                                                            >  1:985996/1‑140 (MQ=255)
                     ttACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTg                                                                                                        >  2:522285/1‑140 (MQ=255)
                       aCAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCg                                                                                                      >  2:512727/1‑140 (MQ=255)
                       aCAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCg                                                                                                      >  2:1207339/1‑140 (MQ=255)
                         aGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGc                                                                                                    >  2:573343/1‑140 (MQ=255)
                            cTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGcaca                                                                                                 >  2:582031/1‑140 (MQ=255)
                            cTACAAACTGGTTTTCTGCTTCCTGACCGCTAAGCGCGTTTTGTGCGCCGTCCAGCGGCAGTTGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGtaca                                                                                                 >  2:346407/1‑140 (MQ=255)
                             tACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACAt                                                                                                >  2:1172153/1‑140 (MQ=255)
                               cAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTg                                                                                              >  1:873356/1‑140 (MQ=255)
                               cAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTg                                                                                              >  2:1016111/1‑140 (MQ=255)
                               cAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTg                                                                                              >  1:888576/1‑140 (MQ=255)
                                aaaCTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGc                                                                                             <  2:873356/140‑1 (MQ=255)
                                 aaCTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGcc                                                                                            <  1:512727/140‑1 (MQ=255)
                                 aaCTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGcc                                                                                            <  1:774856/140‑1 (MQ=255)
                                 aaCTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGcc                                                                                            <  2:928717/140‑1 (MQ=255)
                                  aCTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCg                                                                                           >  2:1453921/1‑140 (MQ=255)
                                    tGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCt                                                                                         <  2:882832/140‑1 (MQ=255)
                                     ggTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTg                                                                                        <  1:522285/140‑1 (MQ=255)
                                     ggTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTg                                                                                        <  1:1172153/140‑1 (MQ=255)
                                      gTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGa                                                                                       <  1:573343/140‑1 (MQ=255)
                                      gTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGa                                                                                       <  1:1233656/140‑1 (MQ=255)
                                      gTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGa                                                                                       <  1:1214266/140‑1 (MQ=255)
                                      gTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGa                                                                                       >  1:977183/1‑140 (MQ=255)
                                        tttCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATc                                                                                     <  2:242309/140‑1 (MQ=255)
                                          tCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCAc                                                                                   <  2:921309/140‑1 (MQ=255)
                                              cTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATg                                                                               <  1:1353357/140‑1 (MQ=255)
                                                     aCCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCAcgc                                                                        >  1:1010871/1‑140 (MQ=255)
                                                        gTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCAcgccgc                                                                     >  1:564757/1‑140 (MQ=255)
                                                        gTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCAACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCAcgccgc                                                                     >  1:266178/1‑140 (MQ=255)
                                                         ttAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACgccgcc                                                                    >  2:386738/1‑140 (MQ=255)
                                                          tAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGccgccg                                                                   >  2:287687/1‑140 (MQ=255)
                                                          tAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGccgccg                                                                   <  1:1425434/140‑1 (MQ=255)
                                                           aaGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGa                                                                  >  2:1296909/1‑140 (MQ=255)
                                                              cgcgATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGc                                                               <  1:46235/140‑1 (MQ=255)
                                                              cgcgATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGc                                                               <  1:287687/140‑1 (MQ=255)
                                                              cgcgATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGc                                                               <  2:1264668/140‑1 (MQ=255)
                                                                 gATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCtgt                                                            >  1:536075/1‑140 (MQ=255)
                                                                 gATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCtgt                                                            >  2:767744/1‑140 (MQ=255)
                                                                   tttGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTGACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCtgttg                                                          >  2:1141305/1‑140 (MQ=255)
                                                                   tttGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCtgttg                                                          >  1:301758/1‑140 (MQ=255)
                                                                        gcgcCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACg                                                     <  2:301758/140‑1 (MQ=255)
                                                                                      cAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGtt                                       >  1:685593/1‑140 (MQ=255)
                                                                                        gatgatAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAg                                                  <  1:1083242/127‑1 (MQ=255)
                                                                                        gatgatAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAg                                                  >  2:1083242/1‑127 (MQ=255)
                                                                                        gatgatAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTc                                     >  1:111255/1‑140 (MQ=255)
                                                                                        gatgatAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTc                                     >  2:203943/1‑140 (MQ=255)
                                                                                        gatgatAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTc                                     >  1:601220/1‑140 (MQ=255)
                                                                                             tAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCa                                             >  2:726095/1‑127 (MQ=255)
                                                                                             tAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCa                                             <  1:726095/127‑1 (MQ=255)
                                                                                               aaGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCAccacgc                              <  2:536075/140‑1 (MQ=255)
                                                                                                aGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCAccacgcc                             <  1:787226/140‑1 (MQ=255)
                                                                                                 gctgctGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTGACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCAccacgcca                            <  1:1141305/140‑1 (MQ=255)
                                                                                                 gctgctGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCAccacgcca                            <  1:412850/140‑1 (MQ=255)
                                                                                                 gctgctGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCAccacgcca                            <  2:977183/140‑1 (MQ=255)
                                                                                                 gctgctGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCAccacgcca                            <  1:386738/140‑1 (MQ=255)
                                                                                                 gctgctGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCAccacgcca                            <  2:111255/140‑1 (MQ=255)
                                                                                                    gctgcAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGa                         >  2:124719/1‑140 (MQ=255)
                                                                                                     ctgcAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCTACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGaa                        >  1:289280/1‑140 (MQ=255)
                                                                                                     ctgcAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGaa                        >  2:359989/1‑140 (MQ=255)
                                                                                                     ctgcAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGaa                        <  2:888576/140‑1 (MQ=255)
                                                                                                     ctgcAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGaa                        <  2:847278/140‑1 (MQ=255)
                                                                                                     ctgcAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGaa                        >  2:337044/1‑140 (MQ=255)
                                                                                                     ctgcAATACGGTTCCGACCGCGACATTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGaa                        >  1:1126249/1‑140 (MQ=255)
                                                                                                      tgcAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAAc                       <  2:685593/140‑1 (MQ=255)
                                                                                                              ggTTCCGCCCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAg               <  2:272392/140‑1 (MQ=255)
                                                                                                               gTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGc              <  1:1207339/140‑1 (MQ=255)
                                                                                                                ttCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGcc             >  1:1426568/1‑140 (MQ=255)
                                                                                                                ttCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGcc             >  2:1093661/1‑140 (MQ=255)
                                                                                                                  ccGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCa           >  1:585614/1‑140 (MQ=255)
                                                                                                                    gACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCAtt         <  2:1010871/140‑1 (MQ=255)
                                                                                                                    gACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCAtt         <  1:124719/140‑1 (MQ=255)
                                                                                                                     aCCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTg        >  2:769355/1‑140 (MQ=255)
                                                                                                                     aCCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTg        >  2:32994/1‑140 (MQ=255)
                                                                                                                          gACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGtt                 <  2:1308857/126‑1 (MQ=255)
                                                                                                                          gACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGtt                 >  1:1308857/1‑126 (MQ=255)
                                                                                                                          gACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGt   >  2:470564/1‑140 (MQ=255)
                                                                                                                           aCTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGtt  <  1:1453921/140‑1 (MQ=255)
                                                                                                                               |                                                                                                                                       
AGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTT  >  NC_000913/789453‑789715

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: