Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,840,050 T→C S390G (AGT→GGT)  topB ← DNA topoisomerase III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,840,0500TC84.2% 41.5 / 1.9 19S390G (AGT→GGT) topBDNA topoisomerase III
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (8/8);  total (10/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.51e-01

ATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGC  >  NZ_CP009273/1839978‑1840138
                                                                        |                                                                                        
aTTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACCCCGTGCGGTCGGAatga                                                                         <  1:377342/90‑1 (MQ=255)
   gCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACCCCGTGCGGTCGGAATGATgg                                                                      >  2:316543/1‑90 (MQ=255)
                   ggCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACCCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGAc                                                      >  2:489123/1‑90 (MQ=255)
                            gTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACCCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTc                                             <  1:447873/90‑1 (MQ=255)
                               ataGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACCCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATc                                          <  1:489123/90‑1 (MQ=255)
                                     cttcgcttcgTTCTCCGTCAGGTTGATCGCAGAACCCCGTGCGGTCGGa                                                                             >  1:477472/1‑49 (MQ=255)
                                     cttcgcttcgTTCTCCGTCAGGTTGATCGCAGAACCCCGTGCGGTCGGa                                                                             <  2:477472/49‑1 (MQ=255)
                                     cttcgcttcgTTCTCCGTCAGGTTGATCGCAGAACCCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACa                                    >  1:454852/1‑90 (MQ=255)
                                     cttcgcttcgTTCTCCGTCAGGTTGATCGCAGAACCCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACa                                    >  1:216555/1‑90 (MQ=255)
                                       tcgcttcgTTCTCCGTCAGGTTGATCGCAGAACCCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGc                                  >  1:261393/1‑90 (MQ=255)
                                          cttcgTTCTCCGTCAGGTTGATCGCAGAACCCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGt                               <  1:477341/90‑1 (MQ=255)
                                           ttcgTTCTCCGTCAGGTTGATCGCAGAACCCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGtt                              <  2:261393/90‑1 (MQ=255)
                                                 ctcCGTCAGGTTGATCGCAGAACCCCGTGCGGTCGGAATGATGGCgt                                                                   <  1:354885/47‑1 (MQ=255)
                                                 ctcCGTCAGGTTGATCGCAGAACCCCGTGCGGTCGGAATGATGGCgt                                                                   >  2:354885/1‑47 (MQ=255)
                                                        aGGTTGATCGCAGAACCCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGa                 <  2:454852/90‑1 (MQ=255)
                                                         ggTTGATCGCAGAACCCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGAt                >  2:34649/1‑90 (MQ=255)
                                                             gATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATc               >  1:380055/1‑87 (MQ=255)
                                                             gATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATc               <  2:380055/87‑1 (MQ=255)
                                                                       cTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGc  >  2:182825/1‑90 (MQ=255)
                                                                        |                                                                                        
ATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGC  >  NZ_CP009273/1839978‑1840138

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