Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 986,671 A→G R104R (CGT→CGC)  ompF ← outer membrane porin 1a (Ia,b,F)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913986,6710AG98.4% 218.9 / NA 62R104R (CGT→CGC) ompFouter membrane porin 1a (Ia,b,F)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (32/29);  minor base T (0/1);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGACGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATG  >  NC_000913/986554‑986806
                                                                                                                     |                                                                                                                                       
aTCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGa                                                                                                                    >  2:816812/1‑139 (MQ=255)
  caccacCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGt                                                                                                                                    >  2:66934/1‑121 (MQ=255)
  caccacCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGt                                                                                                                                    <  1:66934/121‑1 (MQ=255)
  caccacCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGc                                                                                                                  >  2:133454/1‑139 (MQ=255)
  caccacCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGc                                                                                                                  >  1:62898/1‑139 (MQ=255)
       ccaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAg                                                                                                             >  1:1124451/1‑139 (MQ=255)
                     caTATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAgttg                                                                                               <  1:674338/139‑1 (MQ=255)
                            gtgtAGCCCAGTGCATCATTAACCGCACCGTAGTTACGGCCGTAATCGAAAGTACCAACGTCAGCGTATTTAACACACGCGGATGCCAGTCGCGTTTTGTTACCCGTTTGATCGTCAGCGCCTTCAGAGTTGCTACCCt                                                                                        <  2:393342/139‑1 (MQ=255)
                                      gTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTtata                                                                              >  1:779034/1‑139 (MQ=255)
                                       tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGt                                                                                  >  1:752246/1‑134 (MQ=255)
                                       tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGt                                                                                  <  2:752246/134‑1 (MQ=255)
                                       tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTatat                                                                             >  1:680854/1‑139 (MQ=255)
                                       tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTatat                                                                             >  1:359062/1‑139 (MQ=255)
                                       tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTatat                                                                             >  2:1002830/1‑139 (MQ=255)
                                       tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTatat                                                                             >  1:154749/1‑139 (MQ=255)
                                         catcatAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTg                                                                                                                     <  2:108622/97‑1 (MQ=255)
                                         catcatAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTg                                                                                                                     >  1:108622/1‑97 (MQ=255)
                                            catAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                        <  1:816812/139‑1 (MQ=255)
                                            catAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                        <  2:1116239/139‑1 (MQ=255)
                                            catAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                        <  2:1124451/139‑1 (MQ=255)
                                            catAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                        <  2:978141/139‑1 (MQ=255)
                                              tAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACt                                                                      <  2:779034/139‑1 (MQ=255)
                                                 aCCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAgttg                                                                                               >  1:843780/1‑111 (MQ=255)
                                                 aCCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAgttg                                                                                               <  2:843780/111‑1 (MQ=255)
                                                 aCCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGAc                                                                   <  2:154749/139‑1 (MQ=255)
                                                 aCCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGAc                                                                   <  1:1002830/139‑1 (MQ=255)
                                                      acCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATaa                                                              >  1:121265/1‑139 (MQ=255)
                                                                  gCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCg                                                                                                                 <  2:24346/76‑1 (MQ=255)
                                                                  gCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCg                                                                                                                 >  1:24346/1‑76 (MQ=255)
                                                                  gCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGt                                                                                  >  2:147681/1‑107 (MQ=255)
                                                                  gCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGt                                                                                  <  1:147681/107‑1 (MQ=255)
                                                                  gCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGAc                                                                   >  1:782852/1‑122 (MQ=255)
                                                                  gCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGAc                                                                   <  2:782852/122‑1 (MQ=255)
                                                                    cGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGt                                                         >  1:719238/1‑130 (MQ=255)
                                                                    cGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGt                                                         <  2:719238/130‑1 (MQ=255)
                                                                    cGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGa                                                >  2:1323979/1‑139 (MQ=255)
                                                                      tAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAgttgt                                                                                              >  2:1167056/1‑91 (MQ=255)
                                                                      tAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAgttgt                                                                                              <  1:1167056/91‑1 (MQ=255)
                                                                          cGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGAt                                          <  2:668799/139‑1 (MQ=255)
                                                                                   cAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTc                                 <  2:62898/139‑1 (MQ=255)
                                                                                          aGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTAc                                                                                                                            <  2:155621/41‑1 (MQ=255)
                                                                                          aGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTAc                                                                                                                            >  1:155621/1‑41 (MQ=255)
                                                                                          aGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTa                          >  1:554782/1‑139 (MQ=255)
                                                                                                tttAAGACCTGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGAt                                                    <  2:807502/107‑1 (MQ=255)
                                                                                                tttAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGAt                                                    >  1:807502/1‑107 (MQ=255)
                                                                                                   aaGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGa                 <  2:554782/139‑1 (MQ=255)
                                                                                                   aaGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGa                 <  1:133454/139‑1 (MQ=255)
                                                                                                        ccGCGAATGCCAGGGGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCAc                                                                       <  2:917035/80‑1 (MQ=255)
                                                                                                        ccGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTTGAAGTTATATTCCCAc                                                                       >  1:917035/1‑80 (MQ=255)
                                                                                                          gcgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTc                                                                                                               <  2:1102220/38‑1 (MQ=255)
                                                                                                          gcgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTc                                                                                                               >  1:1102220/1‑38 (MQ=255)
                                                                                                          gcgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCAt          >  2:957608/1‑139 (MQ=255)
                                                                                                          gcgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCAt          >  2:373693/1‑139 (MQ=255)
                                                                                                           cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCt                                                                                        <  1:958808/60‑1 (MQ=255)
                                                                                                           cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCt                                                                                        >  2:958808/1‑60 (MQ=255)
                                                                                                           cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGt                                                                                  >  2:34489/1‑66 (MQ=255)
                                                                                                           cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGt                                                                                  <  1:34489/66‑1 (MQ=255)
                                                                                                           cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGt                                                         <  2:679244/91‑1 (MQ=255)
                                                                                                           cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGt                                                         >  1:679244/1‑91 (MQ=255)
                                                                                                               tGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGa                 <  2:1016372/127‑1 (MQ=255)
                                                                                                               tGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGa                 >  1:1016372/1‑127 (MQ=255)
                                                                                                                  cAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATg  >  1:579180/1‑139 (MQ=255)
                                                                                                                     |                                                                                                                                       
ATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGACGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATG  >  NC_000913/986554‑986806

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