Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 361,532 A→C L71R (CTG→CGG)  lacI ← DNA‑binding transcriptional repressor LacI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273361,5320AC90.4% 159.5 / 9.8 52L71R (CTG→CGG) lacIDNA‑binding transcriptional repressor LacI
Reads supporting (aligned to +/- strand):  ref base A (4/1);  new base C (25/22);  total (29/23)
Fisher's exact test for biased strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.55e-01

AAGATTGTGCACCGCCGCTTTACAGGCTTCGACGCCGCTTCGTTCTACCATCGACACCACCACGCTGGCACCCAGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCAGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCC  >  NZ_CP009273/361402‑361666
                                                                                                                                  |                                                                                                                                      
aagaTTGTGCACCGCCGCTTTACAGGCTTCGACGCCGCTTCGTTCTACCATCGACACCACCACGCTGGCACCCAGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAg                                                                                                                                <  1:126938/139‑1 (MQ=255)
       tgCACCGCCGCTTTACAGGCTTCGACGCCGCTTCGTTCTACCATCGACACCACCACGCTGGCACCCAGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCaa                                                                                                                         >  1:482063/1‑139 (MQ=255)
         cACCGCCGCTTTACAGGCTTCGACGCCGCTTCGTTCTACCATCGACACCACCACGCTGGCACCCAGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACg                                                                                                                       >  1:299624/1‑139 (MQ=255)
             gccgcTTTACAGGCTTCGACGCCGCTTCGTTCTACCATCGACACCACCACGCTGGCACCCAGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCaa                                                                                                                   <  1:305678/139‑1 (MQ=255)
             gccgcTTTACAGGCTTCGACGCCGCTTCGTTCTACCATCGACACCACCACGCTGGCACCCAGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCaa                                                                                                                   >  1:80074/1‑139 (MQ=255)
                 cTTTACAGGCTTCGACGCCGCTTCGTTCTACCATCGACACCACCACGCTGGCACCCAGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAg                                                                                                               >  2:651756/1‑139 (MQ=255)
                         gCTTCGACGCCGCTTCGTTCTACCATCGACACCACCACGCTGGCACCCAGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGTGCGTGCAGGGCCCGATTGGAGGCGGCAACGCCAATCAGCAACGACt                                                                                                       >  1:856305/1‑139 (MQ=255)
                               aCGCCGCTTCGTTCTACCATCGACACCACCACGCTGGCACCCAGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGACCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGttgcc                                                                                                 >  1:337883/1‑136 (MQ=255)
                                    gcTTCGTTCTACCATCGACACCACCACGCTGGCACCCAGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGcc                                                                                            <  2:212809/139‑1 (MQ=255)
                                     cTTCGTTCTACCATCGACACCACCACGCTGGCACCCAGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGc                                                                                             >  1:197829/1‑137 (MQ=255)
                                     cTTCGTTCTACCATCGACACCACCACGCTGGCACCCAGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGc                                                                                             >  2:805510/1‑137 (MQ=255)
                                          ttCTACCATCGACACCACCACGCTGGCACCCAGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCaa                                                                                                                   <  1:752900/110‑1 (MQ=255)
                                          ttCTACCATCGACACCACCACGCTGGCACCCAGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCaa                                                                                                                   >  2:752900/1‑110 (MQ=255)
                                                 aTCGACACCACCACGCTGGCACCCAGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCa                                                                               <  2:795293/139‑1 (MQ=255)
                                                  tCGACACCACCACGCTGGCACCCAGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCAGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCAc                                                                              >  1:66070/1‑139 (MQ=255)
                                                  tCGACACCACCACGCTGGCACCCAGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCAGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCAc                                                                              >  1:636834/1‑139 (MQ=255)
                                                       accaccacGCTGGCACCCAGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCa                                                                                                                    <  1:504809/96‑1 (MQ=255)
                                                       accaccacGCTGGCACCCAGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCa                                                                                                                    >  2:504809/1‑96 (MQ=255)
                                                        ccgccacGCTGGCACCCAGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGtt                                                                        <  1:161668/136‑1 (MQ=255)
                                                         caccacGCTGGCACCCAGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTg                                                                       >  1:396974/1‑139 (MQ=255)
                                                          accacGCTGGCACCCAGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTgg                                                                      <  1:117279/139‑1 (MQ=255)
                                                          accacGCTGGCACCCAGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTgg                                                                      <  1:87403/139‑1 (MQ=255)
                                                                   gCACCCAGTTGATCGGCGCGCGATTTAATCGCCGCGACACTTTGCGACGGCGCGTGCAGGGCCCGACTGGCGGTGGCAACGCCAATCAGCAACGACTGCTGGCCCGCCAGTTGTTGTGCCACCCGGTTGGGACTGTact                                                             >  2:806212/1‑137 (MQ=255)
                                                                   gCACCCAGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAAt                                                             >  2:64846/1‑139 (MQ=255)
                                                                         aGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGcc                                                                                            >  2:764733/1‑102 (MQ=255)
                                                                         aGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGcc                                                                                            <  1:764733/102‑1 (MQ=255)
                                                                                cGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCAt                                                >  1:856713/1‑139 (MQ=255)
                                                                                cGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGAAATTCAGCTCCGCCAt                                                >  2:833530/1‑139 (MQ=255)
                                                                                   cgcgAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGCAACGGCTGTTTGCCCGcc                                                                                            <  2:27651/92‑1 (MQ=255)
                                                                                   cgcgAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGcc                                                                                            >  1:27651/1‑92 (MQ=255)
                                                                                   cgcgAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATcgc                                             >  1:841477/1‑139 (MQ=255)
                                                                                     cgAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATcgccg                                           >  2:647983/1‑139 (MQ=255)
                                                                                       agaTTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCt                                         >  2:641484/1‑139 (MQ=255)
                                                                                        gaTTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCtt                                        <  1:348405/139‑1 (MQ=255)
                                                                                             aaTCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACt                                   >  1:31983/1‑139 (MQ=255)
                                                                                                  ccgcGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGc                                                                                                              <  2:388766/59‑1 (MQ=255)
                                                                                                  ccgcGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGc                                                                                                              >  1:388766/1‑59 (MQ=255)
                                                                                                  ccgcGACAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTc                              >  1:338426/1‑139 (MQ=255)
                                                                                                        cAATTTGCGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGt                        <  2:651281/139‑1 (MQ=255)
                                                                                                               cGACGGCGCGTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCa                 <  2:841477/139‑1 (MQ=255)
                                                                                                                     cgcgTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCAcg                                                                             <  1:478710/73‑1 (MQ=255)
                                                                                                                     cgcgTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCAcg                                                                             >  2:478710/1‑73 (MQ=255)
                                                                                                                        gTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGtgg        <  1:641484/139‑1 (MQ=255)
                                                                                                                        gTGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGtgg        <  2:338426/139‑1 (MQ=255)
                                                                                                                         tGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGtggc       <  2:724674/139‑1 (MQ=255)
                                                                                                                         tGCAGGGCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGtggc       <  2:358847/139‑1 (MQ=255)
                                                                                                                              ggCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGaaa             >  1:187116/1‑128 (MQ=255)
                                                                                                                              ggCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGaaa             <  2:187116/128‑1 (MQ=255)
                                                                                                                              ggCCCGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGcc  <  1:651756/139‑1 (MQ=255)
                                                                                                                              ggCCAGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGcc  >  1:550466/1‑139 (MQ=255)
                                                                                                                               gCCAGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCAcgc                                                                            >  2:752390/1‑64 (MQ=255)
                                                                                                                               gCCAGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCAcgc                                                                            <  1:752390/64‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                      
AAGATTGTGCACCGCCGCTTTACAGGCTTCGACGCCGCTTCGTTCTACCATCGACACCACCACGCTGGCACCCAGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCAGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCC  >  NZ_CP009273/361402‑361666

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