Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,509,958 Δ1 bp intergenic (‑195/+314) yhfZ ← / ← trpS hypothetical protein/tryptophan tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,509,9560T.98.6% 310.0 / ‑3.3 70intergenic (‑193/+316)yhfZ/trpShypothetical protein/tryptophan tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (14/55);  minor base A (1/0);  total (15/55)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAA  >  NC_000913/3509881‑3510023
                                                                           |                                                                   
cTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCct                                                                     >  2:3045396/1‑76 (MQ=255)
  cATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCcttc                                                                   <  2:772010/76‑2 (MQ=255)
      tATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCtttt                                                              <  2:2764360/76‑1 (MQ=255)
        tGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCtttttt                                                            <  1:5054252/76‑1 (MQ=255)
         ggTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTg                                                           <  1:3589416/76‑1 (MQ=255)
          gTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGa                                                          >  2:4891903/1‑76 (MQ=255)
          gTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGa                                                          >  1:5278612/1‑76 (MQ=255)
           ttttCAAGCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGAt                                                         <  2:5458196/76‑1 (MQ=255)
            tttCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATa                                                        <  1:335441/76‑1 (MQ=255)
            tttCAATCGATTACCGCAATTTTCAACCAACAAATAACGGGCCAGCGTCGTGATGACGACGCATTCTTGTTTGATa                                                        >  1:4496207/1‑76 (MQ=255)
             ttCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATaa                                                       <  1:3514960/76‑1 (MQ=255)
             ttCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATaa                                                       >  2:3956156/1‑76 (MQ=255)
              tCAATCGATTACCGCGCTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATaaa                                                      <  2:3374242/76‑1 (MQ=255)
              tCAATCGATTACCGCAGTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATaaa                                                      <  2:5099834/76‑1 (MQ=255)
              tCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATaaa                                                      <  2:4632780/76‑1 (MQ=255)
               cAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAAc                                                     <  2:4547560/76‑1 (MQ=255)
               cAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAAc                                                     <  1:1726544/76‑1 (MQ=255)
                aaTCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACg                                                    <  1:3702672/76‑1 (MQ=255)
                aaTCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACg                                                    <  1:3991962/76‑1 (MQ=255)
                 aTCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGt                                                   <  2:5176224/76‑1 (MQ=255)
                  tCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGtt                                                  <  1:5405125/76‑1 (MQ=255)
                   cGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTc                                                 <  1:4296664/76‑1 (MQ=255)
                   cGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTc                                                 <  1:5584524/76‑1 (MQ=255)
                     aTTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGa                                               <  2:44708/76‑1 (MQ=255)
                      ttACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATACACGTTCGAc                                              >  2:5417815/1‑76 (MQ=255)
                      ttACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGAc                                              <  1:3600680/76‑1 (MQ=255)
                      ttACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGAc                                              <  2:4502783/76‑1 (MQ=255)
                       tACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACg                                             >  2:5005608/1‑76 (MQ=255)
                       tACCGCACTTTGCACCCAACAAATAACTGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACg                                             <  1:3248862/76‑1 (MQ=255)
                        tccGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGc                                            <  1:5682238/75‑1 (MQ=255)
                          cGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCat                                          <  1:4143239/76‑1 (MQ=255)
                          cGCACGTTTCCCCCAACAAACAACGGGCCAGCGGCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGCCGCat                                          <  2:4858350/76‑1 (MQ=255)
                              cTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCataata                                      <  2:4806870/76‑1 (MQ=255)
                              cTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCataata                                      >  2:4052777/1‑76 (MQ=255)
                              cTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCataata                                      <  2:3787642/76‑1 (MQ=255)
                               ttttCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCAtaataa                                     <  2:5012032/76‑1 (MQ=255)
                                  tCACCCAACAAATAACGGTCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAAtt                                  <  1:3336196/76‑1 (MQ=255)
                                  tCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAAtt                                  >  1:2173935/1‑76 (MQ=255)
                                      ccAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCa                              <  2:401436/76‑1 (MQ=255)
                                       caacaaATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAg                             <  1:3125369/76‑1 (MQ=255)
                                         acaaATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTa                           <  1:3918895/76‑1 (MQ=255)
                                          caaATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTaa                          <  2:469531/76‑1 (MQ=255)
                                          caaATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTaa                          <  1:2472476/76‑1 (MQ=255)
                                          caaATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTaa                          <  1:3056712/76‑1 (MQ=255)
                                           aaaTAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTaaa                         <  2:2097747/76‑1 (MQ=255)
                                            aaTAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAAt                        <  2:4477269/76‑1 (MQ=255)
                                            aaTAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGATCGACGCATAATAAATTTTCAGTAAAt                        <  2:2943591/76‑1 (MQ=255)
                                             aTAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAAtt                       >  1:3623376/1‑76 (MQ=255)
                                              tAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTg                      >  2:3920674/1‑76 (MQ=255)
                                               aaCGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGa                     <  1:2295198/76‑1 (MQ=255)
                                                 cGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAat                   >  1:5098441/1‑76 (MQ=255)
                                                 cGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAat                   <  1:5417085/76‑1 (MQ=255)
                                                 cGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAat                   >  1:3532045/1‑76 (MQ=255)
                                                  gggCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAANNAATTTTCAGTAAATTGAata                  <  1:3523881/76‑1 (MQ=255)
                                                    gCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATAtt                <  1:968482/76‑1 (MQ=255)
                                                      cAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTta              <  2:4289571/76‑1 (MQ=255)
                                                       aGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCTACGCATAATAAATTTTCAGTAAATTGAATATTtat             <  1:4731756/76‑1 (MQ=255)
                                                       aGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTtat             <  2:1741983/76‑1 (MQ=255)
                                                       aGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTtat             <  1:846318/76‑1 (MQ=255)
                                                            cgtGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCtt        <  2:1772535/76‑1 (MQ=255)
                                                            cgtGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCtt        <  2:2220136/76‑1 (MQ=255)
                                                             gtGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTc       >  1:1738320/1‑76 (MQ=255)
                                                             gtGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTc       <  2:228312/76‑1 (MQ=255)
                                                              tgatgaCGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCa      <  2:2173935/76‑1 (MQ=255)
                                                              tgatgaCGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCa      >  1:918014/1‑76 (MQ=255)
                                                               gatgaCGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAg     <  2:5278611/76‑1 (MQ=255)
                                                               gatgaCGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAg     <  1:105481/76‑1 (MQ=255)
                                                               gatgaCGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAg     >  2:3290226/1‑76 (MQ=255)
                                                               gatgaCGACACCTGCTTTGTGGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATGGATATCTTcct     <  2:5321604/76‑3 (MQ=255)
                                                                atgaCGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAgg    <  2:2888727/76‑1 (MQ=255)
                                                                  gacgacGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGaa  <  2:4496206/76‑1 (MQ=255)
                                                                  gacgacGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGaa  <  2:3058811/76‑1 (MQ=255)
                                                                           |                                                                   
CTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAA  >  NC_000913/3509881‑3510023

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