Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 987,197 A→T 8.7% intergenic (‑215/+388) ompF ← / ← asnS outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913987,1970AT 8.7% 131.8 / 5.5 46intergenic (‑215/+388)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (20/22);  new base T (2/2);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.61e-01

TTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATC  >  NC_000913/987051‑987344
                                                                                                                                                  |                                                                                                                                                   
tttACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAg                                                                                                                                                   >  2:858945/1‑149 (MQ=255)
     ggAACTATTGATGGGGGTTTGGTGTCGTTATGTGTCTGGAGGCATCTTTCCATTCAGACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACt                                                                                                                                              <  2:12598/149‑1 (MQ=255)
     ggAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACt                                                                                                                                              >  2:278177/1‑149 (MQ=255)
       aaCTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGGACttt                                                                                                                                            >  1:828338/1‑149 (MQ=255)
       aaCTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACttt                                                                                                                                            <  1:547744/149‑1 (MQ=255)
         cTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTtgt                                                                                                                                          >  1:238822/1‑149 (MQ=255)
          tATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAACAGAAAAAATATGTAACCAAAAGTAAAAATTAAGGGACtgtgtg                                                                                                                                         >  1:268404/1‑149 (MQ=255)
           aTTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGa                                                                                                                                        <  1:382980/149‑1 (MQ=255)
            atgatgaAAATTTGGTGTCTTTATGTGTATACAGGCATCTTTCCAACCAAACTAACCATTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGAAAAAAAAAAAAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGaa                                                                                                                                       <  2:426746/148‑1 (MQ=255)
                   gagTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTc                                                                                                                                <  1:693170/149‑1 (MQ=255)
                      tttCGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAACGTAAAATTTAAGGATCTTTGTGAACACCGTCata                                                                                                                             <  2:618728/149‑1 (MQ=255)
                         ggTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATCttt                                                                                                                          >  1:206508/1‑149 (MQ=255)
                             tCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCAt                                                                                                                      <  1:738108/149‑1 (MQ=255)
                               tttATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAATTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCTTag                                                                                                                    >  1:17576/1‑149 (MQ=255)
                               tttATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATag                                                                                                                    >  1:854737/1‑149 (MQ=255)
                               tttATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCAATag                                                                                                                    >  1:778054/1‑149 (MQ=255)
                                         tGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCCTTGAGACGtgata                                                                                                          >  2:168833/1‑148 (MQ=255)
                                         tGCAGGCATCTTTCCATTCAAACTAACGGTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCTTAGAGACGtgatg                                                                                                          >  1:505640/1‑149 (MQ=255)
                                         tGCAGGCATCTTTCCATTCAAAATAACGTTACGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGtgatg                                                                                                          <  2:10533/149‑1 (MQ=255)
                                                aTCTTTCCATTCAAAATAACGTTTCGCTACCGTGAAAGTGCTACAAAGCTAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTa                                                                                                   <  1:97075/149‑1 (MQ=255)
                                                  ctctCCATTCAAACTACCGTTGCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACa                                                                                                 <  2:455730/146‑1 (MQ=255)
                                                                   cGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTTAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTAttt                                                                                >  2:836178/1‑149 (MQ=255)
                                                                   cGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTTAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTAttt                                                                                >  2:508274/1‑149 (MQ=255)
                                                                     tttCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTtat                                                                              >  1:321815/1‑149 (MQ=255)
                                                                           tACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACGTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGAtt                                                                        <  1:775740/149‑1 (MQ=255)
                                                                           tACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTACAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGAtt                                                                        >  2:568565/1‑149 (MQ=255)
                                                                                tGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAGGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAAGCTATTTATATGATTTCCtt                                                                   >  1:748888/1‑149 (MQ=255)
                                                                                   aaGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTtata                                                                <  2:560179/149‑1 (MQ=255)
                                                                                   aaGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTATAAACTATTTATATGATTTCCTtata                                                                >  1:276849/1‑149 (MQ=255)
                                                                                         tACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTaaa                                                          <  2:854737/149‑1 (MQ=255)
                                                                                          aCAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAAt                                                         >  1:733006/1‑149 (MQ=255)
                                                                                           cAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAAtt                                                        >  2:523168/1‑149 (MQ=255)
                                                                                                  aaagatTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAAGTTAAGGAACTTTGTGAACACCGCCATATTTCCATAGACACGTGATGATATTTACAGCAATTTTAATCTATTGATATGATTTCCTTATATTTAAATTAACTaaa                                                 <  1:182597/149‑1 (MQ=255)
                                                                                                      atTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGa                                             >  2:115611/1‑149 (MQ=255)
                                                                                                            tttCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAAtttt                                       <  1:738578/149‑1 (MQ=255)
                                                                                                                     aaaaaaaTATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGAt                              <  2:346388/148‑1 (MQ=255)
                                                                                                                       aaaaaTATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATgg                            >  2:798471/1‑149 (MQ=255)
                                                                                                                        aaaaTATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATgta                           >  2:398266/1‑147 (MQ=255)
                                                                                                                           atatGTAACCAAAAGTAAAATGTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACCGAAATTTTGTTTCTGGTGGAAAc                        <  2:621663/149‑1 (MQ=255)
                                                                                                                               gTAACCAAAAGTAAAATTTTAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTa                    <  1:836178/149‑1 (MQ=255)
                                                                                                                               gTAACCAAAAGTAAAATTTTAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTa                    <  1:508274/149‑1 (MQ=255)
                                                                                                                                        aGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTgg           <  1:523168/149‑1 (MQ=255)
                                                                                                                                             acgttAAGGCACTTCGTGAATACCGTAATGTTTCCATAGAGACGTGCTGATATTTACAGCAA‑TTTACTCTATTTATATGATTTCTGTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACaaa     <  2:268404/146‑1 (MQ=255)
                                                                                                                                              aTTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACaaa     <  2:886419/149‑1 (MQ=255)
                                                                                                                                                 tAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAAtc  <  2:687731/149‑1 (MQ=255)
                                                                                                                                                 tAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTACATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAAtc  <  1:602384/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATC  >  NC_000913/987051‑987344

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑