Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,590,223 T→G 26.1% T230P (ACC→CCC)  ugpA ← sn‑glycerol‑3‑phosphate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,590,2230TG26.1% 62.7 / 2.0 56T230P (ACC→CCC) ugpAsn‑glycerol‑3‑phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (10/26);  new base G (15/3);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 1.46e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGACGGTTCTCAATCATTGGTAACGCACCTTGCTTTCAACATAGCGGAACTGCACCACCGTCAGCACGATGACGAGGAACATCAACACCACCGACTGTGCGGCAGACGAAGCCAGATCCAGTCCGGTAAAACCTTCGCGGTAGATCTTATAAATCAGCGTCGTGGTGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGTGTCGAAGAAGGCATACACCAGGTTCACTACCAGCAGGAAGAAACTCACCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAGCGGCGAATCGGCCCTGCACCGTCGATGGCTGCGGCTTCGATCAACGAACGGGGAATGGATTGCAGCGCGGCATAGAAGAACAGGAAGTTGTAGCTGATTTGCTTCCATACTGAGGCAAACACCACCAGAAACATTGCCTGA  >  NC_000913/3590009‑3590446
                                                                                                                                                                                                                      |                                                                                                                                                                                                                                
cgggcggtttccaactcttGGTAACCCCCCTGTATTTCAAACAAGCGGAACGGCACCACCGTCAGCACGATGACGAGGAACATCAACACCACCGACTGTGCGGCAGACGAAGCCAGATCCAGTCCGGTAAAACCTTCGCGGTAGATCTTATAAATCAGCGTCGTGGTGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGTGTCGAAGAAGGCa                                                                                                                                                                                                                     <  1:80632/211‑1 (MQ=255)
       ttCTCAATCATTGGTAACGCACCTTGCTTTCAACATAGCGGAACTGCACCACCGTCAGCACGATGACGAGGAACATCAACACCACCGACTGTGCGGCAGACGAAGCCAGATCCAGTCCGGTAAAACCTTCGCGGTAGATCTTATAAATCAGCGTCGTGGGGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCAACGGGAAGGGGGCGAAGAAGGCGTACACCa                                                                                                                                                                                                              >  2:631178/1‑228 (MQ=255)
           caacacttGGTACCGCCACTTGCTTTAAACAGAGCGGAACTGCACCACCGTCAGCACGATGACGAGGAACATCAACACCACCGACGGTGCGGCAGACGAAGCCAGATCCGGTCCGGTAAAACCTTCGCGGTAGATCTTATAAATCAGCGTCGTGGTGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGTGTCGAAGAAGGCATACACCAGGtt                                                                                                                                                                                                          <  1:234650/222‑1 (MQ=255)
                 ttGGTAACGCACCTTGCTTTCAACATAGCGGAACTGCACCACCGTCAGCACGATGACGAGGAACATCAACACCACCGACTGTGCGGCAGACGAAGCCAGATCCAGTCCGGTAAAACCTTCGCGGTAGATCTTATAAATCAGCGTCGTGGTGGCCTGAACCGGCCCCCCGGACGCGGCGGCGTCGATCACCGGGCAGGGGGCGCAGAAGGCATACTACCGGTTCACTAc                                                                                                                                                                                                    >  2:95846/1‑228 (MQ=255)
                 ttGGTAACGCACCTTGCTTTCAACATAGCGGAACTGCACCACCGTCAGCACGATGACGAGGAACATCAACACCACCGACTGTGCGGCAGACGAAGCCAGATCCAGTCCGGTAAAACCTTCGCGGTAGATCTTATAAATCAGCGTCGTGGTGGCCTGAACCGACCCGCCGGACGTGGCGTGGTCCGTCTCCGGGAGGGGATCGAGAAGGGCAGCCCCCAGGtcgcccac                                                                                                                                                                                                    >  1:637350/1‑221 (MQ=255)
                 ttGGTAACGCACCTTGCTTTCAAAATAGCGGAACTGCACCACCGTCAGCACGATGACGAGGAACATCAACACCACCGACTGCGCGGCAGACGAAGCCAGATCCAGTCCGGTACAACCTTCGCGGTAGATCTTATAAATCAGCGCCGCGGTGGCCTGAACCGGCCCGCCGGACGTGGCGTCGTCGGTCCCCGGGACGGAGTTGAGCAAGGAGCACACAAGGtcacgtcc                                                                                                                                                                                                    >  1:66666/1‑221 (MQ=255)
                            ccTTGCTTTCAACATAGCGGAACTGCACCACCGTCAGCACGATGACGAGGAACATCAACACCACCGACTGTGCGGCAGACGAAGCCAGATCCAGTCCGGTAAAACCTTCGCGGTAGATCTTATAAATCAGCGTCGTGGTGGCCTGAACCGGCCCGCCGGCCGTGGCCGAGTCGGTCACCGGGACGGGGCAGAAAAAGGGATAACCCAGGGTCCCACCCAACAGgaaga                                                                                                                                                                                        >  1:451118/1‑228 (MQ=255)
                             gcggCTTTAAACATAGCGGAACTGCCCCCCCGTCAGCACGATGACGAGGAACATCAACACCACCGACTGGGCGGCAGACGAAGCCAGATCCAGTCCGGTAAAACCTTCGCGGTAGATCTTATAAATCAGCGTCGTGGTGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGTGTCGAAGAAGGCATACACCAGGTTCACTACCAGCAGgaagaa                                                                                                                                                                                       <  1:138271/225‑1 (MQ=255)
                             cTTGCTTTCAACATAGCGGAACTGCACCACCGTCAGCACGATGACGAGGAACATCAACACCACCGACTGTGCGGCAGACGAAGCCAGATCCAGTCCGGTAAAACCTTCGCGGTAGATCTTATAAATCAGCGTCGTGGTGGCCTGAACCGGCCCGCCGGACGGGGCGGCCTCCATCACCGTGAAGGGGTCGAAGAAGGCATAACCCAGGTCCACTACca                                                                                                                                                                                                  >  2:350508/1‑218 (MQ=255)
                                    tCAACATAGCGGAACTGCACCACCGTCAGCACGATGACGAGGAACATCAACACCACCGACTGTGCGGCAGACGAAGCCAGATCCAGTCCGGTAAAACCTTCGCGGTAGATCTTATAAATCAGCGTCGTGTTGGCCTGAACCGGACCGCCGGACGGGGCGGTGTCGATCCCCGGGAAGGTGTCGGAGAAGCCATACACCAGGTTACCTACCAGCAGGCAGaaaatacac                                                                                                                                                                                >  2:128720/1‑222 (MQ=255)
                                       aCATAGCGGAACTGCACCACCGTCAGCACGATGACGAGGAACATCAACACCACCGACTGTGCGGCAGACGAAGCCAGATCCAGTCCGGTAAAACCTTCGCGGTCGATCTTATAAATCAGCGTCGGGGTGGCCTGAACCGGCCCGCCGGACGTGGCGGCCGCGCTCCCCGGGGCGGGGTCGCAGCAGGCGGACCCCCGGTTCCCTCCCCCCAGGCAAGAAGCCACCggg                                                                                                                                                                             >  1:474901/1‑228 (MQ=255)
                                           aGCGGAACTGCACCACCGTCAGCACGATGACGAGGAACATCAACACCACCGACTGTGCGGCAGACGAAGCCAGATCCAGTCCGGTAAAACCTTCGCGGTAGATCTTATAAATCAGCGTCGTGGTGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGGGTCGAAGAAGGCATACACCAGGTTCACTACCAGCAGGAAGAAACTCACCGGGGCgg                                                                                                                                                                         >  2:189738/1‑227 (MQ=255)
                                                aaCTGCACCACCGTCAGCACGATGACGAGGAACATCAACACCACCGACTGTGCGGCAGACGAAGCCAGATCCAGTCCGGTAAAACCTTCGCGGTAGATCTTATAAATCAGCGTCGTGGTGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGGGTCGAAGAAGGCATACCCCAGGTTCACTACCAGCAGGAAGGAACGCACCGGGGCGGTAAg                                                                                                                                                                     >  2:492483/1‑227 (MQ=255)
                                                aaCTGCACCACCGTCAGCACGATGACGAGGAACATCAACACCACCGACTGTGCGGCAGACGAAGCCAGATCCAGTCCGGTAAAACCTTCGCGGTAGATCTTATAAATCAGCGTCGTGGTGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGGGTCGAAGAAGGCATACACCAGGTTCACTACCCGCCGGAAGAAACTCACCGGGGCGATAAGc                                                                                                                                                                    >  2:109991/1‑228 (MQ=255)
                                                          ccGTCAGCACGATGACGAGGAACATCAACACCACCGACTGTGCGGCAGACGAAGCCAGATCCAGTCCGGTAAAACCTTCGCGGTAGATCTTATAAATCAGCGTCGTGGTGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGTGTCgaag                                                                                                                                                                                                                           >  2:571960/1‑164 (MQ=255)
                                                           cGTCAGCACGATGACGAGGAACATCAACACCACCGACTGTGCGGCAGACGAAGCCAGATCCAGTCCGGTAAAACCTTCGCGGTAGATCTTATAAATCAGCGTCGTGGTGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGTGTCgaaga                                                                                                                                                                                                                          <  1:571960/164‑1 (MQ=255)
                                                               aGCACGATGACGAGGAACATCAACACCACCGACTGTGCGGCAGACGAAGCCAGATCCAGTCCGGTAAAACCTTCGCGGTAGATCTTATAAATCAGCGTCGTGGTGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCCCCGGGACGGGGTCGAAGAAGCCATACACCGGGGTCCCTACCCACAGGAAGAAACTCACCGGGGCGAAAAACgggcggccaaacttt                                                                                                                                                     >  1:69495/1‑215 (MQ=255)
                                                                  aCGATGACGAGGAACATCAACACCACCGACTGTGCGGCAGACGAAGCCAGATCCAGTCCGGTAAAACCTTCGCGGAAGATCTTATAAATCAGCGTCGTGGTGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGCAGGGGTCGAGGCAGGCATACACCACGCTCACTACCACGAGGAAGAAACTCACCAGGGCGACaacacgcaccccaacctcaaa                                                                                                                                                   >  1:1377/1‑208 (MQ=255)
                                                                    gATGACGAGGAACATCAACACCACCGACTGTGCGGCAGACGAAGCCAGATCCAGTCCGGTAAAACCTTCGCGGTAGATCTTATAAATCAGCGTCGGGGTGGCCTGAACCGGCCCGCCGGACGTGGAGGCGTCGACCACCGGGACGGGGTCGCAGAGGGCATACACCAGGTTCCCTACCAGCCGGAAGAAACGCACCGGGGGGAAAAACGGGAGAGCaaatgtaaaaaa                                                                                                                                                >  1:545641/1‑218 (MQ=255)
                                                                             gacatccAACCCCCCCCACTGGTCGGCAGACGAAGCCGGGTCCAGTCCGGTAAAACCTTCGCGGTAGATCTTATAAATCAGCGTCGGGGTGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGTGTCGAAGAAGGCATACACCAGGTTCACTACCAGCAGGAAGAAACTCACCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAGCGGCGAAt                                                                                                                                       <  1:189738/222‑1 (MQ=255)
                                                                                      caccacCGACTGTGCGGCAGACGAAGCCAGATCCAGTCCGGTAAAACCTTCGCGGTAGATCTTATAAATCAGCGTCGTGGTGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGTGTCGAAGAAGGCATACACCAGGTTCACTACCAGCAGGAAGAAACTCACCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAAGGGCGAATCGGCCCTGc                                                                                                                              >  2:699926/1‑228 (MQ=255)
                                                                                        ccaccGACTGTGCGGCAGACGAAGCCAGATCCAGTCCGGTAAAACCTTCGCGGTAGATCTTATAAATCAGCGTCGTGGTGGCATGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGGGTCGAAGAAGGCATACACCAGGTTCACTACAAGCAGGAAGAAACTCACCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAACGGCGAATCGGCCCTg                                                                                                                               >  2:63265/1‑225 (MQ=255)
                                                                                         caccGACTGTGCGGCAGACGAAGCCAGATCCAGTCCGGTAAAACCTTCGCGGTAGATCTTATAAATCAGCGTCGTGGTGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGGGTCGAAGAAGGCATACACCAGGTTCACTACCAGCAGGAAGAAACTCACCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAGCGGCGAATCGGCCCTGCAc                                                                                                                            >  2:282660/1‑227 (MQ=255)
                                                                                            cGACTGTGCGGCAGACGAAGCCAGATCCAGTCCGGTAAAACCTTCGCGGTAGATCTTATAAATCAGCGTCGTGGTGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGTGTCGAAGAAGGCATACACCAGGTTCACTACCAGCAGGAAGAAACTCACCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAGCGGCGAATCGGCCCTGCACCGt                                                                                                                         >  2:91976/1‑227 (MQ=255)
                                                                                              aCTGTGCGGCAGACGAAGCCAGATCCAGTCCGGTAAAACCTTCGCGGTAGATCTTATAAATCAGCGTCGTGGTGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGCGTCGAAGAAGGCATACACCAGGTTCACTACCAGCAGGAAGAAACTCACCGGGGCGATAACCGGCAGCGCAATCTTAAAAAAACAGCGAATCGGCCCCGCAccgccg                                                                                                                       >  2:595702/1‑227 (MQ=255)
                                                                                               cTGTGCGGCAGACGAAGCCAGATCCAGTCCGGTAAAACCTTCGCGGTAGATCTTATAAATCAGCGTCGTGGTGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGTGTCgaag                                                                                                                                                                                                                           <  2:234233/127‑1 (MQ=255)
                                                                                                 gtgCGGCAGACGAAGCCAGATCCAGTCCGGTAAAACCTTCGCGGTAGATCTTATAAATCAGCGTCGTGGTGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGTGTCGAAGAAGGCATACACCAGGTTCACTACCAGCCGGAAGAAACTCACCCGGGCGATAAGCGGCAGCGCAATCTTAAAGAAACGGCGAATCGGCCCTGAACCGTCGACgg                                                                                                                   >  2:264175/1‑228 (MQ=255)
                                                                                                  tgCGGCAGACGAAGCCAGATCCAGTCCGGTAAAACCTTCGCGGTAGATCTTATAAATCAGCGTCGTGGTGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGTGTCGAAGAAGGCATACACCAGGTTCACTACCAGCAGGAAGAAACTCACCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAGCGGCGAATCGGCCCTGAACCGTCGCTgg                                                                                                                   >  2:467194/1‑227 (MQ=255)
                                                                                                                      ccAGTCCGGTAAAACCTTCGCGGTAGATCTTATAAATCAGCGTCGTGGTGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGTGTCGAAGAAGGCATACACCAGGTTCACTACCAGCAGGAAGAAACTCACCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAGCGGCGAATCGGCCCTGCACCGTCGATGGCTGCGGCTTCGATCaacgaac                                                                                              <  2:69495/228‑1 (MQ=255)
                                                                                                                         ggccGGCAAAACCGTCGCGGTAGATCTTATAAATCAGCGGCGTGGTGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGTGTCGAAGAAGGCATACACCAGGTTCACTACCAGCAGGAAGAAACTCACCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAGCGGCGAATCGGCCCTGCACCGTCGATGGCTGCGGCTTCGATCAACGAACggg                                                                                           <  2:617874/226‑1 (MQ=255)
                                                                                                                          gctggTAAAACCTTCGCGTGAGTTCTTATAAAGCGGCGTCGGGGGGGCCTGAACCGGCCCGCCGGACGTGGCGGCGCCGATCACCGGGAAGGTGTCGAAGAAGGCATACACCAGGTTCACTACCAGCAGGAAGAAACGCCCCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAGCGGCGAATCGGCCCTGCACCGTCGATGGCTGCGGCTTCGATCAACGAACggg                                                                                           <  1:461387/224‑1 (MQ=255)
                                                                                                                           ccGGTAAAACCTGCGCGGTAGATCTTATAAATCAGCGTCGTGGTGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGTGTCGAAGAAGGCATACACCAGGTTCACTACCAGCAGGAAGAAACTCACCgggg                                                                                                                                                                            <  2:162142/145‑1 (MQ=255)
                                                                                                                                aaaaCCTTCGCGGTAGATCGTATAAATCAGCGGCGTGGTGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGTGGCGAAGAAGGCATACACCAGGTTCACTACCAGCAGGAAGAAACTCACCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAGCGGCGAATCGGCCCTGCACCGTCGATGGCTGCGGCTTCGATCAACGAACGGGGAATg                                                                                      <  2:498019/226‑1 (MQ=255)
                                                                                                                                        cggggTAGTTCTTAGAAATCAGCGTCGGGGGGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGTGTCGAAGAAGGCATACACCAGGTTCACTACCAGCAGGAAGAAACTCACCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAGCGGCGAATCGGCCCTGCACCGTCGATGGCTGCGGCTTCGATCAACGAACGGGGAATGGATTGCAgcg                                                                            <  1:631178/225‑1 (MQ=255)
                                                                                                                                        cgcgGTAGACCTTATAAATCAGCGTCGTGGTGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGTGTCGAAGAAGGCATACACCAGGTTCACTACCAGCAGGAAGAAACTCACCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAGCGGCGAATCGGCCCTGCACCGTCGATGGCTGCGGCTTCGATCAACGAACGGGGAATGGATTGCAgcg                                                                            <  2:276564/228‑1 (MQ=255)
                                                                                                                                        cgcgGGGGATCTTATAAGTCGGGGTGCTGGTGGCCTGAACCGGCCCGGCGGACGTGGAGGCGTCGATCACCGGGAAGGTGTCGAAGAAGGCATACACCAGGTTCACTACCAGCAGGAAGAAACTCACCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAGCGGCGAATCGGCCCTGCACCGTCGATGGCTGCGGCTTCGATCAACGAACGGGGAATGGATTGCAgcg                                                                            <  1:467194/228‑1 (MQ=255)
                                                                                                                                         gcgGTAGATCTTATAAATCAGCGTCGTGGTGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGTGTCGAAGAAGGCATACACCAGGTTCACTACCAGCAGGAAGAAACTCACCGGGGCGATAAGCGGCAGGGCAAACTTAAAGAAGCGGCGAATCGGCCCTCAACCGGCGATGGCGGCTGCTTCGATCAACGAAAGGGGGATGGATTGCAgcgg                                                                           >  2:699979/1‑227 (MQ=255)
                                                                                                                                                cttttATAGATCGCCGTCTGGGTGGCGTGAGCCGGCCCGCCGGACGTGGCGGCGTCGATCCCCGGGAAGGTGCCGAAGAAGGCATACACCAGGTTCACTACCAGCAGGAAGAAACTCACCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAGCGGCGAATCGGCCCTGCACCGTCGATGGCTGCGGCTTCGATCAACGAACGGGGAATGGATTGCAGCGCGGCATag                                                                    <  1:282660/225‑1 (MQ=255)
                                                                                                                                                             gCGTCGTGGTGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGGGTCGAAGAAGGCATACACCAGGGTCACTACCAGCAGGAAGAGACTCACCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAGAGGCGAATCGGCCCCGCACCGGCGCGGGCTGCGGCTTCGATCAACGAACGGGGGATGGAATGCGGCGCGGCAACGAACAACAGGAAga                                                       >  1:685302/1‑227 (MQ=255)
                                                                                                                                                              cgtcgtGGTGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGTGTCGAAGAAGGCATACACCAGGTTCACTACCAGCAGGAAGAAACTCACCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAGCGGCGAATCGGCCCTGAACCGTCGATGGCTGCGGCTTCGGTCAACGAACGGGGAATGGATTGCAGCGCGGCATggaggcggctcttttt                                                      >  2:15663/1‑215 (MQ=255)
                                                                                                                                                                   tggtggCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGTGTCGAAGAAGGCATACACCAGGTTCACTACCAGCAGGAAGAAACTCACCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAGCGGCGAATCGGCCCTGCACCGTCGATGGCTGCGGCTTCGATCAACGAACGGGGAATGGATTGCAGCGCGGCATAGAAGAACAGGAACTTGTAg                                                  >  2:148295/1‑227 (MQ=255)
                                                                                                                                                                        tcataaaCCCGCCCGCCGGACGGTGCGGAGTCGATCACCGGGAAGGGGTCGAAGAAGGCATACACCAGGTTCACTACCAGCAGGAAGAAACTCACCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAGCGGCGAATCGGCCCTGCACCGTCGATGGCTGCGGCTTCGATCAACGAACGGGGAATGGATTGCAGCGCGGCATAGAAGAACAGGAAGTTGTAGCTGAtt                                            <  2:637350/223‑1 (MQ=255)
                                                                                                                                                                        gtgggAACCGCCCCGCCGGACGGGGGGGCGTCCGCCAACGGGAAGGTGTCGAAGAAGGCATGCACCGGGTTCAATACCAGCAGGAAGAACCTCACCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAGCGGCGAATCGGCCCTGCACCGTCGATGGCTGCGGCTTCGATCAACGAACGGGGAATGGATTGCAGCGCGGCATAGAAGAACAGGAAGTTGTAGCTGAtt                                            <  1:95846/224‑1 (MQ=255)
                                                                                                                                                                        ggcTGGACCGGCCCGCCGGGCGGGGCGGCGTCGATCACCGGGAAGGGCGCGAAGGAGGCATGCACCAGGCTCACTACCAGCAGGAAGCAACTCACCGGGGCGGAAAGCGGCAGCGCAATCTTAAAGAAGCGGCGAATCGGCCCTGCACCGTCGATGGATGCGGCTTCGATCAACGAACGGGGAATGGATTGCAGCGCGGCATAGAAGAACAGGAAGTTGTAGCTGAtt                                            <  1:304309/226‑1 (MQ=255)
                                                                                                                                                                         ccTGAATCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGTGTAGAAGAAGGCATACACCAGGTTCACTACCAGCAGGAAGAAACTCACCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAGCGGCGAATCGGCCCTGCACCGTCGATGGCTGCGGCTTCGATCAACGAACGGGGAATGGATTGCAGCGCGGCATAGAAGAACAGGAAGTTGTAGCTGAtt                                            <  2:545641/227‑1 (MQ=255)
                                                                                                                                                                              gccGGGCCGCCGGACTTGGCGGCGTCGATCGCCGGGAAGGTGTCGAAGAAGGCATACACCAGGTTCACTACCAGCAGGAAGAAACTCACCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAGCGGCGAATCGGCCCTGCACCGTCGATGGCTGCGGCTTCGATCAACGAACGGGGAATGGATTGCAGCGCGGCATAGAAGAACAGGAAGTTGTAGCTGATTTGCtt                                       <  2:222008/226‑1 (MQ=255)
                                                                                                                                                                                   cccGCCGGACCGGGCGGCGTCGATCCCCGGGAAGGTGTCGAAGAAGGCATACACCAGGTTCACTAACAGCAGGAAGAAACTCACCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAGCGGCGAATCGGCCCTGCACCGTCGATGGCTGCGGCTTCGATCAACGAACGGGGAATGGATTGCAGCGCGGCATAGAAGAACAGGAAGTTGTAGCTGATTTGCTTCCATa                                  <  2:622427/227‑1 (MQ=255)
                                                                                                                                                                                     cgccgGACGGGGGGGCGTGGATCACCGGGAAGGTGTCAAAGAAGGCAAACACCAGGTTCACTACCAGCAGGAAGAAACTCACCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAGCGGCGAATCGGCCCTGCACCGTCGATGGCTGCGGCTTCGATCAACGAACGGGGAATGGATTGCAGCGCGGCATAGAAGAACAGGAAGTTGTAGCTGATTTGCTTCCATACt                                <  1:595702/227‑1 (MQ=255)
                                                                                                                                                                                         gggcGTGGCGGGGTCGAACCCCGGGAAGGTGTAGGAAGAGGCATACACCAGGTTCCCTACCAGCAGGAAGAAACTCACCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAGCGGCGAATCGGCCCTGCACCGTCGATGGCTGCGGCTTCGATCAACGAACGGGGAATGGATTGCAGCGCGGCATAGAAGAACAGGAAGTTGTAGCTGATTTGCTTCCATACTGAgg                            <  1:91976/224‑1 (MQ=255)
                                                                                                                                                                                            ctggggggcggcGATCCGCGGGAAGGTTGCGAAGAAGGCCTCCACCAGGTTCACTCCCAGCAGGAAGCAACTCCCCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAGCGGCGAATCGGCCCTGCACCGTCGATGGCTGCGGCTTCGATCAACGAACGGGGAATGGATTGCAGCGCGGCATAGAAGAACAGGAAGTTGTAGCTGATTTGCTTCCATACTGAGGCaaa                        <  1:699979/222‑1 (MQ=255)
                                                                                                                                                                                            agTGGCGGCGTCGATACGCGGTAAGGTGTCCAAGAAGGCCTACAACAGGTTCACTACCAGCAGGAAGCAACTCACCGGGGCGATAAGCGGCAGCGCCATCTTAAAGAAGCGGCGAATCGGCCCTGCACCGTCGATGGCTGCGTCTTCGATCAACGAACGGGGAATGGATTGCAGCGCGGCATAGAAGAACAGGAAGTTGTAGCTGATTTGCTTCCATACTGAGGCaaa                        <  2:486554/227‑1 (MQ=255)
                                                                                                                                                                                                    cgTCGATCACCGGGAAGGTGTCGAAGAAGGCATACACCAGGTTCACTACCAGCAGGAAGAAACTCACCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAGCGGCGAATCGGCCCTGCACCGTCGATGGCTGCGGCTTCGATCAACGAACGGGGAATGGATTGCAGCGCGGCATAGAAGAACAGGAAGTTGTAGCTGATTTGCTTCCTTACTGAGGCAAAcaccac                  >  2:471626/1‑226 (MQ=255)
                                                                                                                                                                                                       cgacccccGGGAAGGGGTCGAAGAAGGCATACACCGGGTTCACTACCAGCAGGAAGAAACTCACCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAGCGGCGAATCGGCCCTGCACCGTCGATGGCTGCGGCTTCGATCAACGAACGGGGAATGGATTGCAGCGCGGCATAGAAGAACAGGAAGTTGTAGCTGATTTGCTTCCATACTGAGGCAAACACCACCAGa              <  1:471626/221‑1 (MQ=255)
                                                                                                                                                                                                         aacACCGGGAAGGTGTCGAAGAAGGCATACACCAGGTTCAATACCAGCAGGAAGAAACTCCCCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAGCGGCGAATCGGCCCTGCACCGTCGATGGCTGCGGCTTCGATCAACGAACGGGGAATGGATTGCAGCGCGGCATAGAAGAACAGGAAGTTGTAGCTGATTTGCTTCCATACTGAGGCAAACACCACCAGaaa            <  2:474901/225‑1 (MQ=255)
                                                                                                                                                                                                                  aaGGTGTCGAAGAAGGCAATCACCAGGTCCACTACGCGCAGGAAGAAACGCACCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAGCGGCGAATCGGCCCTGCACCGTCGATGGCTGCGGCTTCGATCAACGAACGGGGAATGGATTGCAGCGCGGCATAGAAGAACAGGAAGTTGTAGCTGATTTGCTTCCATACTGAGGCAAACACCACCAGAAACATTGCCTGa  <  2:629716/228‑1 (MQ=255)
                                                                                                                                                                                                                  aaGGTGGCGAAGAGGGCATTCCCCAGGTTCACTACCAGAAGGAAGAAGCTCACCGGGGCGATAAGGGGCAGCGCAATCGTAAAGAAGCGGCGAAGCGGCACTGCACCGGCGATGGCTGCGGCTTCGATCAACGAACGGGGAATGGATTGCAGCGCGGCATAGAAGAACAGGAAGTTGTAGCTGATTTGCTTCCATACTGAGGCAAACACCACCAGAAACATTGCCTGa  <  1:667167/228‑1 (MQ=255)
                                                                                                                                                                                                                      |                                                                                                                                                                                                                                
CGGACGGTTCTCAATCATTGGTAACGCACCTTGCTTTCAACATAGCGGAACTGCACCACCGTCAGCACGATGACGAGGAACATCAACACCACCGACTGTGCGGCAGACGAAGCCAGATCCAGTCCGGTAAAACCTTCGCGGTAGATCTTATAAATCAGCGTCGTGGTGGCCTGAACCGGCCCGCCGGACGTGGCGGCGTCGATCACCGGGAAGGTGTCGAAGAAGGCATACACCAGGTTCACTACCAGCAGGAAGAAACTCACCGGGGCGATAAGCGGCAGCGCAATCTTAAAGAAGCGGCGAATCGGCCCTGCACCGTCGATGGCTGCGGCTTCGATCAACGAACGGGGAATGGATTGCAGCGCGGCATAGAAGAACAGGAAGTTGTAGCTGATTTGCTTCCATACTGAGGCAAACACCACCAGAAACATTGCCTGA  >  NC_000913/3590009‑3590446

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: