Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,792,826 C→T *286* (TAG→TAA yibB ← YibB family protein, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,792,8260CT98.0% 148.7 / ‑4.2 49*286* (TAG→TAAyibBYibB family protein, function unknown
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (27/21);  total (27/22)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

GTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATA  >  NC_000913/3792698‑3792960
                                                                                                                                |                                                                                                                                      
gTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTg                                                                                                                              >  1:551978/1‑139 (MQ=255)
gTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTg                                                                                                                              >  2:74229/1‑139 (MQ=255)
                           attgattTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATtctc                                                                                                                >  2:109296/1‑126 (MQ=255)
                           attgattTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATtctc                                                                                                                <  1:109296/126‑1 (MQ=255)
                                         aTTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAgtg                                                                                     >  1:1493539/1‑139 (MQ=255)
                                           tGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGtt                                                                                   <  1:1090584/139‑1 (MQ=255)
                                                     aGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAgtg                                                                                     >  2:657861/1‑127 (MQ=255)
                                                     aGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAgtg                                                                                     <  1:657861/127‑1 (MQ=255)
                                                                 gCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACaaa                                                             >  2:214151/1‑139 (MQ=255)
                                                                     aaGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTCTTACTCCTGAAGCAGCGAAACAAATCGa                                                         >  1:920646/1‑139 (MQ=255)
                                                                                  ccTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAgtgt                                                                                    <  1:904153/99‑1 (MQ=255)
                                                                                  ccTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAgtgt                                                                                    >  2:904153/1‑99 (MQ=255)
                                                                                  ccTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCctct                                            >  1:714704/1‑139 (MQ=255)
                                                                                     aTCTCAGTTGGTATCAAAACGCTACCCACTATTGGAGGGAAATTTATTTTCTGAAAAGGAAAATTCTCAGCCCCGGATTTTTTCCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCAAACCATTTTCCTCTccc                                         >  1:880695/1‑139 (MQ=255)
                                                                                            ttGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAgc                                                                                                             <  1:1241753/64‑1 (MQ=255)
                                                                                            ttGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAgc                                                                                                             >  2:1241753/1‑64 (MQ=255)
                                                                                             tGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGt                                 <  1:214151/139‑1 (MQ=255)
                                                                                             tGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTACTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGt                                 <  2:920646/139‑1 (MQ=255)
                                                                                                tATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCAttttt                                                                                                 <  1:1006433/72‑1 (MQ=255)
                                                                                                tATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCAttttt                                                                                                 >  2:1006433/1‑72 (MQ=255)
                                                                                                      aaCGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGa                        >  2:537033/1‑139 (MQ=255)
                                                                                                          cTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAag                                                                         >  1:296475/1‑86 (MQ=255)
                                                                                                          cTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAag                                                                         <  2:296475/86‑1 (MQ=255)
                                                                                                            aCCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTg                                                                                                        >  2:1381171/1‑53 (MQ=255)
                                                                                                            aCCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTg                                                                                                        <  1:1381171/53‑1 (MQ=255)
                                                                                                            aCCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGcccc                                                                                            <  1:896766/65‑1 (MQ=255)
                                                                                                            aCCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGcccc                                                                                            >  2:896766/1‑65 (MQ=255)
                                                                                                               cACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCgg               >  2:944778/1‑139 (MQ=255)
                                                                                                                aCTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCggg              <  2:714704/139‑1 (MQ=255)
                                                                                                                  tAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCAtttt                                                                                                  <  1:1337180/53‑1 (MQ=255)
                                                                                                                  tAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCAtttt                                                                                                  >  2:1337180/1‑53 (MQ=255)
                                                                                                                  tAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGc            <  1:348819/139‑1 (MQ=255)
                                                                                                                     tGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGtt         <  2:1493539/139‑1 (MQ=255)
                                                                                                                       gAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTg       >  1:10692/1‑139 (MQ=255)
                                                                                                                       gAGGGAAACTTAATTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGaa                                                        <  1:249793/90‑1 (MQ=255)
                                                                                                                       gAGGGAAACTTAATTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGaa                                                        >  2:249793/1‑90 (MQ=255)
                                                                                                                         gggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAgtgt                                                                                    <  1:1033830/60‑1 (MQ=255)
                                                                                                                         gggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAgtgt                                                                                    >  2:1033830/1‑60 (MQ=255)
                                                                                                                         gggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTc                                                >  2:1519208/1‑96 (MQ=255)
                                                                                                                         gggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTc                                                <  1:1519208/96‑1 (MQ=255)
                                                                                                                           gAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTa                                                                                          >  2:1477585/1‑52 (MQ=255)
                                                                                                                           gAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTa                                                                                          <  1:1477585/52‑1 (MQ=255)
                                                                                                                           gAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCctct                                            <  1:613805/98‑1 (MQ=255)
                                                                                                                           gAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGt                                 <  2:363445/109‑1 (MQ=255)
                                                                                                                           gAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGt                                 >  1:363445/1‑109 (MQ=255)
                                                                                                                           gAAACTTATTTTCTTAAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACTATTTTCctct                                            >  2:613805/1‑98 (MQ=255)
                                                                                                                            aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCctc                                             <  1:1070228/96‑1 (MQ=255)
                                                                                                                            aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCctc                                             >  2:1070228/1‑96 (MQ=255)
                                                                                                                            aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATa  >  1:1053907/1‑139 (MQ=255)
                                                                                                                                |                                                                                                                                      
GTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATA  >  NC_000913/3792698‑3792960

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: