Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 789,580 C→A E134* (GAA→TAA)  galK ← galactokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913789,5800CA94.4% 187.3 / ‑2.1 71E134* (GAA→TAA) galKgalactokinase
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base A (26/41);  total (29/42)
Fisher's exact test for biased strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GTACGTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACG  >  NC_000913/789330‑789824
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    
gTACGTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTc                                                                                                                                                                                                                                                      >  1:100893/1‑251 (MQ=255)
    gTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGAGACTTACAGt                                                                                                                                                                                                                                                  >  2:36936/1‑251 (MQ=255)
     tttGAAGATACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTg                                                                                                                                                                                                                                                 >  1:242847/1‑251 (MQ=255)
                       gatgaCGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAAc                                                                                                                                                                                                                                <  1:36936/250‑1 (MQ=255)
                       gagtccgacAGCCACCCCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCTGTTACAGCCTGCAAACTGGTTTTCTGCGTCCTGACCGTTAAGCGCGATTTGTGCGCCTTCCAGCGGCAGGTGATAAATCTGCTGCAATACGTTTCTGACCGCTACTTACAGTGAAGCGGAAGAACTTAAcc                                                                                                                                                                                                                               <  2:672769/246‑1 (MQ=255)
                           acgacAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTACTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGc                                                                                                                                                                                                                           >  1:300958/1‑251 (MQ=255)
                             tacAGCCACACCTTTGGGCATGGAACCTGCTGTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGCTTTGTGCGCCGTCCAGCGGCAGCTGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGGGACGCGGAAGAACTTAACCCGGCa                                                                                                                                                                                                                          >  2:148298/2‑250 (MQ=255)
                                    acacCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCgg                                                                                                                                                                                                                  <  1:595932/251‑1 (MQ=255)
                                    acacCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCgg                                                                                                                                                                                                                  <  1:306251/251‑1 (MQ=255)
                                           tGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTGAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACAtt                                                                                                                                                                                                            <  1:29680/250‑1 (MQ=255)
                                            gggCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTg                                                                                                                                                                                                           <  1:316062/250‑1 (MQ=255)
                                             ggcaagttatcTGCTTTGGTCCCCAGTGAGCTGCAATCGATCAGCAAGGCATGATCTGTCTTGCCGAGCGCGTAAATTAGCTGATCCATGATCCCGCCGATACAGGCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGcc                                                                                                                                                                                                         <  2:607113/241‑1 (MQ=255)
                                             ggCATGGAAACTGCTTTGGGCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGcc                                                                                                                                                                                                         <  1:206658/251‑1 (MQ=255)
                                                   gAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGAt                                                                                                                                                                                                   >  1:275938/1‑251 (MQ=255)
                                                   gAAACTGCTTTGGTCCCCAGTCAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAACTTAGCTGATCCATGATACCGCAGTTACAGCCTACAAACTGGATTTCTGCTTCCTGACCGTTAAACGCGATTTGTGCGACGTCCAGCGGCAGTTGCTAAAGCTGCTGCAATACGGTTCCGACGGCGCATTACAGTGACGAGGAAGCACTTAACCCGGCACCCTGCGGCACATTGCCGCTGAt                                                                                                                                                                                                   >  2:239152/1‑251 (MQ=255)
                                                        tGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACCGCCTACCAACTGGCTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGTTAAAGCTGCTGCAATTCGGTCTCGACCGCGACTTACGGTTAAGCGGGAGAACTTAACCCGGCACCCTGAGGCACATTGTCGCTGATCAca                                                                                                                                                                                               >  2:262473/1‑249 (MQ=255)
                                                         gCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTATCCCGGCACCCTGCGGCACATTGCCGCTGATCACCa                                                                                                                                                                                              <  1:507150/250‑1 (MQ=255)
                                                             tGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTc                                                                                                                                                                                          <  1:573318/250‑1 (MQ=255)
                                                                      gTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTTCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCCTGTCCACGccgccg                                                                                                                                                                                >  1:555168/1‑251 (MQ=255)
                                                                                         cAAGGCATGATCTTTCGTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGAAGTTACCACCTACAAACTTGTTTTCTGCTGCCTGACCGTGAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTTCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCTGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTTTTACGCAGtt                                                                                                                                                             <  2:621466/251‑1 (MQ=255)
                                                                                             gctggatatttcctGCCGAGAGCGGAACTTGGCTGGTCCTGTATCCCGTAGTTACAGCCTACAAACTGGTTTGCTGCTTCCTTACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGGGGAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCGCCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCaa                                                                                                                                                         <  2:461399/243‑2 (MQ=255)
                                                                                              cATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCACTACGGTTCCGACCGCGACTTACAGTGAAGCGGCAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGttcccg                                                                                                                                                         >  2:235459/1‑246 (MQ=255)
                                                                                               aTGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCa                                                                                                                                                                                                                          >  1:68016/1‑184 (MQ=255)
                                                                                               aTGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGc                                                                                                                                                           >  1:134260/1‑247 (MQ=255)
                                                                                               aTGATCTTTCTTGCCGAGCGCGGAAATTAGCGGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCa                                                                                                                                                                                                                          <  2:68016/184‑1 (MQ=255)
                                                                                                      ttcttGCTGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCTCCGAAGCTGTTGTTACGCAGTTGCAGATGttt                                                                                                                                                   <  1:245520/248‑1 (MQ=255)
                                                                                                            ccGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCAccacgcc                                                                                                                                          <  1:670262/251‑1 (MQ=255)
                                                                                                             cGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCAccacgccc                                                                                                                                         >  1:62561/1‑250 (MQ=255)
                                                                                                               aGCGCGGAACTTAGCTGAACCATGACCCCGCAGTAACAGCCTACAAACTGGTTTTCTGCTCCCTGGCCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTGCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCTTGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCAccacgccacg                                                                                                                                       <  1:92033/251‑1 (MQ=255)
                                                                                                                    ggAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCTACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGt                                                                                                                                  <  1:100493/251‑1 (MQ=255)
                                                                                                                    ggAAATTAGCCGATCCATGATCCCGCATTTACAGCCTACAAACTGGTTTTCTGCTCCCTGCCCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCGGATGATAAAGCGTCGGCAATACGGTCCCGACCGCTACTTACAGTGAAGCGGAAGAACTTAGCCCGGCCCCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACTAACGt                                                                                                                                  <  2:242847/251‑1 (MQ=255)
                                                                                                                     gAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGt                                                                                                                                  >  1:201908/1‑250 (MQ=255)
                                                                                                                       aaTTAGCTGCTCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGt                                                                                                                               <  1:574926/251‑1 (MQ=255)
                                                                                                                        aTTAGCTGTTCCATGATCCCGCCGTTACAGCCTAAAAACGGGTTTTCTGCTCCCTGACCGTGAAGCGCGATTTGTGCGCCGCCCAGCGGCTCATGATAAAGCTGCTGCATTACGGTTCCGTCCGCGCCTTCCAGTGAAGCGGAAGAACTTACCCCGGCCCCCTGCGGCACATTGCCGCTGATCTCCAGGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACTTAGtt                                                                                                                              <  2:593619/251‑1 (MQ=255)
                                                                                                                                    aTGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGTTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGAAGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATa                                                                                                                  <  1:242210/251‑1 (MQ=255)
                                                                                                                                    aTGAGCCCGCAGGTACAGCCTACAAACGGGTTGTCTGCTTCCTGACCTGTAAGCGCGATTGGTGCGCCTCCCAGCGGCAGAGGATAAAGCTGCTGTAATACGGTTCCGACCGTGACTTACAGTGAAGCGTAAGAACTTAACCAGTCACCCTGCGGCACATTGCCGCTGATCCCCATGTCCACGCCGCCGAAGCTGTTTTTACGCAGTTGCAGATGTTTCACCACGCCACGCACGTAGTTAGCCCATTg                                                                                                                     <  2:297608/248‑1 (MQ=255)
                                                                                                                                              cAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCAtg                                                                                                         >  1:329944/1‑251 (MQ=255)
                                                                                                                                                        cTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATg                                                                                               >  1:195883/1‑251 (MQ=255)
                                                                                                                                                          aCAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATgg                                                                                              <  1:577539/250‑1 (MQ=255)
                                                                                                                                                           cAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGgc                                                                                            >  2:162235/1‑251 (MQ=255)
                                                                                                                                                              aCTGGTTTTCTGCTTACTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCCCGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGCCAATGGgcgc                                                                                          >  1:504297/1‑250 (MQ=255)
                                                                                                                                                                    tttcgacTTCCTGACCGTTAAGTGCGATTTGTGGGCCGCCCAGCGGAAGAAGATAAATCTTCTGCAATACGGTTCCGACCGCGGCTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCAGGTCCACGCCGTCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGGATCGAg                                                                                    <  2:415676/244‑1 (MQ=255)
                                                                                                                                                                    tttCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGCACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCCTGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTCGTTAGCCCATTGATATTTTTCATGTGCGACAATGGGCGCATCGAgg                                                                                   >  1:426557/1‑251 (MQ=255)
                                                                                                                                                                         gCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACAGTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGaaaa                                                                              <  1:371161/251‑1 (MQ=255)
                                                                                                                                                                            tCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGccc                                                                                                                         >  1:16367/1‑205 (MQ=255)
                                                                                                                                                                               tGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCCGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGCGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGACGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGACCGTAGttcgc                                                                                                                           >  2:396408/1‑197 (MQ=255)
                                                                                                                                                                               tGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGc                                                                                                                           <  1:396408/200‑1 (MQ=255)
                                                                                                                                                                                gacctgtaagcccgagttcggggcCGTCCAGCGGCAGATGATAAAGCTGCTGCAATCCGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACAGTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCAGCGAGGGAAAACtcgt                                                                         <  2:329944/227‑1 (MQ=255)
                                                                                                                                                                                            ggtgttGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGAACACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGGAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCa                                                           <  2:581415/247‑1 (MQ=255)
                                                                                                                                                                                             gATTTATGCGCCGTCCAGCGGCAGATGATAAGACTGCTGCAATACGGATCCGACCGCGCCCTACAGTGAAGCGGCAGAACTTAACCCGGCACCCTGCGGCACATTGCCGATGATCACCATGTCCACGCCGCCGAAGCTGTTGTGACGCAGTAGCAGAGGTTTAACCACGCAACTGACAGAGTTATCCCATGGATAGTTTTCATATGCGACAATGGAAGAGTCAAGGGAAAACCAGTCGAGCTGATTTTCAt                                                          >  2:578207/1‑251 (MQ=255)
                                                                                                                                                                                                ttgggcgcCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGGGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCCCGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATaa                                                        <  2:275938/246‑1 (MQ=255)
                                                                                                                                                                                                    gcgcCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCgg                                                    <  2:591896/250‑1 (MQ=255)
                                                                                                                                                                                                      gcCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTg                                                 <  1:269262/251‑1 (MQ=255)
                                                                                                                                                                                                             aGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGACCGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGACGAGCTGATTTTCATCATCGGCTGCCATCAc                                          >  2:43659/1‑251 (MQ=255)
                                                                                                                                                                                                             aGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGtt                                                                                                                                                    <  2:308011/145‑1 (MQ=255)
                                                                                                                                                                                                             aGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGtt                                                                                                                                                    >  1:308011/1‑145 (MQ=255)
                                                                                                                                                                                                              gCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTTAAGCGGAAGAACTGAACCCGGCTCCCGGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTTCGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTGTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCAcg                                         <  1:624531/251‑1 (MQ=255)
                                                                                                                                                                                                              gCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCAcg                                         <  1:490294/251‑1 (MQ=255)
                                                                                                                                                                                                                 gCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGCAGAACTTAACCCGGCACACTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTATTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGAAGTTAGCCCATTGATAGTTTTCATGTGCGAACATGGGAGCATCGAAGGAAAACTCGTCGAGCTGCTTTTCATCATCGGCTGACATCAcgcg                                       >  2:560515/1‑250 (MQ=255)
                                                                                                                                                                                                                   aGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATAGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAAc                                    >  1:308329/1‑251 (MQ=255)
                                                                                                                                                                                                                           aaGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGt                            <  1:60202/251‑1 (MQ=255)
                                                                                                                                                                                                                                   gcAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTATCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGAGTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATcgc                      <  1:595462/249‑1 (MQ=255)
                                                                                                                                                                                                                                      aTACGGTTCCGACCGCGACTTACAGTGAAGCGGACGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGtt                                                                                                                                                             >  2:602364/1‑111 (MQ=255)
                                                                                                                                                                                                                                      aTACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGtt                                                                                                                                                             <  1:602364/111‑1 (MQ=255)
                                                                                                                                                                                                                                       tACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCgtggtg                >  1:483209/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aCGGATCCGACCGCGACTTACAGTGAATCGTAAGAACTTAACACGGCTCCCTGCGGCACATTTCCGCTGAGAACCATGTCCACGCCGCCGAAGCTGTTGTTACTCAGTTGCAGATGTTTCACCACGCCACGAACGCAGTTAGCCCATTGATAGGTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGTCATCACGCGAACTTTACGGTCATCGCgtggtg                <  2:325872/250‑1 (MQ=255)
                                                                                                                                                                                                                                         cGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCTGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGca              <  1:6523/251‑1 (MQ=255)
                                                                                                                                                                                                                                           gtgccgcccGCGATTTACAGTGAAGCTGAAGAAATTAACCCGGCACCCTGCGGCCCATTGCCTGTGATCACCATGTCCCCGCCTGCGAAGTTGTTGTAACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCAGGTGCTACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGcaca            <  2:212638/244‑1 (MQ=255)
                                                                                                                                                                                                                                           gTTCCGACCGCGACTTACAGTGAAGCGGAAGCACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGCGGGAAAACTCATCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGcac             >  2:331309/1‑250 (MQ=255)
                                                                                                                                                                                                                                              ccgcccGCGACTTACAGGGAAGCGGAAAAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCTTGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAAc          <  2:195883/246‑1 (MQ=255)
                                                                                                                                                                                                                                               ccaCCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCCGCACCTTGCCGCTGATCCCCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCCCGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACGTTACGGTCATCGCGTGGTGCACAACTg        <  1:168266/249‑1 (MQ=255)
                                                                                                                                                                                                                                                     cgACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGCACTTTACGGTCATCGCGTGGTGCACAACTGATCACg  >  1:478708/1‑251 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    
GTACGTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACG  >  NC_000913/789330‑789824

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: