Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,101,777 A→C intergenic (‑7/+86) hslV ← / ← ftsN ATP‑dependent protease peptidase subunit/essential cell division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,101,7770AC96.2% 203.5 / ‑1.6 53intergenic (‑7/+86)hslV/ftsNATP‑dependent protease peptidase subunit/essential cell division protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (21/30);  total (22/31)
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.30e-01

ACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGG  >  REL606/4101704‑4101847
                                                                         |                                                                      
aCGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGccc                                                                      >  1:1310091/1‑76 (MQ=255)
    ggCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGCCCCCtt                                                                  <  2:1441077/76‑1 (MQ=255)
     gCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGCCCCCTTg                                                                 <  2:817426/76‑1 (MQ=255)
     gCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGCCCCCTTg                                                                 <  1:1416273/76‑1 (MQ=255)
     gCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGCCCCCTTg                                                                 <  1:389488/76‑1 (MQ=255)
        tgcccATCGCCAGCGATTGCCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGCCCCCTTGGtt                                                              <  2:166827/73‑1 (MQ=255)
         gACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGCCCCCTTGGTTa                                                             >  2:993561/1‑76 (MQ=255)
         gACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGCCCCCTTGGTTa                                                             >  2:440300/1‑76 (MQ=255)
           ccATCGCCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGCCCCCTTGGTTACg                                                           <  2:601333/76‑1 (MQ=255)
                       aTGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTAc                                               <  1:1136503/76‑1 (MQ=255)
                       aTGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTAc                                               >  1:541191/1‑76 (MQ=255)
                         gACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACaa                                             <  1:1525199/76‑1 (MQ=255)
                          aCCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACaaa                                            <  2:571400/76‑1 (MQ=255)
                            cacaTGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAAcc                                          <  1:942170/76‑1 (MQ=255)
                              caTGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGCCCCCTTGGTTACGAATACCGAGTACAAAcccc                                        >  2:437100/1‑76 (MQ=255)
                                tGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGt                                      <  1:276209/76‑1 (MQ=255)
                                 ggCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTa                                     >  2:1459320/1‑76 (MQ=255)
                                  gCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTAc                                    <  1:963175/76‑1 (MQ=255)
                                  gCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTAc                                    <  1:868447/76‑1 (MQ=255)
                                   ccGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACa                                   >  1:357316/1‑76 (MQ=255)
                                        aCGGCGGGCGCTTACTATAGTTGTCACGAGCTGCCCCCGTGGTTACGAATACAGAGTACAAACCACGTACAAAAGt                              <  1:101247/76‑1 (MQ=255)
                                         cGGCGTACGCTTACTATAGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTa                             <  2:1278548/76‑1 (MQ=255)
                                         cGGCGTACGCTTACTATAGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTa                             >  2:1643358/1‑76 (MQ=255)
                                             gTACGCTTACTATAGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACggg                         <  2:338913/76‑1 (MQ=255)
                                              tACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACgggg                        >  1:1228049/1‑76 (MQ=255)
                                              tACCCTTACTATAGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACgggg                        <  2:807371/76‑1 (MQ=255)
                                                cGCTTACTATAGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCa                      >  1:1247968/1‑76 (MQ=255)
                                                 gCTTACTATAGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCAt                     <  2:1580858/76‑1 (MQ=255)
                                                   ttACTATAGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATaa                   >  1:1697171/1‑76 (MQ=255)
                                                   ttACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATaa                   <  1:1563469/76‑1 (MQ=255)
                                                    tACTGTAGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTGCAAAAGTACGGGGCATAAt                  <  2:668883/76‑1 (MQ=255)
                                                    tACTATAGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAAt                  >  1:1360460/1‑76 (MQ=255)
                                                    tACTATAGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAAt                  <  1:453823/76‑1 (MQ=255)
                                                    tACTATAGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGAACGGGGCATAAt                  >  1:255639/1‑76 (MQ=255)
                                                    tACTATAGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGCACAAAAATACGGGGCATAAt                  >  1:1595953/1‑76 (MQ=255)
                                                    tACCATAGTTGTCACGCGCTGCCCCCTTGGTTACGAATACAGAGTACAACCCCCGTACAAAAGTACGGGGCATAAt                  <  2:581408/76‑1 (MQ=255)
                                                     aCTATAGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATg                 >  1:552449/1‑76 (MQ=255)
                                                     aCTATAGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATg                 <  1:1075491/76‑1 (MQ=255)
                                                     aCTATAGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATg                 <  1:1222836/76‑1 (MQ=255)
                                                     aCTATAGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATg                 <  1:194602/76‑1 (MQ=255)
                                                     aCTATAGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATg                 <  2:347917/76‑1 (MQ=255)
                                                     aCTATAGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATg                 >  1:754962/1‑76 (MQ=255)
                                                         taGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAAt             >  1:441361/1‑76 (MQ=255)
                                                          aGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAAtt            <  2:1013100/76‑1 (MQ=255)
                                                          aGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAAtt            >  1:879838/1‑76 (MQ=255)
                                                          aGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAAtt            >  1:947938/1‑76 (MQ=255)
                                                          aGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAAtt            >  2:1550987/1‑76 (MQ=255)
                                                          aGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAAAACGGGGCATAAGGCAAtt            >  1:687333/1‑76 (MQ=255)
                                                           gTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATta           >  1:1381464/1‑76 (MQ=255)
                                                           gTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATta           <  2:848111/76‑1 (MQ=255)
                                                            ttGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATtat          >  2:233525/1‑76 (MQ=255)
                                                            ttGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATtat          <  2:520884/76‑1 (MQ=255)
                                                            ttGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATtat          <  2:554729/76‑1 (MQ=255)
                                                              gTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAg        <  2:1009737/76‑1 (MQ=255)
                                                                 aggAGCTGCCCCCCTCGTTACGAAAACAGAGAACAAGCGCAGTAGAAAAGGCAAGAGCATAATGCAATTATATATg     <  2:498227/74‑1 (MQ=255)
                                                                  cGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATgg    >  2:1655113/1‑76 (MQ=255)
                                                                    aGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATgggg  >  1:290327/1‑76 (MQ=255)
                                                                         |                                                                      
ACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGG  >  REL606/4101704‑4101847

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