Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 789,604 C→T G126R (GGG→AGG)  galK ← galactokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913789,6040CT97.6% 114.8 / NA 41G126R (GGG→AGG) galKgalactokinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (20/20);  minor base G (0/1);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCG  >  NC_000913/789359‑789840
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            
gacAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCTGGCAc                                                                                                                                                                                                                                         <  1:362270/251‑1 (MQ=255)
          ccTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGTACTTAACCTGGCACCCTGCGGcac                                                                                                                                                                                                                               >  2:123124/1‑251 (MQ=255)
           cTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCACGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGTTTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCAGACCGCGACTTCAAGTGAAGCGGAAGAACTTAACATGGCACCCTGCGGcca                                                                                                                                                                                                                               >  2:171062/1‑248 (MQ=255)
              tGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCTGGCACCCTGCGGCACAtt                                                                                                                                                                                                                            <  1:175832/250‑1 (MQ=255)
               gggCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGACGCGGAAGAACTTAACCTGGCACCCTGCGGCACATTGc                                                                                                                                                                                                                          >  2:409784/1‑251 (MQ=255)
                ggCATGGAAACTGCTTTGTGCCCCAGGGTGCGGCAATCGATCAGCAAGGCATGATCTTTCTTTCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAATCTGGTGTTCTGCTACCTGACCGCTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCTGGCACCCTGCGGCACATTGcc                                                                                                                                                                                                                         <  2:214502/251‑1 (MQ=255)
                ggCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTAAACCTGGCACCCTGCGGAACATTGcc                                                                                                                                                                                                                         >  2:7545/1‑251 (MQ=255)
                   aTGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACAGCGACTTCCCGTGACGAGGAAGAACTTAACCTGGCACCCTGCGGCACATTGCCGCt                                                                                                                                                                                                                      >  2:149910/1‑251 (MQ=255)
                        aacgtcTTTTGTCCCCAGTGATCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGGGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGAGTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTACGGCCGCGACTGCCAGTGAAGCGGAAGAACTTAACCTGGCACCCTGCGGCACATTGCCGCTGATCa                                                                                                                                                                                                                 <  2:115194/246‑1 (MQ=255)
                               ttGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCTGGCACCCTGCGGCACATTGCCGCTGATCACCATGt                                                                                                                                                                                                           <  1:15221/250‑1 (MQ=255)
                                    ccccAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCACTACGGTTCCGACCGCGCCTTCCAGTGAAGCGGAAGAACTTAACCTGGCACCCTGCGGCACATTGCCGCTGATCaacaatcccacgc                                                                                                                                                                                                     >  2:158417/1‑239 (MQ=255)
                                    ccccAGTGAGCGGCAATCGAGCAGCAAGGCATGATCTTTCTGGCCGTGCGCGGAAATTAGTTGATCCATGATCCCGCAGTTACAGCCTACAAACGTGTTTTCGGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCTGGCTCCCTGCGGCACATTGCCGATGATCACCATGTCCAcg                                                                                                                                                                                                      <  2:376033/250‑1 (MQ=255)
                                        aGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCTGGCCCCCTGCGGCACATTGCCGCTGATCACCATGTCCACgccgcc                                                                                                                                                                                                 <  1:149910/251‑1 (MQ=255)
                                                      gATCAGCAAGGCATGATCTTTCTTGCCGCACGCGGAAATTAGCTGATCCATGGTCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGAAATACGGTTCCGACCGAGACTTCCAGTGAAGCGGAAGAACTTAACCTGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTAc                                                                                                                                                                                   >  1:124347/1‑251 (MQ=255)
                                                                    gATCGTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCTGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGt                                                                                                                                                                     <  2:97081/251‑1 (MQ=255)
                                                                         ttcttGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCTGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCGGTTGTTACGCAGTTGCAGATGTTTcac                                                                                                                                                                >  2:55009/1‑251 (MQ=255)
                                                                                cGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCTGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCCGCTGTTTCAccacgcc                                                                                                                                                          >  2:286920/1‑250 (MQ=255)
                                                                                cGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCTGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCAccacgccc                                                                                                                                                         >  2:16934/1‑250 (MQ=255)
                                                                                   gcgcGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCTGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGa                                                                                                                                                      <  2:114073/251‑1 (MQ=255)
                                                                                     gcgGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCTGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAAc                                                                                                                                                    <  1:38062/251‑1 (MQ=255)
                                                                                          aaTTAGCTGATCCATGATCCCGCAGTTACAGCCTAAAAACTGGTTTTCTGCGTCCTGACCGTTAAGCGCGTTTTCTGCGCCGTCCAACGGCAGATGATAAAGCTGCTGAAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCTGGCACCCGAAGGCACATGGCCGCTGATCCCCATCTCCACGCCGCCGAAGCTGGTGTTACGCAATTGCAGATATTTCAACACCCCAAGAACGTaat                                                                                                                                               >  1:392549/1‑249 (MQ=255)
                                                                                                   aTCCATGATCCCGCAGTTACAGCCTACAACCTGGTTTTCTGCTTCCTGACCGTTCAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCTGGCACCCTGCGGCACCTTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCCCGCCACGAACGTAGTTAGcccctt                                                                                                                                      >  2:191769/1‑248 (MQ=255)
                                                                                                               gCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGTCCGCGACTTCCAGTGAAGCGGAAGAACTTAACCTGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGGCGTTGAAGATGTTTCACCACGCCACGAACGAAGTTAGCCCATTGATAGTTTTCAt                                                                                                                          >  2:94351/1‑251 (MQ=255)
                                                                                                                                  tGGTTTTCGGCTTCCGGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCATTGAAGCGGAAGAACTTAACCTGGCACCCTGCTGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGa                                                                                                                                                                        <  2:16519/186‑1 (MQ=255)
                                                                                                                                  tGGTTTTATGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCTGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGa                                                                                                                                                                        >  1:16519/1‑186 (MQ=255)
                                                                                                                                  gggTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGGCAAGAGGCAGATGATAAAGCTGCTGCATTACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCTGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATTTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTGTCATGTGCGACAATGGGCGCAt                                                                                                        <  2:251621/249‑1 (MQ=255)
                                                                                                                                   ggTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCTGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCg                                                                                                      >  2:303905/1‑251 (MQ=255)
                                                                                                                                    gTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCTGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCgg                                                                                                     >  1:131756/1‑250 (MQ=255)
                                                                                                                                           gCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCTGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGaaa                                                                                               <  1:286920/250‑1 (MQ=255)
                                                                                                                                                      gTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCTGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCt                                                                                   >  1:306940/1‑251 (MQ=255)
                                                                                                                                                                        gccttCCAGCGGCAGATGAGAAAGCTGCTGCAATCCGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCTGGCACCCTGCGGCACATTGCCGCTGAGCACCATGTCCACGCCGCCGAAGCTTTTGTTACGCAGTTGCAGATGTTTCACCACGCCACTAACGTAGTTAGCCCAATGATAGTGTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTg                                                                 <  2:387773/247‑1 (MQ=255)
                                                                                                                                                                            tCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCTGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGcc                                                               <  1:16934/249‑1 (MQ=255)
                                                                                                                                                                                  ggCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCTGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGCTTTTCATAATCGGCTGCCATCAcgcg                                                       >  1:103332/1‑251 (MQ=255)
                                                                                                                                                                                  ggCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCTGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCAcgcg                                                       <  1:55009/251‑1 (MQ=255)
                                                                                                                                                                                          acaaaGCTGCTGCACTCCGGGACCGACCGCGACTACCAGGGAAGCGGAAGACCTTACCCTTGCATCCTGCGGCACTTTTCCGCTGAACACCAGGTCCACGCCGCCGAAGTGGTTTTTACGCAGTTGATGATGTTTCCCCACGCCACGAATGTAGTTAGCCCATTGATAGTTTGCATGTGCGACAATGGGCGCATCGAGGGAAAACTTGACGAGCTGCTTTTCATAATCGGCTGCCATCAcgtgatttgtac                                               <  2:11366/249‑11 (MQ=255)
                                                                                                                                                                                              aGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCTGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTc                                           <  1:409784/251‑1 (MQ=255)
                                                                                                                                                                                                ctgctgCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCTGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCAt                                         >  1:4109/1‑251 (MQ=255)
                                                                                                                                                                                                        atccGGTTCCGACAGCGACTGCCAGTGAAGCGGAAGATCTTAGCGTTGCTCCCTGCGGCTCATTGCCGCTGAGCCCCCTGTCCACGCCGCCGAAGCTGTTGTTCCGCAGTTGCAGATGTTGCACCACGCCACGAACGGAGTTAGCCCATTGAGAGTTTTCAGGTGCTACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTTATTTTCGTAATCGGCTGACATCACGCGAACTTGACTGTTTTCGCgtggt                                 <  2:193786/248‑1 (MQ=255)
                                                                                                                                                                                                                 cGACCGCGACTTCCAGTGAAGAGGAAGAACTTAACCTGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTg                        >  1:334776/1‑251 (MQ=255)
                                                                                                                                                                                                                                 ggacgCGGAAGAACTTAACCGGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGGCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCTACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCTCGGTTTGATaat        <  2:394531/247‑1 (MQ=255)
                                                                                                                                                                                                                                       ggAAGAACTTAACCTGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATcg  <  1:107625/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            
GACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCG  >  NC_000913/789359‑789840

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑