Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,952,380	G→A	L60L (CTG→CTA)	ilvM →	acetolactate synthase 2 small subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,952,380	0	G	A	97.5%	267.8 / ‑2.2	81	L60L (CTG→CTA)	ilvM	acetolactate synthase 2 small subunit
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (40/39);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01

TTTCCACGTCTGCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAGCGCAAAAGGAA  >  NC_000913/3952284‑3952476                                                                                                 |                                                                                                 tttCCACGTCTGCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTAttt                                                                                               <  2:228045/100‑1 (MQ=255) tttCCACGTCTGCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTAttt                                                                                               <  1:392082/100‑1 (MQ=255)   tCCACGTCTGCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAg                                                                                             <  2:1086710/100‑1 (MQ=255)   tCCACGTCTGCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAg                                                                                             >  1:422359/1‑100 (MQ=255)   tCCACGTCTGCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAg                                                                                             <  2:343635/100‑1 (MQ=255)     cACGTCTGCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTc                                                                                           <  2:922434/100‑1 (MQ=255)       cGTCTGCTCACTGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGGCCACGGTCGGTCGACTTACTATTTAGTCAg                                                                                         <  1:1191045/100‑1 (MQ=255)        gTCTGCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAg                                                                                         >  1:1533634/1‑99 (MQ=255)         tCTGCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGtt                                                                                       <  1:1431824/100‑1 (MQ=255)          cTGCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTta                                                                                      <  1:78417/100‑1 (MQ=255)          cTGCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTta                                                                                      >  2:1005819/1‑100 (MQ=255)           tGCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTta                                                                                      <  2:225220/99‑1 (MQ=255)            gCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTtaaa                                                                                    >  2:1264070/1‑100 (MQ=255)            gCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTtaaa                                                                                    <  2:19887/100‑1 (MQ=255)            gCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTtaaa                                                                                    >  2:1709728/1‑100 (MQ=255)             cTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTtaaat                                                                                   <  1:1360215/100‑1 (MQ=255)               cAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTtaaataa                                                                                 >  1:36437/1‑100 (MQ=255)                aaTGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTtaaataaa                                                                                >  2:1115511/1‑100 (MQ=255)                 aTGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAAc                                                                               <  1:1305696/100‑1 (MQ=255)                     atatGGCCGCCGCCTGCGATGCATAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACtgg                                                                            <  2:132196/99‑1 (MQ=255)                     atatGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACtggt                                                                           <  1:1042738/100‑1 (MQ=255)                      tatGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACtggtg                                                                          <  2:32280/100‑1 (MQ=255)                       atGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACtggtgg                                                                         <  1:26403/100‑1 (MQ=255)                       atGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACtggtgg                                                                         <  1:1014176/100‑1 (MQ=255)                       atGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACtggtg                                                                          >  1:772580/1‑99 (MQ=255)                          gccgccgccAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACg                                                                      >  1:1654651/1‑100 (MQ=255)                          gccgccgccAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACg                                                                      >  1:288429/1‑100 (MQ=255)                           ccgccgccAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGt                                                                     <  1:1412863/100‑1 (MQ=255)                             gccgccAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCg                                                                   <  2:1675658/100‑1 (MQ=255)                               cgccAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGca                                                                 <  1:1335643/100‑1 (MQ=255)                               cgccAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGca                                                                 >  2:778434/1‑100 (MQ=255)                                gccAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGcac                                                                <  2:1017268/100‑1 (MQ=255)                                   aGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACg                                                             >  2:1690074/1‑100 (MQ=255)                                       aTGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGc                                                         >  1:1134914/1‑100 (MQ=255)                                       aTGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGc                                                         >  2:17319/1‑100 (MQ=255)                                        tGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGcc                                                        <  1:205780/100‑1 (MQ=255)                                        tGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGcc                                                        <  2:422361/100‑1 (MQ=255)                                        tGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGAACACGTTGcc                                                        >  1:39267/1‑100 (MQ=255)                                         gCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCa                                                       <  1:617025/100‑1 (MQ=255)                                         gCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCa                                                       <  2:1677837/100‑1 (MQ=255)                                          cacaAAATATAAATATCGAATTGACCGTTGACAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCAGACGTTGCCAt                                                      >  1:1192262/1‑100 (MQ=255)                                           acaAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATc                                                     >  2:1671287/1‑100 (MQ=255)                                           acaAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATc                                                     >  1:1262637/1‑100 (MQ=255)                                             aaaaTATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTg                                                   <  2:1630079/100‑1 (MQ=255)                                             aaaaTATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTg                                                   <  1:996996/100‑1 (MQ=255)                                                 tataaatatCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCag                                               <  2:790453/100‑1 (MQ=255)                                                      atatCGAATTGACCGTTGCCAGCCAACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGcac                                          >  1:1264810/1‑100 (MQ=255)                                                        atCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACaa                                        >  2:922403/1‑100 (MQ=255)                                                         tCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAAc                                       <  1:221437/100‑1 (MQ=255)                                                         tCGAATTGACCGTTGCCACCCCACGGTCGGTCGACTTACTATTGAGTCAGTTACATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAAc                                       <  1:1162525/100‑1 (MQ=255)                                                          cGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAAcc                                      <  1:773123/100‑1 (MQ=255)                                                          cGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAAcc                                      >  2:1126399/1‑100 (MQ=255)                                                           gAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACca                                     <  1:778429/100‑1 (MQ=255)                                                           gAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACca                                     >  2:1814746/1‑100 (MQ=255)                                                             aTTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACcaca                                   >  2:1301580/1‑100 (MQ=255)                                                                gACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATca                                >  2:991554/1‑100 (MQ=255)                                                                gACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATca                                <  2:1262648/100‑1 (MQ=255)                                                                 aCCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATcac                               >  1:1468243/1‑100 (MQ=255)                                                                  ccGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATcac                               >  1:134309/1‑99 (MQ=255)                                                                    gTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCacaac                            <  2:288431/100‑1 (MQ=255)                                                                     ttGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCacaac                            <  1:1054118/99‑1 (MQ=255)                                                                      tGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCAcaacaa                          >  2:1454080/1‑100 (MQ=255)                                                                      tGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCAcaacaa                          >  2:1767695/1‑100 (MQ=255)                                                                         cAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAAt                        >  1:366126/1‑99 (MQ=255)                                                                             ccccgatcggtcgACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGCGCACAACCACATCACAAGAAATCcg                     <  1:613814/92‑1 (MQ=255)                                                                             ccACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCcgcg                   >  2:1484397/1‑100 (MQ=255)                                                                              cACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGAGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCgcgc                  <  1:1277297/100‑1 (MQ=255)                                                                               aCGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGcc                 >  1:1747479/1‑100 (MQ=255)                                                                                cggtcggtCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGGTGCCATCTGCCATAGCACAACCACATCACAACAAATCCGCGCCt                >  2:989733/1‑100 (MQ=255)                                                                                 ggtcggtcGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTg               <  1:922397/100‑1 (MQ=255)                                                                                  gtcggtcgACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTg               <  1:567762/99‑1 (MQ=255)                                                                                    cggtcgACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAgc            >  2:1708375/1‑100 (MQ=255)                                                                                     ggtcgACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAgcg           >  2:643625/1‑100 (MQ=255)                                                                                         gACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAGCGCaaa       >  1:145142/1‑100 (MQ=255)                                                                                           cTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAGCGCAAAAg     >  2:1594852/1‑100 (MQ=255)                                                                                           cTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAGCGCAAAAg     >  2:19392/1‑100 (MQ=255)                                                                                           cTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAGCGCAAAAg     >  2:1508289/1‑100 (MQ=255)                                                                                           cTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAGCGCAAAAg     <  2:1533645/100‑1 (MQ=255)                                                                                            ttACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAGCGCAAAAgg    <  2:1042382/100‑1 (MQ=255)                                                                                             tACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAGCGCAAAAGGa   >  1:1088293/1‑100 (MQ=255)                                                                                              aCTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAGCGCAAAAGGaa  >  1:541768/1‑100 (MQ=255)                                                                                                 |                                                                                                 TTTCCACGTCTGCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAGCGCAAAAGGAA  >  NC_000913/3952284‑3952476

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