Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 85,337 G→T 7.5% intergenic (+25/‑293) leuO → / → ilvI global transcription factor/acetolactate synthase 3 large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091385,3370GT 7.5% 184.8 / 5.9 53intergenic (+25/‑293)leuO/ilvIglobal transcription factor/acetolactate synthase 3 large subunit
Reads supporting (aligned to +/- strand):  ref base G (23/26);  new base T (2/2);  total (25/28)
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.35e-01

GTATTACCGCTGCCGTTAAAACAAAACAGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCA  >  NC_000913/85190‑85454
                                                                                                                                                   |                                                                                                                     
gTATTACCGCTGCCGTTAAAACAAAACAGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATg                                                                                                                       <  2:340782/148‑1 (MQ=255)
         cTGCCGTTAAAACAAAACAGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTgg                                                                                                              >  4:97746/1‑148 (MQ=255)
               ttAAAACAAAACAGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGaa                                                                                                                     <  4:269527/135‑1 (MQ=255)
               ttAAAACAAAACAGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGaa                                                                                                                     >  3:269527/1‑135 (MQ=255)
                tAAAACAAAACAGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTc                                                                                                       <  3:97746/148‑1 (MQ=255)
                  aaacaaaacaGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGTCGCGATAAAGGCCATCAGTGAATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTg                                                                                                     <  1:392069/148‑1 (MQ=255)
                   aacaaaacaGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGc                                                                                                    >  1:56779/1‑148 (MQ=255)
                     caaaacaGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCAt                                                                                                  >  1:95251/1‑148 (MQ=255)
                       aaacaGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGaa                                                                                                                     >  3:248680/1‑127 (MQ=255)
                       aaacaGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGaa                                                                                                                     <  4:248680/127‑1 (MQ=255)
                       aaacaGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCAttt                                                                                                <  3:821622/148‑1 (MQ=255)
                         acaGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGGTAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCAttttt                                                                                              >  2:406907/1‑148 (MQ=255)
                          cagcagAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCAtttttt                                                                                             <  2:440654/148‑1 (MQ=255)
                            gcagAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTc                                                                                                                 <  3:371565/126‑1 (MQ=255)
                            gcagAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTc                                                                                                                 >  4:371565/1‑126 (MQ=255)
                               gAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGtt                                                                                                            <  2:1125294/128‑1 (MQ=255)
                               gAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGtt                                                                                                            >  1:1125294/1‑128 (MQ=255)
                                      ttATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGt                                                                                                             <  4:228481/120‑1 (MQ=255)
                                      ttATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGt                                                                                                             >  3:228481/1‑120 (MQ=255)
                                       tATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGTATTTTTTATTGAATGTAGaa                                                                                >  4:461652/1‑148 (MQ=255)
                                         tctctcCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCAtttt                                                                                               >  1:138703/1‑131 (MQ=255)
                                         tctctcCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCAtttt                                                                                               <  2:138703/131‑1 (MQ=255)
                                           tctcCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAAtttt                                                                            <  2:95251/148‑1 (MQ=255)
                                            ctcCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTTCATTTTTTATTGAATGTAGAATTTTa                                                                           >  3:112568/1‑148 (MQ=255)
                                             tcCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGTATTTTTTATTGAATGTAGAATTTTAt                                                                          <  3:461652/148‑1 (MQ=255)
                                               cTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTc                                                                        <  2:56779/148‑1 (MQ=255)
                                                      gAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGt                                                                                                                  <  1:943502/99‑1 (MQ=255)
                                                      gAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGt                                                                                                                  >  2:943502/1‑99 (MQ=255)
                                                         gCTGCCGGGCGCGGTAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAAtgtgtg                                                              <  1:406907/148‑1 (MQ=255)
                                                         gCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTTCATTTTTTATTGAATGTAGAATTTTATTCTGAAtgtgtg                                                              <  4:112568/148‑1 (MQ=255)
                                                                  cgcgATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGt                                                                                                             >  3:862946/1‑92 (MQ=255)
                                                                  cgcgATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGt                                                                                                             <  4:862946/92‑1 (MQ=255)
                                                                      aTAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGa                                                                                 >  1:575220/1‑116 (MQ=255)
                                                                      aTAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGa                                                                                 >  3:1071450/1‑116 (MQ=255)
                                                                      aTAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGa                                                                                 <  4:1071450/116‑1 (MQ=255)
                                                                      aTAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGa                                                                                 <  2:575220/116‑1 (MQ=255)
                                                                      aTAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAAACCAGACATTCATGTCTGACCTATTCTGCaatc                                                                                                >  4:71176/1‑98 (MQ=18)
                                                                      aTAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAAACCAGACATTCATGTCTGACCTATTCTGCaatc                                                                                                <  3:71176/101‑4 (MQ=18)
                                                                        aaaGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGctc                                                         >  1:947223/1‑138 (MQ=255)
                                                                        aaaGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGctc                                                         <  2:947223/138‑1 (MQ=255)
                                                                               aTCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTa                                                                                            <  3:90710/96‑1 (MQ=255)
                                                                               aTCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTa                                                                                            >  4:90710/1‑96 (MQ=255)
                                                                                tCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTAtt                                                                                                        <  1:975916/83‑1 (MQ=255)
                                                                                tCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTAtt                                                                                                        >  2:975916/1‑83 (MQ=255)
                                                                                                        ctcAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGa                                                                     <  2:259140/94‑1 (MQ=255)
                                                                                                        ctcAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGa                                                                     >  1:259140/1‑94 (MQ=255)
                                                                                                        ctcAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAAtgtgtg                                                              <  2:414000/101‑1 (MQ=255)
                                                                                                        ctcAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAAtgtgtg                                                              >  1:414000/1‑101 (MQ=255)
                                                                                                            aTTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGa                                                                                                                      >  2:154829/1‑41 (MQ=255)
                                                                                                            aTTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGa                                                                                                                      <  1:154829/41‑1 (MQ=255)
                                                                                                                                            cAGACATTCATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCa  >  1:794510/1‑125 (MQ=255)
                                                                                                                                            cAGACATTCATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCa  <  2:794510/125‑1 (MQ=255)
                                                                                                                                                cATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTAtt                                                   <  3:28465/72‑1 (MQ=255)
                                                                                                                                                cATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTAtt                                                   >  4:28465/1‑72 (MQ=255)
                                                                                                                                                   |                                                                                                                     
GTATTACCGCTGCCGTTAAAACAAAACAGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCA  >  NC_000913/85190‑85454

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: