Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 231,042 A→C G76G (GGT→GGG gloB ← hydroxyacylglutathione hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273231,0420AC97.9% 149.8 / ‑5.6 47G76G (GGT→GGGgloBhydroxyacylglutathione hydrolase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (18/28);  total (19/28)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAATATGG  >  NZ_CP009273/230956‑231127
                                                                                      |                                                                                     
aTTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGCCCa                                                                                    <  2:589437/90‑1 (MQ=255)
  tACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGCCCATa                                                                                  >  2:216547/1‑90 (MQ=255)
   acacTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGCCCATac                                                                                 <  1:1557194/90‑1 (MQ=255)
             cATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAAttt                                                                       >  1:25824/1‑90 (MQ=255)
              aTGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTtg                                                                      >  1:1252714/1‑90 (MQ=255)
              aTGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTtg                                                                      >  1:872184/1‑90 (MQ=255)
              aTGCCCCAAAACGAAGACAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTtg                                                                      >  1:1187243/1‑90 (MQ=255)
                  cccAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGa                                                                  >  2:149634/1‑90 (MQ=255)
                   ccAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGaa                                                                 <  1:199305/90‑1 (MQ=255)
                    cAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGaaa                                                                <  1:839203/90‑1 (MQ=255)
                    cAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGaaa                                                                <  2:1169276/90‑1 (MQ=255)
                    cAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGaaa                                                                >  2:1152461/1‑90 (MQ=255)
                     aaaaCGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGAAAc                                                               >  2:923135/1‑90 (MQ=255)
                           aaGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGAAACTTTTcc                                                         >  1:1580186/1‑90 (MQ=255)
                           aaGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGAAACTTTTcc                                                         >  2:795093/1‑90 (MQ=255)
                            aGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGAAACTTTTcca                                                        >  1:80158/1‑90 (MQ=255)
                              gCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGAAACTTTTccacc                                                      <  1:1213689/90‑1 (MQ=255)
                                         aTCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGAAACTTTTCCCCCAGTTCTTTTAc                                           <  1:1350687/90‑1 (MQ=255)
                                           cTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTAcgc                                         <  2:1456265/90‑1 (MQ=255)
                                             ttGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTAcgccg                                       <  1:85412/90‑1 (MQ=255)
                                             ttGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTAcgccg                                       <  1:449681/90‑1 (MQ=255)
                                                aCTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCa                                    <  1:1440792/90‑1 (MQ=255)
                                                  tACCTGTGTTGTTACCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAAc                                  <  2:918727/90‑1 (MQ=255)
                                                   aCCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGccc                                                                                     <  1:580866/38‑1 (MQ=255)
                                                   aCCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACg                                 >  1:727303/1‑90 (MQ=255)
                                                   aCCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACg                                 >  2:702934/1‑90 (MQ=255)
                                                    ccTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGt                                <  1:180441/90‑1 (MQ=255)
                                                        tgttgtTCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATc                            <  2:772661/90‑1 (MQ=255)
                                                        tgttgtTCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATc                            <  2:312763/90‑1 (MQ=255)
                                                        tgttgtTCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATc                            <  2:1520382/90‑1 (MQ=255)
                                                         gttgttCCCTTATCTTGTGTCTCTTGTGGCCCATAcaccac                                                                            >  1:1128799/1‑41 (MQ=255)
                                                         gttgttCCCTTATCTTGTGTCTCTTGTGGCCCATAcaccac                                                                            <  2:1128799/41‑1 (MQ=255)
                                                          ttgttCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCgt                          <  1:1013181/90‑1 (MQ=255)
                                                          ttgttCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCgt                          <  2:1214625/90‑1 (MQ=255)
                                                           tgttCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATc                            >  1:706687/1‑87 (MQ=255)
                                                           tgttCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATc                            <  2:706687/87‑1 (MQ=255)
                                                           tgttCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCg                           >  2:466534/1‑88 (MQ=255)
                                                           tgttCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCg                           <  1:466534/88‑1 (MQ=255)
                                                            gttCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCgtgg                        <  1:865346/90‑1 (MQ=255)
                                                            gttCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCgtgg                        <  2:1574301/90‑1 (MQ=255)
                                                            gggcccTTATCCTGTGGCCCTTGGGGCCCAAACACCCCAAATTATGGAAACCTTTCCACCAGGTCTTTTACGCCGCCAACGTGATCgtgg                        <  2:916561/87‑1 (MQ=255)
                                                              tCCCTTATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCgtggtg                      <  2:214509/90‑1 (MQ=255)
                                                                   tATCTTGTGTCTCTTGTGGCCCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGAtggt                 <  1:1414075/90‑1 (MQ=255)
                                                                        tgtgTCTCTTGTGGCCCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATggtgggtg            <  1:996552/90‑1 (MQ=255)
                                                                         gtgtCTCTTGTGGCCCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTga           >  1:855422/1‑90 (MQ=255)
                                                                         gtgtCTCTTGTGGCCCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTga           >  1:581512/1‑90 (MQ=255)
                                                                                  gtgGCCCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAATATgg  >  2:994323/1‑90 (MQ=255)
                                                                                   tgGCCCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTAcgc                                         <  2:829023/50‑1 (MQ=255)
                                                                                      |                                                                                     
ATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAATATGG  >  NZ_CP009273/230956‑231127

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