Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,596,405 C→T P139L (CCG→CTG)  lgoD → L‑galactonate oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,596,4050CT90.7% 134.3 / 3.4 54P139L (CCG→CTG) lgoDL‑galactonate oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base T (32/17);  total (35/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

TGTTTTAGGCCATGAAATATGTGGGGAGATTGTTGGGCTGGGTAAAAATATTGCTGATCTTAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGGCGTGTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGGGCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGA‑TACCAGTCCGGCGCGCCGTGAACA‑TGTGG‑CAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGATTTTGACGCGCA  >  NC_000913/4596163‑4596651
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                         
tgtTTTAGGCCATGAAATATGTGGGGAGATTGTTGGGCTGGGTAAAAATATTGCTGATCTTAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGGCGTGTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTTAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCTGCAggcgg                                                                                                                                                                                                                                                   >  1:10268/1‑251 (MQ=255)
 gtTTTAGGCCATGAAATATGTGGGGAGATTGTTGGGCTGGGTAAAAATATTGCTGATCTTAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGGCGTGTAAAAGAGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGGGGCGAACCTTTTGCCCGCAGAAGGTATTGACCTGCAggcggc                                                                                                                                                                                                                                                  >  2:589245/1‑251 (MQ=255)
              aaaTATGTGGGGAGATTGTTGGGCTGGGTAAAAATATTGCTGATCTTAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGGCGTGTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTTAGCGTGCCGGTGGCGAACCTTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCGGCattttttga                                                                                                                                                                                                                                      >  2:177862/1‑244 (MQ=255)
              aaaTATGTGGGGAGATTGTTGGGCTGGGTAAAAATATTGCTGATCTTAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGGCGTGTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCAGCAttgattga                                                                                                                                                                                                                                      >  2:176605/1‑250 (MQ=255)
               aaTATGTGGGGAGATTGTTGGGCTGGGTAAAAATATTGCTGATCTTAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGGCGTGTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGGAGCATTGATTGAAc                                                                                                                                                                                                                                    >  2:184133/1‑251 (MQ=255)
                  atGTGGGGAGATTGTTGGGCTGGGTAAAAATATTGCTGATCTTAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGGCGTGTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCAGCCTTGATTGAACCtt                                                                                                                                                                                                                                 >  2:169203/1‑251 (MQ=255)
                         gagaTTGTTGGGCTGGGTAAAAATATTGCTGATCTTAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGGCGTGTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTAt                                                                                                                                                                                                                          <  1:223182/251‑1 (MQ=255)
                           gaTTGTTGGGCTGGGTAAAAATATTGCTGATCTTAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGGCGTGTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTAtt                                                                                                                                                                                                                         >  1:228626/1‑250 (MQ=255)
                                 tgGGCTGGGTAAAAATATTGCTGATCTTAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGGCGTGTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCt                                                                                                                                                                                                                   <  1:20292/250‑1 (MQ=255)
                                        ggTAAAAATATTGCTGATCTTAAAAATGGGCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGGCGTGTAAAAGCGGGCGTACCAATTGCTGTGAAAAAAGTTCAGTCATTGTCGTGCATCAGGATGGCGGTTTGAGTGAGTATTTGAGCGTTCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCAGCATTGAGTGAACCTTTCGCTATTAGCGCTCATGCGGt                                                                                                                                                                                                           <  2:597141/251‑1 (MQ=255)
                                        ggTAAAAATAGTGCTGATCTTAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGGCGTGTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGt                                                                                                                                                                                                           <  2:119904/251‑1 (MQ=255)
                                         gTAAAAATATTGCTGATCTTAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGGCGTGTAAAAGCGGGCGTACCAATTGCTGTGAAAAACTTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGCGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCAGGCGGTg                                                                                                                                                                                                          >  1:388822/1‑251 (MQ=255)
                                         gTAAAAATATTGCTGATCTTAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGGCGTGTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTg                                                                                                                                                                                                          >  2:52977/1‑251 (MQ=255)
                                            aaaaTATTGCTGATCTTAAAACTGGTCAGCAAGTTGCTGTGATCCATTATGTTGCCTGTAAGCAATGCCCGGCGTGTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGGTTTAGGGAGTATTTGAGCGTGCCGGGGGCGAACATTTTGCCCGAAGACGGTATTGACCTGCAGGAGGCAGCATTGATTGTACCTTTCGCTATTAGCGCTCATGCgggtcgg                                                                                                                                                                                                       >  2:431559/1‑246 (MQ=255)
                                                   tGCTGATCTTAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGGCGTGTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCAGCAATGATTGAACCCTTCGCTATTAGCGCTCATGCGGTGCGTCgcgc                                                                                                                                                                                                  >  1:37156/1‑249 (MQ=255)
                                                           ttAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGGCGTGTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCt                                                                                                                                                                                        >  1:249251/1‑251 (MQ=255)
                                                            tAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGGCGTGTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTtctc                                                                                                                                                                                       >  1:412254/1‑251 (MQ=255)
                                                            tAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGGCGTGTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCt                                                                                                                                                                                        <  1:151178/250‑1 (MQ=255)
                                                                aaTGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGGCGTGTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCgg                                                                                                                                                                                   >  1:2136/1‑251 (MQ=255)
                                                                       aGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGGCGTGTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGc                                                                                                                                                                              >  1:283654/1‑249 (MQ=255)
                                                                       aGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGGCGTGTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGc                                                                                                                                                                              >  2:223934/1‑249 (MQ=255)
                                                                       aGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGGCGTGTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGGATTGACCTGCAGGCGGCAGCATTGATTGAAACTTTCGCTATTAGCGCTCATGAGGGGCGTCGCGCAGCCATTGCTCACGGCGAGc                                                                                                                                                                              >  1:144426/1‑249 (MQ=255)
                                                                         cAAGTTGCTGTGATCCCTTATGTGGCCTGTCAGCAATGCCCGGCGTGTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGt                                                                                                                                                                          <  2:116621/251‑1 (MQ=255)
                                                                                     aTCCCTTATGTTGCCTGTCAGCAATGCCCGGCGTGTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTTGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGGGGTcgg                                                                                                                                                              >  2:620787/1‑251 (MQ=255)
                                                                                                tGCCTGTCCGCAATGCCCGGCGTGTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGCACATTTTGCCCGCAGACGGTCTTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCC‑GGCGAGCAGGTGCTGGGGGTCGGGGCGGGGGCaa                                                                                                                                                   >  2:238529/1‑250 (MQ=255)
                                                                                                    tGTCAGCAATGCCCGGCGTGTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGGTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAAAATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTCTTAGCGCTCATGGGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGGGCGGGGGCAAtagg                                                                                                                                               >  2:550612/1‑248 (MQ=255)
                                                                                                           aaTGCCCGGCGTGTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGCACATTTTGCCCGCCGACGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTCTTAGCGCTCATGCGGGGCGTCGCGCAGCCATTGCTCCCGGCGGGCAGGGGCTGGTGGTcggggcgggggcacacggtgtgggg                                                                                                                                        >  2:269419/1‑236 (MQ=255)
                                                                                                                              cGGGCGTACCAATTGCGGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCGTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGGGCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAg                                                                                                                     <  2:2136/251‑1 (MQ=255)
                                                                                                                                                         ttCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGTAGACGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCGTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGGGCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCg                                                                                           <  1:571117/250‑1 (MQ=255)
                                                                                                                                                               cATTGGCGTGCGTCCGTATGGCGTTTTTAGTGAGTGTTTGCGCGTGCCGGTGGCGAACATTTTGCCCGCCGACGTTTTTGCCCTGCAGGCGGCTGCATTGATTGCCCCTTTCGCTATTATCTCTCACGCGTTGCGTCTCGCAGCCTTTGCTCCCGGCGATCAGGTGCTGGTGGTCGGGGCGGGGCCACTCGGTCTGGGCGCGGCGGCACTCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCTGA‑TACCa                                                                                   <  2:144426/251‑1 (MQ=255)
                                                                                                                                                                  tggtgTGCATCAGGATGGGGGTTTGAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCAGCAGTGATTGAACCTGTCGCTATTAGCGCTCATGCGGTGGGGCGCGCAGCGATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGGGCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGA‑TACCAGt                                                                                 <  2:310530/246‑1 (MQ=255)
                                                                                                                                                                  tGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGTGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCCGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGGGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGGGGTCGGGGCGGGGGCAAACGGGCTTGGCGCGGGGGGAATCGCTTAAACCGGTGGGGCCCGCGGGGGGGTGGCGGA‑TAcccctc                                                                                >  1:22610/1‑247 (MQ=255)
                                                                                                                                                                                      gTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGGGCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGA‑TACCAGTCCGGCGCGCCGTGAACA‑tgtg                                                           <  1:48454/251‑1 (MQ=255)
                                                                                                                                                                                       ttttAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGGGCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCCCAGGGGGGGGTGGCGGA‑TACCAGTCCGGCGCGCCggggaaatgggg                                                           >  1:22290/1‑240 (MQ=255)
                                                                                                                                                                                         ttAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGGGCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCCCAGGGGGGGGTGGCGGG‑TACCCGTCCGGCGCGCCGTGAACC‑TGTGG‑c                                                        >  1:528473/1‑250 (MQ=255)
                                                                                                                                                                                           aGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGGGCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGCTGGCGCACAGGGGGGGGTGGCGGA‑TACCAGTACGGCGCGCCCGGAACC‑TGTGG‑CAAc                                                     >  1:317777/1‑251 (MQ=255)
                                                                                                                                                                                                         cGTGCCGGTGGCGAACCTTTTGCCCGCAGACGGTATTGCCCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGGGCGGGGCCAAATGGTCTGGGCCCGGCGGCAATCGCTAAAACCGATGGGGCCAAGGTGGGTGTGGCGGA‑TACCCGGCCGGGCCGCCGTGA‑CACTATGG‑CAACGCGTCTGGGATTAc                                       >  2:553479/1‑251 (MQ=255)
                                                                                                                                                                                                          gTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGGGCGGGGCCAATCGGTCTGGGCGCGGCGGCAATAGCTAAATCCGCTGGCGCACAGGTGGTGGTGGCGGA‑TACCAGTCCGGCGCGCCGTGAACA‑TGTGG‑CAACGCGTCTGGGATTAcc                                      >  1:6099/1‑251 (MQ=255)
                                                                                                                                                                                                            gCCGGTGGCTAACATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGGGCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGA‑TACCAGTCCGGCGCGCCGTGAACA‑TGTGG‑CAACGCGTCTGGAATTACCtt                                    <  1:512079/251‑1 (MQ=255)
                                                                                                                                                                                                               ggTGGCGAACATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGGGCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGGGGTGGTGGCGGA‑TACCAGTCCGGCGCGCCGTGAACA‑TGGGG‑CAACGCGGCTGGAATTACCTTTAc                                 >  1:33181/1‑251 (MQ=255)
                                                                                                                                                                                                                 tGGCGAACATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGGGCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGCTGGCGCACAGGTGGTGGGGGCGGG‑TACCCGTCCGGCGCGCCGGGGACA‑TGTGG‑CCAAGCGTCTGGAATTAACTTTACTg                               >  1:386649/1‑251 (MQ=255)
                                                                                                                                                                                                                     gAACATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTCTTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGGGCGGGGCCAATCGGTACGGGCGCGGCGGCAAACGCTAAAACCGAAGGCGCACGGGGGGTGGGGGCGGA‑TATCACTCCGGCGCGCCTGGTACA‑AGTGGACAA‑GCGACTGGGATTACCTTTACTgggac                           >  2:162955/1‑248 (MQ=255)
                                                                                                                                                                                                                      aaCATTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGGGCGGGGCCAATCGGGCTGGGCGCGGCGGCAATCGCTAAAGCCGCTGGCGCACAGGTGGTGGGGGCGGA‑TACCAGTCCGGCGCGCCGGGAACA‑TGTGG‑CAACGCGGCTgggaattccttttcttggacc                          >  1:170537/1‑232 (MQ=255)
                                                                                                                                                                                                                         aTTTTGCCCGCAGACGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGGGCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCCCTAAAGCCGGTGGCGCACAGGTGGGGGTGGCGGG‑CACCAGGCCCGCGCGCCGGGAAAA‑TGTGG‑CAACGCGGCTGGAATTACCTTTTCTTGACCCggc                       >  1:437627/1‑249 (MQ=255)
                                                                                                                                                                                                                            ttGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGGGCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGA‑TACCAGTCCGGCGCGCCGTGAACA‑TGTGG‑CAACGCGTCTGGAATTACCTTTACTGGACCCGTCAg                     <  1:177862/250‑1 (MQ=255)
                                                                                                                                                                                                                                ccGCAGACGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGGGCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGA‑TACCAGTCCGGCGCGCCGTGAACA‑TGTGG‑CAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGaa                <  1:184133/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cGCAGACGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGGGCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGA‑TACCAGTCCGGCGCGCCGTGAACA‑TGTGG‑CAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAg               <  1:549041/251‑1 (MQ=255)
                                                                                                                                                                                                                                    aGACGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGGGCGGGGGCCCTCGGTCTGGGCGCGGGGGCCATCGCCAAAGCCCGTGGCGCCAAGGTGGGGGGGGCGGGGTAC‑ACTCCGGCCCGCCGTGGACG‑TGGGGGCCA‑GGGGCTGGAATTACCTTTTCTGGACCCGGCACCCGGAGAtt            >  2:332152/1‑251 (MQ=255)
                                                                                                                                                                                                                                      aCGGTATTGCCCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGGGCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGA‑TACCAGTCCGGCGCGCCGTGAACA‑TGTGG‑CAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGAtttt          <  1:433035/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aCGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGGGCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGA‑TACCAGTCCGGCGCGCCGTGAACA‑TGTGG‑CAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGAtttt          <  2:530748/251‑1 (MQ=255)
                                                                                                                                                                                                                                       cGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGGGCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGA‑TACCAGTCCGGCGCGCCGTGAACA‑TGTGG‑CAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGcc                   <  2:111719/241‑1 (MQ=255)
                                                                                                                                                                                                                                       cGGTATTGACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGGGCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGGGGGGGTGGCGGA‑TACCAGTCCGGCGCGCCGGGAACA‑TGGGG‑CAACGCGTCTGGAATTACCTTTACTGGACCCGGCAGcc                   >  1:111719/1‑241 (MQ=255)
                                                                                                                                                                                                                                           aTTGACCTGCAGGCGGCAGCATGGATTGAATCTTTCGCTAGTAGCGCTCAGGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGGGCAGGTGCTGGTGGTCGGGGCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGA‑TACCAGTCCGGCGCGCCGTGAACA‑TGTGG‑CAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGATTTTGAcg      <  2:400709/250‑1 (MQ=255)
                                                                                                                                                                                                                                              gACCTGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGGGCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGA‑TACCAGTCCGGCGCGCCGGGAACA‑TGGGG‑CAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGATTTTGACgcgcc  >  1:223655/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                         
TGTTTTAGGCCATGAAATATGTGGGGAGATTGTTGGGCTGGGTAAAAATATTGCTGATCTTAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGGCGTGTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGGGCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGA‑TACCAGTCCGGCGCGCCGTGAACA‑TGTGG‑CAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGATTTTGACGCGCA  >  NC_000913/4596163‑4596651

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