Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,646,854 C→T intergenic (+117/+320) flxA → / ← intK Qin prophage; uncharacterized protein/pseudogene, integrase fragment, Qin prophage;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,646,8540CT98.3% 219.4 / ‑4.2 59intergenic (+117/+320)flxA/intKQin prophage; uncharacterized protein/pseudogene, integrase fragment, Qin prophage;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (34/24);  total (34/25)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

ATACCGCTGGCACCATTGATATTTACGTCTAAGTGACAGCCGTATTGTGGCCCTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATT  >  NC_000913/1646706‑1646998
                                                                                                                                                    |                                                                                                                                                
aTACCGCTGGCACCATTGATATTTACGTCTAAGTGACAGCCGTATTGTGGCCCTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGc                                                                                                                                                >  1:383500/1‑151 (MQ=255)
 tACCGCTGGCACCATTGATATTTACGTCTAAGTGACAGCCGTATTGTGGCCCTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGc                                                                                                                                                >  2:289276/1‑150 (MQ=255)
  aCCGCTGGCACCATTGATATTTACGTCTAAGTGACAGCCGTATTGTGGCCCTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCAc                                                                                                                                              >  2:662411/1‑151 (MQ=255)
   ccGCTGGCACCATTGATATTTACGTCTAAGTGACAGCCGTATTGTGGCCCTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACt                                                                                                                                             <  2:535760/151‑1 (MQ=255)
       tGGCACCATTGATATTTACGTCTAAGTGACAGCCGTATTGTGGCCCTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTccc                                                                                                                                          >  2:478628/1‑150 (MQ=255)
        ggCACCATTGATATTTACGTCTAAGTGACAGCCGTATTGTGGCCCTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTccc                                                                                                                                          <  2:947574/149‑1 (MQ=255)
             cATTGATATTTACGTCTAAGTGACAGCCGTATTGTGGCCCTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGt                                                                                                                                   <  1:617147/151‑1 (MQ=255)
                 gATATTTACGTCTAAGTGACAGCCGTATTGTGGCCCTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTAt                                                                                                                               <  1:330884/151‑1 (MQ=255)
                        aCGTCTAAGTGACAGCCGTATTGTGGCCCTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAAt                                                                                                                         <  1:736191/150‑1 (MQ=255)
                            cTAAGTGACAGCCGTATTGTGGCCCTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGcc                                                                                                                     >  1:129053/1‑150 (MQ=255)
                              aaGTGACAGCCGTATTGTGGCCCTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACg                                                                                                                  >  2:233941/1‑151 (MQ=255)
                                  gACAGCCGTATTGTGGCCCTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGt                                                                                                               <  1:977968/150‑1 (MQ=255)
                                   aCAGCCGTATTGTGGCCCTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGttt                                                                                                             >  2:20378/1‑151 (MQ=255)
                                     aGCCGTATTGTGGCCCTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTa                                                                                                            >  1:194739/1‑150 (MQ=255)
                                      gCCGTATTGTGGCCCTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAAc                                                                                                          <  1:748910/151‑1 (MQ=255)
                                        cGTATTGTGGCCCTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACt                                                                                                         <  2:153480/150‑1 (MQ=255)
                                                     tCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATc                                                                                           >  1:536147/1‑151 (MQ=255)
                                                      cATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGAt                                                                                            <  2:443469/149‑1 (MQ=255)
                                                         cGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAAtt                                                                                       <  2:870884/151‑1 (MQ=255)
                                                          gggCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACGACCCATGATCAATTa                                                                                      <  1:133940/151‑1 (MQ=255)
                                                          gggCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAAtt                                                                                       >  1:718688/1‑150 (MQ=255)
                                                           ggCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTAt                                                                                     >  2:93386/1‑151 (MQ=255)
                                                            gCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATTAATTATg                                                                                    <  2:140134/151‑1 (MQ=255)
                                                            gCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATg                                                                                    >  2:133480/1‑151 (MQ=255)
                                                                 ttttCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTg                                                                               >  2:143421/1‑151 (MQ=255)
                                                                 ttttCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTg                                                                               >  2:54838/1‑151 (MQ=255)
                                                                  tttCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAAtt                                                                                <  2:726581/149‑1 (MQ=255)
                                                                  tttCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTg                                                                               <  1:928423/150‑1 (MQ=255)
                                                                  tttCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTg                                                                               <  1:698374/150‑1 (MQ=255)
                                                                  tttCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGc                                                                              >  2:643575/1‑151 (MQ=255)
                                                                     cGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAAc                                                                           >  1:419835/1‑151 (MQ=255)
                                                                          tCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTAttt                                                                      >  2:463804/1‑151 (MQ=255)
                                                                           cAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTc                                                                     <  1:478679/151‑1 (MQ=255)
                                                                           cAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTc                                                                     <  1:131989/151‑1 (MQ=255)
                                                                             gCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTg                                                                   >  2:756860/1‑151 (MQ=255)
                                                                              ccTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGt                                                                  <  1:994927/151‑1 (MQ=255)
                                                                                 tttCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGt                                                               >  1:17539/1‑151 (MQ=255)
                                                                                     tttAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCACt                                                           >  1:103007/1‑151 (MQ=255)
                                                                                      ttAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCAc                                                            <  1:530601/149‑1 (MQ=255)
                                                                                      ttAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCACtt                                                          >  1:214571/1‑151 (MQ=255)
                                                                                      ttAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCACt                                                           <  2:239145/150‑1 (MQ=255)
                                                                                        aaaGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCACtttt                                                        >  1:588061/1‑151 (MQ=255)
                                                                                         aaGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCACTTTtg                                                       <  1:21066/151‑1 (MQ=255)
                                                                                          aGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCACTTTtgt                                                      >  2:50262/1‑151 (MQ=255)
                                                                                           gACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCACTTTtgtg                                                     >  1:296622/1‑151 (MQ=255)
                                                                                            aCATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCACTTTTGTgg                                                    <  2:394340/151‑1 (MQ=255)
                                                                                              atatTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCACTTTTGTgggg                                                  >  1:215992/1‑151 (MQ=255)
                                                                                                attatCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCACTTTTGTGGGGa                                                 >  2:545013/1‑150 (MQ=255)
                                                                                                 ttatCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCACTTTTGTGGGGAc                                                >  1:224265/1‑150 (MQ=255)
                                                                                                  tatCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCACTTTTGTGGGGACa                                               <  2:816263/150‑1 (MQ=255)
                                                                                                           aTCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCACTTTTGTGGGGACAGTCCACaa                                       >  1:446857/1‑149 (MQ=255)
                                                                                                             cATTTCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCACTTTTGTGGGGACAGTCCACAAAACt                                   >  1:949274/1‑151 (MQ=255)
                                                                                                                 tCTGATAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCACTTTTGTGGGGACAGTCCACAAAACTGcc                                >  2:290369/1‑150 (MQ=255)
                                                                                                                     aTAGTTAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACtt                           >  2:959080/1‑151 (MQ=255)
                                                                                                                         ttAACATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGc                       >  1:541008/1‑151 (MQ=255)
                                                                                                                            aCATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGcttc                    <  1:626905/151‑1 (MQ=255)
                                                                                                                             cATTACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGcttct                   >  1:961909/1‑151 (MQ=255)
                                                                                                                                tACAAGATATAAGTAATGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTg                 <  1:259104/150‑1 (MQ=255)
                                                                                                                                              aaTGGATGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATAGCAACTATTTCTGTAGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACAtt  >  2:561391/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
ATACCGCTGGCACCATTGATATTTACGTCTAAGTGACAGCCGTATTGTGGCCCTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATT  >  NC_000913/1646706‑1646998

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