Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,928,260 T→A I15N (ATC→AAC)  rho → transcription termination factor Rho

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,928,2600TA98.0% 192.0 / NA 52I15N (ATC→AAC) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base T (1/0);  major base A (21/29);  minor base G (0/1);  total (22/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCT  >  REL606/3928187‑3928332
                                                                         |                                                                        
tttACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACa                                                                        >  2:567639/1‑76 (MQ=255)
 ttACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAacac                                                                       >  1:41742/1‑76 (MQ=255)
        aaGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGc                                                                <  1:277349/76‑1 (MQ=255)
        aaGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTACTGAGCTGAACACTCTCGGc                                                                <  1:9554/76‑1 (MQ=255)
         aGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCg                                                               <  2:775747/76‑1 (MQ=255)
           tttAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaa                                                             >  1:810689/1‑76 (MQ=255)
           tttAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaa                                                             >  2:85496/1‑76 (MQ=255)
             tAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaaaa                                                           >  2:1036376/1‑76 (MQ=255)
             tAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaaaa                                                           >  1:340090/1‑76 (MQ=255)
             tAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaaaa                                                           <  1:220437/76‑1 (MQ=255)
              aaGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAat                                                          <  1:319563/76‑1 (MQ=255)
                     cacacaACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGc                                                   >  2:380726/1‑76 (MQ=255)
                     cacaCCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGc                                                   >  2:1319128/1‑76 (MQ=255)
                      acacCACTATGAATCTTACCGAATTAAAGAATNCGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGNGCt                                                  >  2:1333231/1‑76 (MQ=255)
                      acacCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCt                                                  >  2:1031037/1‑76 (MQ=255)
                      acacCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCt                                                  >  1:819231/1‑76 (MQ=255)
                         ccacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGNGCTGAACACTCTCGGCGAAAATATGGGGCTGGa                                               <  2:1211749/76‑1 (MQ=255)
                          cacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGaa                                              <  2:1123517/76‑1 (MQ=255)
                           acTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGaaa                                             <  1:1392110/76‑1 (MQ=255)
                           acTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGaaa                                             <  1:755578/76‑1 (MQ=255)
                                gAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTg                                        <  2:13435/76‑1 (MQ=255)
                                   tCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCt                                     <  2:915497/76‑1 (MQ=255)
                                    cTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGcac                                    >  2:39335/1‑74 (MQ=255)
                                    cTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTc                                    <  1:1204231/76‑1 (MQ=255)
                                      tACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTAGGTCGt                                  <  1:85496/76‑1 (MQ=255)
                                       nccGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGNAAACCTGGCTCGTa                                 >  2:203603/2‑76 (MQ=255)
                                        ccGAATTAAAGAGTACGCCGGTTTCTGAGCTGNACACTCTCNGCGAAAATATGGGGCTGGGNAACCTGGCTCGtct                                >  2:1170447/1‑74 (MQ=255)
                                           aaTTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCg                             <  2:19255/76‑1 (MQ=255)
                                           aaTTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCg                             <  2:844628/76‑1 (MQ=255)
                                            aTTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGNAAACCTGGCTCGTATGCGt                            <  1:1325760/76‑1 (MQ=255)
                                            aTTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGt                            >  1:1064933/1‑76 (MQ=255)
                                            aTTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGt                            <  1:512438/76‑1 (MQ=255)
                                            aTTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGt                            <  2:903946/76‑1 (MQ=255)
                                             ttAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTa                           <  2:1369518/76‑1 (MQ=255)
                                             ttAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTa                           >  2:299323/1‑76 (MQ=255)
                                              tAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGCAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTaa                          <  1:134904/76‑1 (MQ=255)
                                               aaaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAg                         <  1:369248/76‑1 (MQ=255)
                                                aaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGc                        <  1:756293/76‑1 (MQ=255)
                                                aaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGc                        <  2:326897/76‑1 (MQ=255)
                                                aaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTANGCGTAAGc                        >  1:1252865/1‑76 (MQ=255)
                                                      aCGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACa                  <  1:252363/76‑1 (MQ=255)
                                                       cGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACat                 <  2:562647/76‑1 (MQ=255)
                                                         ccGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACatta               <  1:1202049/76‑1 (MQ=255)
                                                         ccGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACatta               >  1:203567/1‑76 (MQ=255)
                                                          cGGTTTCTGAGCTGAGCACTCTGGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACattat              <  2:46816/76‑1 (MQ=255)
                                                          cGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACattat              <  1:867162/76‑1 (MQ=255)
                                                          cGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACattat              <  1:1229162/76‑1 (MQ=255)
                                                               tCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTg         >  2:1020797/1‑76 (MQ=255)
                                                                cTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGc        <  2:799577/76‑1 (MQ=255)
                                                                 tGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGcc       >  1:348926/1‑76 (MQ=255)
                                                                  gAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCa      >  2:256350/1‑76 (MQ=255)
                                                                  gAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCa      >  1:838203/1‑76 (MQ=255)
                                                                   aGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCAt     <  1:1231483/76‑1 (MQ=255)
                                                                      tGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCt  >  2:359899/1‑76 (MQ=255)
                                                                      tGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCt  >  2:1023071/1‑76 (MQ=255)
                                                                         |                                                                        
TTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCT  >  REL606/3928187‑3928332

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