Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,704,800 A→G 7.0% F348L (TTT→CTT)  crtQ ← glycosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,704,8000AG 7.0% 155.2 / 3.1 57F348L (TTT→CTT) crtQglycosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (30/23);  new base G (2/2);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

CATACTAGTCCTCCTATATTGAAATGCCATTACCTTGATTTATTGTTCTTTAGATATAGAATATTGACGACCTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATGCAGATAAATAAATTGTGTAGTATAGCTCAAGTAAAGTGCTACCATTTTTCTTACAGACATTTGGCGATATTTTAAAC  >  CP000730/2704666‑2704940
                                                                                                                                      |                                                                                                                                            
cATACTAGTCCTCCTATATTGAAATGCCATTACCTTGATTTATTGTTCTTTAGATATAGAATATTGACGACCTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACat                                                                                                                                        >  2:32460/1‑141 (MQ=255)
 aTACTAGTCCTCCTATATTGAAATGCCATTACCTTGATTTATTGTTCTTTAGATATAGAATATTGACGACCTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACata                                                                                                                                       <  2:121556/141‑1 (MQ=255)
   aCTAGTCCTCCTATATTGAAATGCCATTACCTTGATTTATTGTTCTTTAGATATAGAATATTGACGACCTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACataaa                                                                                                                                     >  2:407726/1‑141 (MQ=255)
    cTAGTCCTCCTATATTGAAATGCCATTACCTTGATTTATTGTTGTGTAGATATAGAATATTGACGACCTTTCCATTCAACTACACCATACCGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACataaat                                                                                                                                    <  2:17721/141‑1 (MQ=255)
    cTAGTCCTCCTATATTGAAATGCCATTACCTTGATTTATTGTTCTTTAGATATAGAATATTGACGACCTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACataaat                                                                                                                                    <  2:409925/141‑1 (MQ=255)
       gTCCTCCTATATTGAAATGCCATTACCTTGATTTATTGTTCTTTAGATATAGAATATTGACGACCTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAAc                                                                                                                                 <  1:443671/141‑1 (MQ=255)
          ctcctATATTGAAATGCCATTACCTTGATTTATTGTTCTTTAGATATAGAATATTGACGACCTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAAt                                                                                                                              <  1:94490/141‑1 (MQ=255)
              tataTTGAAATGTCATTCCCTCGATTTATTGTTTTTAAGATATAGAATATTTACGACCTTTCCATTCAACTACATCATAACGATGCGTTTGTTTCCAGGATTGCATAAAAATTTTAGTAAAAATCATAAATAACAATGGAt                                                                                                                          >  2:384879/1‑141 (MQ=255)
               atatTGAAATGCCATTACCTTGATTTATTGTTCTTTAGATATAGAATATTGACGACCTTTCCATTCAACTACACCATAACGATGCGTTTTTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATg                                                                                                                         <  1:175929/141‑1 (MQ=255)
                 atTGAAATGCCATTACCTTGATTTATTGTTCTTTAGATATAGAATATTGACGACCTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGAc                                                                                                                       <  1:19791/141‑1 (MQ=255)
                  ttGAAATGCCATTACCTTGATTTATTGTTCTTTAGATATAGAATATTGACGACCTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACa                                                                                                                      >  1:10643/1‑141 (MQ=255)
                    gAAATGCCATTACCTTGATTTATTGTTCTTTAGATATAGAATATTGACGACCTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATa                                                                                                                    >  2:217013/1‑141 (MQ=255)
                     aaaTGCCATTACCTTGATTTATTGTTCTTTAGATATAGAATATTGACGACCTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGAGGACATAc                                                                                                                   >  2:113490/1‑141 (MQ=255)
                           cATTACCTTGATTTATTGTTCTTTAGATATAGAATATTGACGACCTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGAGGACATACCATtaa                                                                                                             >  1:320773/1‑141 (MQ=255)
                                   tGATTTATTGTTCTTTAGATATAGAATATTGACGACCTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAg                                                                                                     <  2:43696/141‑1 (MQ=255)
                                    gATTTATTGTTCTTTAGATATAGAATATTGACGACCTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGa                                                                                                    >  2:121236/1‑141 (MQ=255)
                                        tATTGTTCTTTAGATATAGAATATTGACGACCTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAAt                                                                                                >  1:238640/1‑141 (MQ=255)
                                          ttGTTCTTTAGATATAGAATATTGACGACCTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTg                                                                                              >  1:219661/1‑141 (MQ=255)
                                          ttGTTCTTTAGATATAGAATATTGACGACCTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTg                                                                                              >  1:261247/1‑141 (MQ=255)
                                          ttGTTCTTTAGATATAGAATATTGACGACCTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTg                                                                                              >  1:20344/1‑141 (MQ=255)
                                               cTTTAGATATAGAATATTGACGACCTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGa                                                                                         <  1:291232/141‑1 (MQ=255)
                                                  tAGATATAGAATATTGACGACCTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGAccc                                                                                      >  2:424395/1‑141 (MQ=255)
                                                  tAGATATAGAATATTGACGACCTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGAccc                                                                                      <  1:92549/141‑1 (MQ=255)
                                                    aatatAGAATATTGACGACCTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCtt                                                                                    <  2:124217/140‑1 (MQ=255)
                                                      tataGAATATTGACGACCTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCg                                                                                  >  1:427875/1‑141 (MQ=255)
                                                         aGAATATTGACGACCTTTCCATTGAACTACACTATAACGAGGTGTTGGTTTCCAAGACTGGATGAAAATTTGAGGAAAAAACATAAATAACAAATGATGACATAACATTAAAAAAATAAAATATGGTCCGACGCTTccgtg                                                                               >  2:59481/1‑137 (MQ=255)
                                                         aGAATATTGACGACCTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAGAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATg                                                                               >  1:424763/1‑141 (MQ=255)
                                                         aGAATATTGACGACCTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAGAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATg                                                                               >  1:424535/1‑141 (MQ=255)
                                                                tGACGACCTTTCCATTCAACTACACCATAACTATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATGCAGataa                                                                        <  2:194957/141‑1 (MQ=255)
                                                                tGACGACCTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATGCAGataa                                                                        <  1:229284/141‑1 (MQ=255)
                                                                     acCTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATGCAGAtaaataaa                                                                   <  2:133516/141‑1 (MQ=255)
                                                                      ccTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATGCAGATaaataaat                                                                  >  2:261052/1‑141 (MQ=255)
                                                                      ccTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATGCAGATaaataaat                                                                  >  1:300623/1‑141 (MQ=255)
                                                                      ccTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATGCAGATaaataaat                                                                  <  1:362812/141‑1 (MQ=255)
                                                                          tCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATGCAGATAAATAAATtgtg                                                              >  2:365791/1‑141 (MQ=255)
                                                                           ccATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATGCAGATAAATAAATTgtgt                                                             >  2:219175/1‑141 (MQ=255)
                                                                               tCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATGCAGATAAATAAATTGTgtagta                                                         <  2:331801/141‑1 (MQ=255)
                                                                                 aaCTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAACTGGCCGACCCTTCGATACAAATAAATAAATTGTGTTGtata                                                       >  2:451542/1‑141 (MQ=255)
                                                                                            aaCGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATGCAGATAAATAAATTGTGTAGTATAGCTCAAGTaaa                                            >  2:82572/1‑141 (MQ=255)
                                                                                              cGATGCGTTTGTTTCCAAGATTGGTTGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGACCGACCCTTCGATGCAGATAAATAAATTGTGTAGTATAGCTCAAGTAAAGt                                          >  2:108424/1‑141 (MQ=255)
                                                                                                       tgtttCCAAGATTGGATGAAAATTTTAGTAAGAAACATAAATATCAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATGCAGATAAATAAATTGTGTAGTATAGCTCAAGTAAAGTGCTACCAtt                                 <  2:424763/141‑1 (MQ=255)
                                                                                                       tgtttCCAAGATTGGATGAAAATTTTAGTAAGAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATGCAGATAAATAAATTGTGTAGTATAGCTCAAGTAAAGTGCTACCAtt                                 <  2:424535/141‑1 (MQ=255)
                                                                                                       tgatgccAAGAGTGGATGAAAATATTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCTACCCTTCGATGCAGATAACTAAATTGTGTAGGATAGCTCAAGTAAAGTGCTACCCtt                                 <  2:259555/136‑1 (MQ=255)
                                                                                                        gtttCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATGCAGATAAATAAATTGTGTAGTATAGCTCAAGTAAAGTGCTACCAttt                                >  2:255845/1‑141 (MQ=255)
                                                                                                        gtttCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATGCAGATAAATAAATTGTGTAGTATAGCTCAAGTAAAGTGCTACCAttt                                >  1:145947/1‑141 (MQ=255)
                                                                                                         tttCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATGCAGATAAATAAATTGTGTAGTATAGCTCAAGTAAAGTGCTACCAtttt                               <  1:121236/141‑1 (MQ=255)
                                                                                                             cAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATGCAGATAAATAAATTGTGTAGTATAGCTCAAGTAAAGTGCTACCATTTTTCtt                           >  1:351840/1‑141 (MQ=255)
                                                                                                              aaGATTGGATGAAAATTTTAGTAAAAAACATAAATAACACTGGATGACATACCATTAATGAATTAGAAAATTGGCCGACCCTTCGATGCAGATAAATAAATTGTGTAGTATAGCTCAAGTAAAGTGCTACCATTTTTCTTa                          <  2:263428/141‑1 (MQ=255)
                                                                                                              aaGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATGCAGATAAATAAATTGTGTAGTATAGCTCAAGTAAAGTGCTACCATTTTTCTTa                          >  1:428926/1‑141 (MQ=255)
                                                                                                                     gATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATGCAGATAAATAAATTGTGTAGTATAGCTCAAGTAAAGTGCTACCATTTTTCTTACAGACAt                   <  1:179296/141‑1 (MQ=255)
                                                                                                                     gATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATGCAGATAAATAAATTGTGTAGTATAGCTCAAGTAAAGTGCTACCATTTTTCTTACAGACAt                   >  2:401421/1‑141 (MQ=255)
                                                                                                                     gATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATGCAGATAAATAAATTGTGTAGTATAGCTCAAGTAAAGTGCTACCATTTTTCTTACAGACAt                   >  1:331083/1‑141 (MQ=255)
                                                                                                                       tGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATGCAGATAAATAAATTGTGTAGTATAGCTCAAGTAAAGTGCTACCATTTTTCTTACAGACAttt                 >  2:64012/1‑141 (MQ=255)
                                                                                                                       tGAAAATTATAGTAAAAAACATAAATAACCATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATGCAGATAAATAAATTGTGTAGTATAGCTCAAGTAAAGTGCTACCATTTTTCTTACAGACAttt                 <  2:118408/141‑1 (MQ=255)
                                                                                                                            aTTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATGCAGATAAATAAATTGTGTAGTATAGCTCAAGTAAAGTGCTACCATTTTTCTTACAGACATTTGGCGa            <  1:156457/141‑1 (MQ=255)
                                                                                                                              tttAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATGCAGATAAATAAATTGTGTAGTATAGCTCAAGTAAAGTGCTACCATTTTTCTTACAGACATTTGGCGata          <  1:436578/141‑1 (MQ=255)
                                                                                                                                tAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGGCGACCCTTCGATGCAGATAAATAAATTGTGTAGTATAGCTCAAGTAAAGTGCTACCATTTTTCTTACAGACATTTGGCGATAtt        >  2:447973/1‑141 (MQ=255)
                                                                                                                                     aaaaaCATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATGCAGATAAATAAATTGTGTAGTATAGCTCAAGTAAAGTGCTACCATTTTTCTTACAGACATTTGGCGATATTTTaaa   <  2:427875/141‑1 (MQ=255)
                                                                                                                                      aaaaCATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATGCAGATAAATAAATTGTGTAGTATAGCTCAAGTAAAGTGCTACCATTTTTCTTACAGACATTTGGCGATATTTTAAAc  <  2:238640/141‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                            
CATACTAGTCCTCCTATATTGAAATGCCATTACCTTGATTTATTGTTCTTTAGATATAGAATATTGACGACCTTTCCATTCAACTACACCATAACGATGCGTTTGTTTCCAAGATTGGATGAAAATTTTAGTAAAAAACATAAATAACAATGGATGACATACCATTAATAAATTAGAAAATTGGCCGACCCTTCGATGCAGATAAATAAATTGTGTAGTATAGCTCAAGTAAAGTGCTACCATTTTTCTTACAGACATTTGGCGATATTTTAAAC  >  CP000730/2704666‑2704940

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