Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,068,952 C→T *302* (TAG→TAA insD1 ← IS2 transposase TnpB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,068,9520CT98.2% 177.6 / ‑4.6 57*302* (TAG→TAAinsD1IS2 transposase TnpB
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (28/28);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCT  >  NC_000913/2068818‑2069075
                                                                                                                                      |                                                                                                                           
aTTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTtata                                                                                                                         >  2:217682/1‑139 (MQ=255)
  tGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATAtt                                                                                                                       >  1:990623/1‑139 (MQ=255)
  tGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATAtt                                                                                                                       >  2:1592702/1‑139 (MQ=255)
       aacaaacaTTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAg                                                                                                                  <  1:1592702/139‑1 (MQ=255)
                tGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGtt                                                                                                         >  2:248566/1‑139 (MQ=255)
                 gCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTa                                                                                                        >  1:148814/1‑139 (MQ=255)
                 gCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTa                                                                                                        >  1:231312/1‑139 (MQ=255)
                          caccaATGCCTCCATAGCCGATAATTTGCAGGACAACCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAc                                                                                               >  1:19858/1‑139 (MQ=255)
                                 gCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                              <  1:418663/133‑1 (MQ=255)
                                 gCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                              >  2:418663/1‑133 (MQ=255)
                                 gCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTAc                                                                                        >  2:1563904/1‑139 (MQ=255)
                                        tAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                                 <  2:243927/139‑1 (MQ=255)
                                             gATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGcc                                                                            <  1:1563904/139‑1 (MQ=255)
                                              aTAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCg                                                                           >  1:1207970/1‑139 (MQ=255)
                                              aTAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATCTCCAGACAt                                                                                                              <  2:726490/104‑1 (MQ=255)
                                              aTAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATCTCCAGACAt                                                                                                              >  1:726490/1‑104 (MQ=255)
                                               tAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGACTGGATTTGCCCTTATATTTCCAGACATCTGTTATCGCTTAACCCATTACAAGCCCGCTGCCGc                                                                          >  1:1468410/1‑139 (MQ=255)
                                                             tCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTgg                                                            <  2:1207970/139‑1 (MQ=255)
                                                              ccccGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGa                                                                                                                 <  1:234474/85‑1 (MQ=255)
                                                              ccccGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGa                                                                                                                 >  2:234474/1‑85 (MQ=255)
                                                               cccGGAGTAGAGATAACTGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATCTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgct                                                          <  2:990623/139‑2 (MQ=255)
                                                                   gAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACAt                                                                                                              >  1:477839/1‑83 (MQ=255)
                                                                   gAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACAt                                                                                                              <  2:477839/83‑1 (MQ=255)
                                                                        gagaTAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACAt                                                                                                              >  1:849314/1‑78 (MQ=255)
                                                                        gagaTAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACAt                                                                                                              <  2:849314/78‑1 (MQ=255)
                                                                        gagaTAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCg                                                                           <  1:875716/113‑1 (MQ=255)
                                                                        gagaTAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCg                                                                           >  2:875716/1‑113 (MQ=255)
                                                                              aCAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgcg                                           >  1:922010/1‑139 (MQ=255)
                                                                                aGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                                   >  2:286744/1‑81 (MQ=255)
                                                                                aGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                                   <  1:286744/81‑1 (MQ=255)
                                                                                aGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCa                                         <  2:148814/139‑1 (MQ=255)
                                                                                aGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCa                                         <  2:231312/139‑1 (MQ=255)
                                                                                  aCAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTtat                                                                                                                          <  1:441516/56‑1 (MQ=255)
                                                                                  aCAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTtat                                                                                                                          >  2:441516/1‑56 (MQ=255)
                                                                                     aCTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGcc                                                                            <  1:1448767/99‑1 (MQ=255)
                                                                                     aCTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGcc                                                                            >  2:1448767/1‑99 (MQ=255)
                                                                                         cTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGa                                <  1:248566/139‑1 (MQ=255)
                                                                                          ttGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACtt                                                                                                  <  2:933829/72‑1 (MQ=255)
                                                                                          ttGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACtt                                                                                                  >  1:933829/1‑72 (MQ=255)
                                                                                             gCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                            <  2:19858/139‑1 (MQ=255)
                                                                                                ggCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCAtt                         <  1:336162/139‑1 (MQ=255)
                                                                                                 gCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgcg                                           >  1:1446707/1‑120 (MQ=255)
                                                                                                 gCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgcg                                           <  2:1446707/120‑1 (MQ=255)
                                                                                                   cGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCa                                                                                                                   <  1:597712/46‑1 (MQ=255)
                                                                                                   cGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCa                                                                                                                   >  2:597712/1‑46 (MQ=255)
                                                                                                     tAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATg                              >  1:1181773/1‑129 (MQ=255)
                                                                                                     tAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATg                              <  2:1181773/129‑1 (MQ=255)
                                                                                                            aTCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGtt                                                                                                         <  1:887277/47‑1 (MQ=255)
                                                                                                            aTCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGtt                                                                                                         >  2:887277/1‑47 (MQ=255)
                                                                                                            aTCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTa                                                                                                        >  1:1531318/1‑48 (MQ=255)
                                                                                                            aTCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTa                                                                                                        <  2:1531318/48‑1 (MQ=255)
                                                                                                             tCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTc            >  2:453184/1‑139 (MQ=255)
                                                                                                                gCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGaa         <  2:922010/139‑1 (MQ=38)
                                                                                                                gCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGaa         <  1:217682/139‑1 (MQ=255)
                                                                                                                    tCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCa                           <  2:1497581/117‑1 (MQ=25)
                                                                                                                    tCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCa                           >  1:1497581/1‑117 (MQ=25)
                                                                                                                       tGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCctct  >  2:1343610/1‑139 (MQ=32)
                                                                                                                                      |                                                                                                                           
ATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCT  >  NC_000913/2068818‑2069075

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