Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,883,147 (A)5→4 intergenic (‑341/+276) atpI ← / ← gidB F0F1 ATP synthase subunit I/glucose‑inhibited division protein; b3740_1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,883,1430A.92.3% 218.2 / 13.8 52intergenic (‑337/+280)atpI/gidBF0F1 ATP synthase subunit I/glucose‑inhibited division protein; b3740_1
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base . (26/22);  total (29/23)
Fisher's exact test for biased strand distribution p-value = 6.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

TGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAACCACCCATTGACATTTTTAATAATGT  >  REL606/3883070‑3883212
                                                                         |                                                                     
tGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTaaa                                                                     <  1:1547583/76‑1 (MQ=255)
tGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTaaa                                                                     <  1:1569436/76‑1 (MQ=255)
 gcgcATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTaaaa                                                                    <  2:1080119/76‑1 (MQ=255)
 gcgcATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTaaaa                                                                    <  2:868577/76‑1 (MQ=255)
  cgcATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTaaaac                                                                   >  2:1155066/1‑75 (MQ=255)
   gcATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTaaaaca                                                                  >  1:1415053/1‑74 (MQ=255)
    cATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTaaaacag                                                                 <  1:1487637/76‑4 (MQ=255)
    cATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTaaaacag                                                                 <  1:640394/76‑4 (MQ=255)
      tGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAACAGAt                                                              >  2:450461/1‑76 (MQ=255)
       gTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAACAGATa                                                             <  2:1539687/76‑1 (MQ=255)
       gTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAACAGATa                                                             <  1:72450/76‑1 (MQ=255)
         tCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAACAGATACt                                                           >  2:145553/1‑76 (MQ=255)
            cgcACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAACAGATACTGtt                                                        >  1:1239404/1‑76 (MQ=255)
              cacaGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAACAGATACTGtttt                                                      <  2:1010796/76‑1 (MQ=255)
              cacaGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAACAGATACTGtttt                                                      <  2:1212080/76‑1 (MQ=255)
               acaGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAACAGATACTGTTTta                                                     <  2:336389/76‑1 (MQ=255)
                caGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAACAGATACTGTTTtaa                                                    <  2:152757/76‑1 (MQ=255)
                 aGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAACAGATACTGTTTtaat                                                   <  1:178931/76‑1 (MQ=255)
                   cATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACAGATACTGTTTtaata                                                  <  2:915993/76‑1 (MQ=255)
                     tatttagtTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAAt                                               <  2:1097823/76‑1 (MQ=255)
                           tttaCTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACAtt                                         <  2:665008/76‑1 (MQ=255)
                            ttaCTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACAttt                                        <  1:571163/76‑1 (MQ=255)
                            ttaCTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACAttt                                        >  1:1687005/1‑76 (MQ=255)
                            ttaCTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACAttt                                        >  2:1032944/1‑76 (MQ=255)
                                 tGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACAGATACTGTTTTAATAAATGACATTTacac                                    >  2:316406/1‑76 (MQ=255)
                                 tGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTAcaca                                   >  1:1707466/1‑76 (MQ=255)
                                 tGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTAcaca                                   >  2:1180186/1‑76 (MQ=255)
                                   gCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACacaac                                 >  2:434528/1‑76 (MQ=255)
                                    cAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACacaaca                                >  2:1728930/1‑76 (MQ=255)
                                    cAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACacaaca                                >  2:960927/1‑76 (MQ=255)
                                      aaTGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACACAACaaa                              <  1:1113453/76‑1 (MQ=255)
                                       atgatgCCTTTGCAAGTTTATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACACAACaaaa                             <  2:1527101/76‑1 (MQ=255)
                                        tgatgCCTTTGCAAGTTTATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACACAACaaaaa                            >  1:1047281/1‑76 (MQ=255)
                                        tgatgCCTTTGCAAGTTTATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACACAACaaaaa                            <  2:546523/76‑1 (MQ=255)
                                          atgCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAac                           >  2:34089/1‑76 (MQ=255)
                                          atgCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAac                           >  2:34057/1‑76 (MQ=255)
                                          atgCCTTTGCAAGTTTATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAacc                          >  1:661246/1‑76 (MQ=255)
                                          atgCCTTTGCAAGTTTATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAacc                          <  1:698348/76‑1 (MQ=255)
                                            gCCTTTGCAAGTTTATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAaccac                        >  2:1671925/1‑76 (MQ=255)
                                             ccTTTGCAAGTTTATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAaccacc                       >  1:1170920/1‑76 (MQ=255)
                                               tttGCAAGTTTATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAACCACCCa                     <  1:99860/76‑1 (MQ=255)
                                                ttGCAAGTTTATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAACCACCCAt                    >  1:1365291/1‑76 (MQ=255)
                                                 tGCAAGTTTATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAACCACCCAtt                   >  1:760012/1‑76 (MQ=255)
                                                   cAAGTTTATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAACCACCCATTGa                 <  2:988549/76‑1 (MQ=255)
                                                     aGTTTATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAACCACCCATTGACa               >  1:83244/1‑76 (MQ=255)
                                                     aGTTTATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAACCACCCATTGACa               <  2:946244/76‑1 (MQ=255)
                                                      gTTTATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAACCACCCATTGACAt              <  2:1447407/76‑1 (MQ=255)
                                                       tttATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAACCACCCATTGACAtt             >  2:567313/1‑76 (MQ=255)
                                                       tttATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAACCACCCATTGACAtt             >  2:82667/1‑76 (MQ=255)
                                                        ttATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAACCACCCATTGACAttt            >  2:609120/1‑76 (MQ=255)
                                                        ttATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAACCACCCATTGACAttt            >  1:1680401/1‑76 (MQ=255)
                                                        ttATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAACCACCCATTGACAttt            >  1:335989/1‑76 (MQ=255)
                                                         tATGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAACCACCCATTGACAtttt           >  2:723155/1‑76 (MQ=255)
                                                          aTGATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAACCACCCATTGACAttttt          <  1:1216381/76‑1 (MQ=255)
                                                            gATATTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAACCACCCATTGACATTTTtaa        <  2:1137478/76‑1 (MQ=255)
                                                             atatTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAACCACCCATTGACATTTTtaat       >  2:1640578/1‑76 (MQ=255)
                                                              tatTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAACCACCCATTGACATTTTtaata      >  1:175394/1‑76 (MQ=255)
                                                              tatTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAACCACCCATTGACATTTTtaata      <  2:768263/76‑1 (MQ=255)
                                                              tatTTCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAACCACCCATTGACATTTTtaata      >  2:101075/1‑76 (MQ=255)
                                                                  tCAGTCTAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAACCACCCATTGACATTTTTAATAATGt  <  2:63116/76‑1 (MQ=255)
                                                                         |                                                                     
TGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAACCACCCATTGACATTTTTAATAATGT  >  REL606/3883070‑3883212

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