Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,596,480 A→G 23.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,4800AG23.9% 67.9 / 5.7 43E164G (GAG→GGG) lgoDL‑galactonate oxidoreductase; L‑gulonate oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (5/26);  new base G (6/5);  total (12/31)
Fisher's exact test for biased strand distribution p-value = 2.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGTTGCCTGTCAGCAATGCCCGGCGT‑GTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGA  >  NC_000913/4596255‑4596701
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                              
aTGTTGCCTGTCAGCAATGCCCGGCGT‑GTAAAAGCGGGCGTACCAATTGCTGTTAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGc                                                                                                                                                                                                                              <  1:184314/228‑1 (MQ=255)
  tgtGCCTGACAGCAATGCCCGGCGT‑GTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAgg                                                                                                                                                                                                                           <  1:184651/227‑1 (MQ=255)
     gCCTGTCAGCAATGCCCGGCGT‑GTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTTAGCGTGCCGGTGGCGAACATTTTGCCCTCAGACGGTATTGACACGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCGCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCCAGAAgggtc                                                                                                                                                                                                                        >  1:512636/1‑226 (MQ=255)
      ccTGTCCGCAATTGCCGGGGT‑TGAAAAGGGGGCCGTCCAATTGGTTGGAAAAAATTTCAGTTATTGGGGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGGGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCt                                                                                                                                                                                                                       <  2:526184/229‑1 (MQ=255)
      ccTGTCAGCAATGCCCGGCGT‑GTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCCGGCGGCAGCATTGGTTGAACCTTTCGCTATTAGCGCCCATCCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCt                                                                                                                                                                                                                       >  1:221687/1‑229 (MQ=255)
      ccTGTCAGCAATGCCCGGCGT‑GTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGGTCAGGTGCt                                                                                                                                                                                                                       >  1:685658/1‑229 (MQ=255)
      ccTGTCAGCAATGCCCGGCGT‑GTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGGGCCGGTGGCGAACATTTTGCCCGCAGACGGAATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGGCCGTGCGGCAGCCATTGCTCCCGCCGGGCAGGGGCt                                                                                                                                                                                                                       >  1:378065/1‑229 (MQ=255)
      ccTGTCAGCAATGCCCGGCGT‑GTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGAGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCCCATTCGGTGCGGCGCGCAGCCATTGCTCCCGCGGTGCCGGTGCt                                                                                                                                                                                                                       >  1:339335/1‑229 (MQ=255)
       ctggtaGCAATGCCCGGGGT‑GTAAAAGCGGGCGTACCAATGGCTTTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGGGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCtg                                                                                                                                                                                                                      <  2:456525/224‑1 (MQ=255)
        tgactctcattcccGGTGT‑TTAAAAAGTGGGGTACCAATTTCTGTGGAAATATTGCAGTCGTTGGGCTTCATCAGGATGGACGGTTTAGTGAGTATTTGAGCGTGCCGGTGTCGAACATTTTGCCCGCAGACGGTATTTACCCGCAGGCGGCAGCATTGATTGAGCCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGGTCTCGGCGAGCAGGTGCtgg                                                                                                                                                                                                                     <  2:502/218‑1 (MQ=255)
        tGTCAGCAATGCCCGGCGT‑GTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGGATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGGGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGGGCtgg                                                                                                                                                                                                                     >  1:526184/1‑229 (MQ=255)
            aagaaTGCCCGGGGG‑TGAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCCTTGGCGGGCATCAGGATGGCGGTTTTAGGGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGGCGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTggtggt                                                                                                                                                                                                                 <  2:463260/226‑1 (MQ=255)
            aGCAATGCCCGGTGT‑TTAAAAGCGGGAGTACCAATTGCTGGGAAAAAATTTCAGTCATTGGGGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTggtggt                                                                                                                                                                                                                 <  2:309373/229‑1 (MQ=255)
             aagaTGGCCGGGGTGGTAAAAGGGGGCGTTCCAATTGCTGTGAAAAAATTTCAGTCATGGGCGTGCCTCAGGATGGCGGTTTTAGTGAGTTTTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTggtggt                                                                                                                                                                                                                 <  2:367686/226‑1 (MQ=255)
                 tGCCCGGCGT‑GTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAACGCTCATGCGGTGCGGCGCGCGGCCATTGCTCCCGGCGGGCCGGTGCTGGTGGTcggg‑g                                                                                                                                                                                                           >  1:217179/1‑229 (MQ=255)
                 gtcccGGTGT‑TTAAAAGGGGGCGTCCCAAGTGATGGAAAAAAATTTCAGTCATTGGGTTGCATCAGGATGGCGGTTTTAGTGAGTATTTGGGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTcggg‑g                                                                                                                                                                                                           <  2:644211/227‑1 (MQ=255)
                    ccGGCGT‑GTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGGATTTGAGCGTGCCGGTGGCGAACCTTTTGCCCCCAGACGGTATTGACCCGCAGGCGGCAGCCTTGGTTGAACCTTTCGCCCTTCGCGCTCATGCGGGGCGTCGCGCGGCCCATGCCCCCGGCCAGCCGGTTCTGGGGGTccgg‑ggg                                                                                                                                                                                                         >  2:389385/1‑226 (MQ=255)
                                     ggCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCGGCATTGATTGAACCTTTCGCTATTAGCGCTCAGTCGGTGCGGCGCGCAGCCATTGTCCCCGGCGGGCAGGGGCTGGTGGCCGGG‑GGGGGGGCagacgggcgtgg                                                                                                                                                                                        >  2:199256/1‑218 (MQ=255)
                                     ggCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCCGGCAGGTGCTGGGGGTCGGG‑GCGGGGGCAAATGGTCCggg                                                                                                                                                                                        >  1:193186/1‑229 (MQ=255)
                                     ggCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGAGCATCAGGAGGGCGGTTTTAGGGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCCCTCATCCGGTGCGCCGGCGAGCCCTTTCTCCCGCCGACCCGGTGGTGGGGGCGGGG‑GGGGGGCCaaaaggcccggg                                                                                                                                                                                        >  2:424944/1‑219 (MQ=255)
                                                tGCTGTGACAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGGGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCCGACGGTATTGACCCGCAGGCGGCGGCCTTGCTTTTACCTTTCGCTATTCGCGCTCATGCGGTGCGTCTCGCAGCCATTACGCCCGGCGGCCAGGTGCTGGGGGGCGGGTGGGGGGCCACCGGGCCGGGGGGGGGCGGCa                                                                                                                                                                              >  2:59814/1‑229 (MQ=255)
                                                 gCTGGGAAAAACTTTGAGGCATGGGCGTGGATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGGGCCGGGGGCGAACATTTTGCCCGCAGACGGTATTGGCCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCaa                                                                                                                                                                             <  1:78787/228‑1 (MQ=255)
                                                              ttgaGTGATGGGCGTGAAGTAGGATGGCGTTTTTAGGGGGTATTTGCGCGTGCCGGTGGCGAACATTTTGCCCGCAGCCGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGa                                                                                                                                                                <  2:338659/225‑1 (MQ=255)
                                                              ttCAGTCCTGGGCGTGGATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCGGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGa                                                                                                                                                                <  1:99515/228‑1 (MQ=255)
                                                                                  aGGGTGGGGGTTTTAGGTAGTTTTTGGGCGGGCCGGGGGCGAACATTTTGCCCGGCGGCGGTATTGCCCCGGAGGCGGCCGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGc                                                                                                                                           <  2:221687/229‑1 (MQ=255)
                                                                                    gATGGCGGTTTTAGTGAGTTTTTGAGCGTGCCGGTGGCGAACATTTGGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCg                                                                                                                                          <  1:613240/228‑1 (MQ=255)
                                                                                                    ggtttttGGGCGTGCCGGGGGGTAACATTTTGCCCGCGGAAGTGTTTCACCGGTGGGGGGCAGCAGTGATTGGCCCGTTCGGTATGAGCGCTCAGGCGGTGCGGCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTACAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGgcgcg                                                                                                                         <  2:343889/225‑1 (MQ=255)
                                                                                                     gttttggggCGTGCCGGTGGCGAACATTTGGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGgcgcg                                                                                                                         <  1:277571/220‑1 (MQ=255)
                                                                                                        tttGTGCGTGGCGGGTGCGATCTTTCTGCCCGCAGACGGTATTGCCCCGGAGGCGCCAGCATTGATTTGACCTTGCGCTATTAGGGCTCATGCGGGGGGTCCCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGG‑GCGGGGCCAAACGGTCTGGGTGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGGCGCGCCGt                                                                                                                     <  2:334457/229‑1 (MQ=255)
                                                                                                                     ttGGCGTACTGTTTGCCCGCAGACCGTATTCACCCGGAGGGTGCCGCAGTTACTGAACCCTCCGCTATTACGGCTCATGCGGTGCGTGGGTGAGCCATTGCTCCCGGCGAGCAGGTGCTGGGGGTCGGG‑GCGGGGCCAAGCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGGACAGGTGGTGGTGGCGGATACCAGTCCGGCGCTCCGTGAACATGTGGCaa                                                                                                        <  2:263665/228‑1 (MQ=255)
                                                                                                                                   gCCCGCAGACGGTATTTACCCGCGGGTGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGGGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAAtt                                                                                           <  1:422387/228‑1 (MQ=255)
                                                                                                                                     ccGCAGGCGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCAGGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAATTAc                                                                                         <  1:114549/228‑1 (MQ=255)
                                                                                                                                                          gggcgggAGCTTTGATTGATTCGTTCGCTGTTGGCGCGCATGGGGGGGGTCGCGCAGCCATTGCTCCCGGCGGGCAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCg                                                                    <  2:339335/227‑1 (MQ=255)
                                                                                                                                                                       gattgaACCTTTCGATATTAGCGGGCATGCGGTGGCGCGGTCGGCCAGTGCTCCCGGGGGGCAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGATTTTGAcgc                                                       <  1:92887/228‑1 (MQ=255)
                                                                                                                                                                                   cGCTTTTAGCGCTCCTGCGGTGCGTCCCGCCACCATTGCCCCCGGCGAGCAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGATTTTGACGCGCAGCTACGGgcg                                          <  2:63167/229‑1 (MQ=255)
                                                                                                                                                                                    gCTATGAGCGCTCATGGGGTGCGGCGCGCAGCTGTTGCTCCGGGCGAGCAGGTGCTGGTGGTCGGG‑GCGGGGCCAAGCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGGCGCGACGTGAACATGTGGCAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGATTTTGACGCGCAGCTACGGgcg                                          <  1:41056/228‑1 (MQ=255)
                                                                                                                                                                                       aTTAGTGCTAATGCGGTGCGCAGCGAGGCCAGTGCGCCCGGCGAGCAGGTGCTGGTGGTAGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGATTTTGACGCGCAGCTACGGGCGCAg                                       <  1:563951/228‑1 (MQ=255)
                                                                                                                                                                                                    gggTGCGTCGGGCAGCCATTGCTCCCGGGGGGCAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCTGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCt                         <  1:123077/228‑1 (MQ=255)
                                                                                                                                                                                                       tGCGTCGGGGCGCCATTGCCCCCGGCGAGGAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGGTACCAGTCCGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCTgg                       <  1:204624/228‑1 (MQ=255)
                                                                                                                                                                                                                        gcccccGGCGAGCAGGCTCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAATTACCTTTACTg                                                                                 <  1:621343/150‑1 (MQ=255)
                                                                                                                                                                                                                         cgcccGGCGAGCGGGTGCTGGGGGTGGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGAcgc    <  1:424944/227‑1 (MQ=255)
                                                                                                                                                                                                                         cTCCCGGCGGGCAGGTGCTGGTGGTCGGG‑GCGGGGCCAAGCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTTGTGGTGGCGGATACCAGTCCGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGAcg     <  1:199256/228‑1 (MQ=255)
                                                                                                                                                                                                                           cccGGCGGGGGGGTGGTGGTGGTGGGG‑GGGGGGCAAATGGGGCTGGGCGCGGCGGCAATCGCTAAAGGCGATGGCGCACGGGTGGTGGTGGCGGATACCAGTCCGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGa  <  2:277318/229‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                              
ATGTTGCCTGTCAGCAATGCCCGGCGT‑GTAAAAGCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCGCAGACGGTATTGACCCGCAGGCGGCAGCATTGATTGAACCTTTCGCTATTAGCGCTCATGCGGTGCGTCGCGCAGCCATTGCTCCCGGCGAGCAGGTGCTGGTGGTCGGG‑GCGGGGCCAATCGGTCTGGGCGCGGCGGCAATCGCTAAAGCCGATGGCGCACAGGTGGTGGTGGCGGATACCAGTCCGGCGCGCCGTGAACATGTGGCAACGCGTCTGGAATTACCTTTACTGGACCCGTCAGCCGAAGATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGA  >  NC_000913/4596255‑4596701

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: ‑