Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,596,481 A→G 22.9% E164G (GAG→GGG)  lgoD → L‑galactonate oxidoreductase, L‑gulonate oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,596,4810AG22.9% 75.9 / 5.7 48E164G (GAG→GGG) lgoDL‑galactonate oxidoreductase, L‑gulonate oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (7/28);  new base G (6/6);  total (14/34)
Fisher's exact test for biased strand distribution p-value = 6.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

ATGTTGCCTGTCAGCAATGCCCGGCGT‑GTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAG‑CAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGA  >  NC_000913/4596256‑4596702
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                               
aTGTTGCCTGTCAGCAATGCCCGGCGT‑GTAAAAGCGGGCGTACCAATTGCTGTTAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAG‑c                                                                                                                                                                                                                              <  1:184314/228‑1 (MQ=255)
  tgtGCCTGACAGCAATGCCCGGCGT‑GTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAG‑CAgg                                                                                                                                                                                                                           <  1:184651/227‑1 (MQ=255)
     gCCTGTCAGCAATGCCCGGCGT‑GTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTTAGCGTGCCGGTGGCGAACATTTTGCCCTCAGACGGTATTGACACGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCGCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCCAG‑AAgggtc                                                                                                                                                                                                                        >  1:512636/1‑226 (MQ=255)
      ccTGTCCGCAATTGCCGGGGT‑TGAAAAGGGGGCCGTCCAATTGGTTGGAAAAAATTTCAGTTATTGGGGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGGGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAG‑CAGGTGCt                                                                                                                                                                                                                       <  2:526184/229‑1 (MQ=255)
      ccTGTCAGCAATGCCCGGCGT‑GTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCCGGCGGCAGCATTGGTTGAACCTTTCGCTATTAGCGCCCATCCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAG‑CAGGTGCt                                                                                                                                                                                                                       >  1:221687/1‑229 (MQ=255)
      ccTGTCAGCAATGCCCGGCGT‑GTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGGT‑CAGGTGCt                                                                                                                                                                                                                       >  1:685658/1‑229 (MQ=255)
      ccTGTCAGCAATGCCCGGCGT‑GTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGGGCCGGTGGCGAACATTTTGCCCGCAGACGGAATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGGCCGTGCGGCAGCCATTGCTCCCGCCGGG‑CAGGGGCt                                                                                                                                                                                                                       >  1:378065/1‑229 (MQ=255)
      ccTGTCAGCAATGCCCGGCGT‑GTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGAGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCCCATTCGGTGCGGCGCGCAGCCATTGCTCCCGCGGTG‑CCGGTGCt                                                                                                                                                                                                                       >  1:339335/1‑229 (MQ=255)
       ctggtaGCAATGCCCGGGGT‑GTAAAAGCGGGCGTACCAATGGCTTTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGGGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAG‑CAGGTGCtg                                                                                                                                                                                                                      <  2:456525/224‑1 (MQ=255)
        tgactctcattcccGGTGT‑TTAAAAAGTGGGGTACCAATTTCTGTGGAAATATTGCAGTCGTTGGGCTTCATCAGGATGGACGGTTTAGTGAGTATTTGAGCGTGCCGGTGTCGAACATTTTGCCCGCAGACGGTATTTACCCGCAGGCGGCAGCATTGATTGAGCCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGGTCTCGGCGAG‑CAGGTGCtgg                                                                                                                                                                                                                     <  2:502/218‑1 (MQ=255)
        tGTCAGCAATGCCCGGCGT‑GTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGGATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGGGCGTCGCGCAGCCATTGCTCCCGGCGAG‑CAGGGGCtgg                                                                                                                                                                                                                     >  1:526184/1‑229 (MQ=255)
            aagaaTGCCCGGGGG‑TGAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCCTTGGCGGGCATCAGGATGGCGGTTTTAGGGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGGCGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAG‑CAGGTGCTggtggt                                                                                                                                                                                                                 <  2:463260/226‑1 (MQ=255)
            aGCAATGCCCGGTGT‑TTAAAAGCGGGAGTACCAATTGCTGGGAAAAAATTTCAGTCATTGGGGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAG‑CAGGTGCTggtggt                                                                                                                                                                                                                 <  2:309373/229‑1 (MQ=255)
             aagaTGGCCGGGGTGGTAAAAGGGGGCGTTCCAATTGCTGTGAAAAAATTTCAGTCATGGGCGTGCCTCAGGATGGCGGTTTTAGTGAGTTTTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAG‑CAGGTGCTggtggt                                                                                                                                                                                                                 <  2:367686/226‑1 (MQ=255)
                 tGCCCGGCGT‑GTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAACGCTCATGCGGTGCGGCGCGCGGCCATTGCTCCCGGCGGG‑CCGGTGCTGGTGGTcggg‑g                                                                                                                                                                                                           >  1:217179/1‑229 (MQ=255)
                 gtcccGGTGT‑TTAAAAGGGGGCGTCCCAAGTGATGGAAAAAAATTTCAGTCATTGGGTTGCATCAGGATGGCGGTTTTAGTGAGTATTTGGGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAG‑CAGGTGCTGGTGGTcggg‑g                                                                                                                                                                                                           <  2:644211/227‑1 (MQ=255)
                    ccGGCGT‑GTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGGATTTGAGCGTGCCGGTGGCGAACCTTTTGCCCCCAGACGGTATTGACCCGCAGGCGGCAGCCTTGGTTGAACCTTTCGCCCTTCGCGCTCATGCGGGGCGTCGCGCGGCCCATGCCCCCGGCCAG‑CCGGTTCTGGGGGTccgg‑ggg                                                                                                                                                                                                         >  2:389385/1‑226 (MQ=255)
                                     ggCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCGGCATTGATTGAACCTTTCGCTATTAGCGCTCAGTCGGTGCGGCGCGCAGCCATTGTCCCCGGCGGG‑CAGGGGCTGGTGGCCGGG‑GGGGGGGCagacgggcgtgg                                                                                                                                                                                        >  2:199256/1‑218 (MQ=255)
                                     ggCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCCGG‑CAGGTGCTGGGGGTCGGG‑GCGGGGGCAAATGGTCCggg                                                                                                                                                                                        >  1:193186/1‑229 (MQ=255)
                                     ggCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGAGCATCAGGAGGGCGGTTTTAGGGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCCCTCATCCGGTGCGCCGGCGAGCCCTTTCTCCCGCCGAC‑CCGGTGGTGGGGGCGGGG‑GGGGGGCCaaaaggcccggg                                                                                                                                                                                        >  2:424944/1‑219 (MQ=255)
                                                tGCTGTGACAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGGGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCCGACGGTATTGACCCGCAGGCGGCGGCCTTGCTTTTACCTTTCGCTATTCGCGCTCATGCGGTGCGTCTCGCAGCCATTACGCCCGGCGGC‑CAGGTGCTGGGGGGCGGGTGGGGGGCCACCGGGCCGGGGGGGGGCGGCa                                                                                                                                                                              >  2:59814/1‑229 (MQ=255)
                                                 gCTGGGAAAAACTTTGAGGCATGGGCGTGGATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGGGCCGGGGGCGAACATTTTGCCCGCAGACGGTATTGGCCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAG‑CAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCaa                                                                                                                                                                             <  1:78787/228‑1 (MQ=255)
                                                              ttgaGTGATGGGCGTGAAGTAGGATGGCGTTTTTAGGGGGTATTTGCGCGTGCCGGTGGCGAACATTTTGCCCGCAGCCGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAG‑CAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGa                                                                                                                                                                <  2:338659/225‑1 (MQ=255)
                                                              ttCAGTCCTGGGCGTGGATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCGGCCATTGCTCCCGGCGAG‑CAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGa                                                                                                                                                                <  1:99515/228‑1 (MQ=255)
                                                                 aGTCATTGGCGTGCCTCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGGGGCGAACATTTTGCCCGCAGACGGGATTGACCCGCAGGCGGCAGCCTTGATTTAAACTTTCGCTATTAGCGCTCATTCGGGGCGTCGCGCACCCCTTTCCCCCGGCGAGCCGGGGGCGGGGGGTCGGC‑GGGGGGCCACCCGCTCGGGGCGcgggggggaccctgaaggcgatggg                                                                                                                                                             >  2:370924/1‑210 (MQ=255)
                                                                                  aGGGTGGGGGTTTTAGGTAGTTTTTGGGCGGGCCGGGGGCGAACATTTTGCCCGGCGGCGGTATTGCCCCGGAGGCGGCCGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAG‑CAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGc                                                                                                                                           <  2:221687/229‑1 (MQ=255)
                                                                                    gATGGCGGTTTTAGTGAGTTTTTGAGCGTGCCGGTGGCGAACATTTGGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAG‑CAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCg                                                                                                                                          <  1:613240/228‑1 (MQ=255)
                                                                                                    ggtttttGGGCGTGCCGGGGGGTAACATTTTGCCCGCGGAAGTGTTTCACCGGTGGGGGGCAGCAGTGATTGGCCCGTTCGGTATGAGCGCTCAGGCGGTGCGGCGCGCAGCCATTGCTCCCGGCGAG‑CAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTACAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGgcgcg                                                                                                                         <  2:343889/225‑1 (MQ=255)
                                                                                                     gttttggggCGTGCCGGTGGCGAACATTTGGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAG‑CAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGgcgcg                                                                                                                         <  1:277571/220‑1 (MQ=255)
                                                                                                        tttGTGCGTGGCGGGTGCGATCTTTCTGCCCGCAGACGGTATTGCCCCGGAGGCGCCAGCATTGATTTGACCTTGCGCTATTAGGGCTCATGCGGGGGGTCCCGCAGCCATTGCTCCCGGCGAG‑CAGGTGCTGGTGGTCGGG‑GCGGGGCCAAACGGTCTGGGTGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGGCGCGCCGt                                                                                                                     <  2:334457/229‑1 (MQ=255)
                                                                                                                     ttGGCGTACTGTTTGCCCGCAGACCGTATTCACCCGGAGGGTGCCGCAGTTACTGAACCCTCCGCTATTACGGCTCATGCGGTGCGTGGGTGAGCCATTGCTCCCGGCGAG‑CAGGTGCTGGGGGTCGGG‑GCGGGGCCAAGCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGGACAGGTGGTGGTGGCGGATACCAGTCCGGCGCTCCGTGAACATGTGGCaa                                                                                                        <  2:263665/228‑1 (MQ=255)
                                                                                                                                  tGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAG‑CAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGGCCGGCGGCCAATGCCAAAGCCGAGTGGGCCAAGGGGGTGCGGGCGGTGACCAGTCCGGCGCGCCCCGGAAACGTGGCaaaccggccggggat                                                                                           >  2:117953/1‑216 (MQ=255)
                                                                                                                                   gCCCGCAGACGGTATTTACCCGCGGGTGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGGGCGTCGCGCAGCCATTGCTCCCGGCGAG‑CAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAAtt                                                                                           <  1:422387/228‑1 (MQ=255)
                                                                                                                                     ccGCAGGCGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCAGGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAG‑CAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAATTAc                                                                                         <  1:114549/228‑1 (MQ=255)
                                                                                                                                                          gggcgggAGCTTTGATTGATTCGTTCGCTGTTGGCGCGCATGGGGGGGGTCGCGCAGCCATTGCTCCCGGCGGG‑CAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCg                                                                    <  2:339335/227‑1 (MQ=255)
                                                                                                                                                                       gattgaACCTTTCGATATTAGCGGGCATGCGGTGGCGCGGTCGGCCAGTGCTCCCGGGGGG‑CAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGATTTTGAcgc                                                       <  1:92887/228‑1 (MQ=255)
                                                                                                                                                                                   cgtttttgggcgtccggCGGGGCGGCGCGCAGCCTTTGCGCCGGGCGGG‑GAGGGGCTGGGGGTCGGG‑GCGGGGCCAATCGGTCTGGGGGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGATTTTGACGCGCAGCTACGGgcg                                          <  2:685657/213‑1 (MQ=255)
                                                                                                                                                                                   cGCTTTTAGCGCTCCTGCGGTGCGTCCCGCCACCATTGCCCCCGGCGAG‑CAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGATTTTGACGCGCAGCTACGGgcg                                          <  2:63167/229‑1 (MQ=255)
                                                                                                                                                                                    gCTATGAGCGCTCATGGGGTGCGGCGCGCAGCTGTTGCTCCGGGCGAG‑CAGGTGCTGGTGGTCGGG‑GCGGGGCCAAGCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGGCGCGACGTGAACATGTGGCAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGATTTTGACGCGCAGCTACGGgcg                                          <  1:41056/228‑1 (MQ=255)
                                                                                                                                                                                       aTTAGTGCTAATGCGGTGCGCAGCGAGGCCAGTGCGCCCGGCGAG‑CAGGTGCTGGTGGTAGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGATTTTGACGCGCAGCTACGGGCGCAg                                       <  1:563951/228‑1 (MQ=255)
                                                                                                                                                                                          gtgcttcttggggcgtggcgcgcgGCCATTGGCCCGGGCGAG‑GAGGTTGGGGTGTGCGGG‑GGGGGGCCAAGCGGTCTGGGGCGGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGATTTTGACGCGCAGCTACGGGCGCAGTTtg                                   <  2:732377/211‑1 (MQ=255)
                                                                                                                                                                                                    gggTGCGTCGGGCAGCCATTGCTCCCGGGGGG‑CAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCTGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCt                         <  1:123077/228‑1 (MQ=255)
                                                                                                                                                                                                    cgggtcggcgcgggtcaaTTGGTCCGGGGGAG‑CGGGGGTTGTGGGTCGGG‑GGGGGGCCAATCGGTCTGGGCGCGGGGGCAATCGGTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCt                         <  2:275931/212‑1 (MQ=255)
                                                                                                                                                                                                       tGCGTCGGGGCGCCATTGCCCCCGGCGAG‑GAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGGTACCAGTCCGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCTgg                       <  1:204624/228‑1 (MQ=255)
                                                                                                                                                                                                                        gcccccGGCGAG‑CAGGCTCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAATTACCTTTACTg                                                                                 <  1:621343/150‑1 (MQ=255)
                                                                                                                                                                                                                         cgcccGGCGAG‑CGGGTGCTGGGGGTGGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGAcgc    <  1:424944/227‑1 (MQ=255)
                                                                                                                                                                                                                         cTCCCGGCGGG‑CAGGTGCTGGTGGTCGGG‑GCGGGGCCAAGCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTTGTGGTGGCGGATACCAGTCCGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGAcg     <  1:199256/228‑1 (MQ=255)
                                                                                                                                                                                                                           cccGGCGGG‑GGGGTGGTGGTGGTGGGG‑GGGGGGCAAATGGGGCTGGGCGCGGCGGCAATCGCTAAAGGCGATGGCGCACGGGTGGTGGTGGCGGATACCAGTCCGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGa  <  2:277318/229‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                               
ATGTTGCCTGTCAGCAATGCCCGGCGT‑GTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAG‑CAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGA  >  NC_000913/4596256‑4596702

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