Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,643,049 G→A 25.4% A341V (GCT→GTT)  accC ← biotin carboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,643,0490GA25.4% 89.7 / 28.4 51A341V (GCT→GTT) accCbiotin carboxylase
Reads supporting (aligned to +/- strand):  ref base G (17/21);  new base A (4/9);  total (21/30)
Fisher's exact test for biased strand distribution p-value = 5.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCTACCATCGAATCATAATATGGCGGTATCGTATAATTAGTATAACATGCTGACTCTATTCGAACACCATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACGGTAACACGTCACCCATAGCAACTTGTAATTGTAATTTAACTAAATCAATTCCTGTTACCATTTCAGTTACAGGA  >  CP000730/1642909‑1643180
                                                                                                                                            |                                                                                                                                   
tCGCTACCATCGAATCATAATATGGCTGTATCGTATAATTAGTATAACATGCTGACTCTATTCGAACACCATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAg                                                                                                                                     <  2:368072/141‑1 (MQ=255)
                atattaTGGCGGTATCGTATAATTAGTATAACATGCTGACTCTATTCGAACACCATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAACATTAATTCTAAATTc                                                                                                                     <  1:42027/141‑1 (MQ=255)
                 taataTGGCGGTATCGTATAATTAGTATAACATGCTGACTCTATTCGAACACCATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCa                                                                                                                    <  2:385700/141‑1 (MQ=255)
                  aataTGGCGGTATCGTATAATTAGTATAACATGCTGACTCTATTCGAACACCATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCaa                                                                                                                   <  1:25812/141‑1 (MQ=255)
                   atatGGCGGTATCGTATAATTAGTATAACATGCTGACTCTATTCGAACACCATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAAt                                                                                                                  <  2:182476/141‑1 (MQ=255)
                    tatGGCGGTATCGTATAATTAGTATAACATGCTGACTCTATTCGAACACCATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAAtt                                                                                                                 <  1:231308/141‑1 (MQ=255)
                    tatGGCGGTATCGTATAATTAGTATAACATGCTGACTCTATTCGAACACCATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAACATTAATTCTAAATTCAAtt                                                                                                                 <  1:35180/141‑1 (MQ=255)
                      tGGCGGTATCGTATAATTAGTATAACATGCTGACTCTATTCGAACACCATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGc                                                                                                               <  2:159621/141‑1 (MQ=255)
                      tGGCGGTATCGTATAATTAGTATAACATGCTGACTCTATTCGAACACCATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGc                                                                                                               <  2:157427/141‑1 (MQ=255)
                         cGGTATCGTATAATTAGTATAACATGCTGACTCTATTCAAACACCATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCgtg                                                                                                            <  2:499665/141‑1 (MQ=255)
                          ggTATCGTATAATTAGTATAACATGCTGACTCTATTCGAACACCATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCgtgt                                                                                                           >  1:336443/1‑141 (MQ=255)
                           gTATCGTATAGTTAGTATAACATGCTGACTCTATTCGAACACCATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTc                                                                                                          >  1:367761/1‑141 (MQ=255)
                              tCGTATAATTAGTATAACATGCTGACTCTATTCGAACACCATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAACATTAATTCTAAATTCAATTGCGTGTCCTg                                                                                                       <  1:183094/141‑1 (MQ=255)
                               cGTATAATTAGTATAACATGCTGACTCTATTCGAACACCATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCCTGt                                                                                                      <  2:471246/141‑1 (MQ=255)
                               cGTATAATTAGTATAACATGCTGACTCTATTCGAACACCATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCCTGt                                                                                                      <  2:471208/141‑1 (MQ=255)
                                    aaTTAGTATAACATGCTGACTCTATTCGAACACCATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCCTGttaatt                                                                                                 >  2:101352/1‑141 (MQ=255)
                                         gTATAACATGCTGACTCTATTCGAACACCATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAata                                                                                            >  1:289034/1‑141 (MQ=255)
                                             aaCATGCTGACTCTATTCGAACACCATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATAtctt                                                                                        >  1:287798/1‑141 (MQ=255)
                                                 tGCTGACTCTATTCGAACACCATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAACATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTg                                                                                    >  1:429712/1‑141 (MQ=255)
                                                       ctctATTCGAACACCATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTtata                                                                              <  2:171336/141‑1 (MQ=255)
                                                       ctctATTCGAACAACATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTtata                                                                              <  1:309774/141‑1 (MQ=255)
                                                         ctATTCGAACACCATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACg                                                                            >  1:150863/1‑141 (MQ=255)
                                                               gAACACCATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATAGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACGGTAaca                                                                      <  2:211807/141‑1 (MQ=255)
                                                                 acacCATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATAtc                                                                                          <  1:200913/119‑1 (MQ=255)
                                                                 acacCATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATAtc                                                                                          >  2:200913/1‑119 (MQ=255)
                                                                 acacCATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAACATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACGGTAACACg                                                                    <  1:149776/141‑1 (MQ=255)
                                                                  cacCATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACGGTAACACGt                                                                   <  2:429407/141‑1 (MQ=255)
                                                                   acCATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACGGTAACACGTc                                                                  >  2:249964/1‑141 (MQ=255)
                                                                      atatCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAACATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACGGTAACACGTCAcc                                                               <  2:147074/141‑1 (MQ=255)
                                                                        atCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACGGTAACACGTCACCCa                                                             >  1:401065/1‑141 (MQ=255)
                                                                        atCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAACATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACGGTAACACGTCACCCa                                                             <  2:460176/141‑1 (MQ=255)
                                                                              cTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACGGTAACACGTCACCCATAGCaa                                                       >  2:36157/1‑141 (MQ=255)
                                                                              cTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAACATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACGGTAACACGTCACCCATAGCaa                                                       >  1:276994/1‑141 (MQ=255)
                                                                                    cAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCTTGTTAATTTAATATCTTCTTGTTTATACGGTAACACGTCACCCATAGCAACttgta                                                 <  2:494927/141‑1 (MQ=255)
                                                                                      aGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAACATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACGGTAACACGTCACCCATAGCAACttgtaat                                               >  1:327640/1‑141 (MQ=255)
                                                                                                   ttttACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAACATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACGGTAACACGTCACCCATAGCAACTTGTAATTGTAATTTAACTa                                  <  1:423384/141‑1 (MQ=255)
                                                                                                        ccTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACGGTAACACGTCACCCATAGCAACTTGTAATTGTAATTTAACTAAATCa                             >  1:404283/1‑141 (MQ=255)
                                                                                                           ggTGATGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACGGTAACACGTCACCCATAGCAACTTGTAATTGTAATTTAACTAAATCAAtt                          >  1:46758/1‑141 (MQ=255)
                                                                                                               aTGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACGGTAACACGTCACCCATAGCAACTTGTAATTGTAATTTAACTAAATCAATTCCTg                      <  2:289034/141‑1 (MQ=255)
                                                                                                               aTGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACGGTAACACGTCACCCATAGCAACTTGTAATTGTAATTTAACTAAATCAATTCCTg                      <  1:372405/141‑1 (MQ=255)
                                                                                                                 ggCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACGGTAACACGTCACCCATAGCAACTTGTAATTGTAATTTAACTAAATCAATTCCTGtt                    <  2:386977/141‑1 (MQ=255)
                                                                                                                  gCATAAAGTTCTTGTAAGGATTTTCAACATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACGGTAACACGTCACCCATAGCAACTTATAATTGTAATTTAACTAAATCAATTTCttttt                   >  2:132945/1‑140 (MQ=255)
                                                                                                                    aTAAAGTTCTTGTAAGGATTTTCAACATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACGGTAACACGTCACCCATAGCAACTTGTAATTGTAATTTAACTAAATCAATTCCTGTTAcc                 <  2:400121/141‑1 (MQ=255)
                                                                                                                       aaGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACGGTAACACGTCACCCATAGCAACTTGTAATTGTAATTTAACTAAATCAATTCCTGTTACCAtt              <  1:445717/141‑1 (MQ=255)
                                                                                                                        aGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACGGTAACACGTCACCCATAGCAACTTGTAATTGTAATTTAACTAAATCAATTCCTGTTACCAttt             >  2:199403/1‑141 (MQ=255)
                                                                                                                        aGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACGGTAACACGTCACCCATAGCAACTTGTAATTGTAATTTAACTAAATCAATTCCTGTTACCAttt             >  2:199402/1‑141 (MQ=255)
                                                                                                                          ttCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACGGTAACACGTCACCCATAGCAACTTGTAATTGTAATTTAACTAAATCAATTCCTGTTACCATTTCa           <  2:43925/141‑1 (MQ=255)
                                                                                                                            cTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACGGTAACACGTCACCCATAGCAACTTGTAATTGTAATTTAACTAAATCAATTCCTGTTACCATTTCAGt         <  2:392695/141‑1 (MQ=255)
                                                                                                                             ttGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACGGTAACACGTCACCCATAGCAACTTGTAATTGTAATTTAACTAAATCAATTCCTGTTACCATTTCAGtt        >  1:245874/1‑141 (MQ=255)
                                                                                                                                 aaGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACGGTAACACGTCACCCATAGCAACTTGTAATTGTAATTTAACTAAATCAATTCCTGTTACCATTTCAGTTACAg    >  1:502118/1‑141 (MQ=255)
                                                                                                                                 aaGGATTTTCAACATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACGGTAACACGTCACCCATAGCAACTTGTAATTGTAATTTAACTAAATCAATTCCTGTTACCATTTCAGTTACAg    <  1:480222/141‑1 (MQ=255)
                                                                                                                                   ggATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACGGTAACACGTCACCCATAGCAACTTGTAATTGTAATTTAACTAAATCAATTCCTGTTACCATTTCAGTTACAGGa  >  2:167779/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                   
TCGCTACCATCGAATCATAATATGGCGGTATCGTATAATTAGTATAACATGCTGACTCTATTCGAACACCATATCCACCTGGTGCAAGATATTGCTCAATTTTACCTGGTGATGGCATAAAGTTCTTGTAAGGATTTTCAGCATTAATTCTAAATTCAATTGCGTGTCCTGTTAATTTAATATCTTCTTGTTTATACGGTAACACGTCACCCATAGCAACTTGTAATTGTAATTTAACTAAATCAATTCCTGTTACCATTTCAGTTACAGGA  >  CP000730/1642909‑1643180

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