Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,005,601 +C coding (237/282 nt) b4568 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,005,6001.C98.4% 323.4 / ‑5.6 125G79A (GGC→GCC) b4568hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  major base C (68/54);  minor base A (2/0);  total (71/54)
Fisher's exact test for biased strand distribution p-value = 5.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 130 total aligned reads displayed.

GGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTG  >  NC_000913/4005525‑4005670
                                                                             |                                                                      
ggATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTgg                                                                         >  2:3237668/1‑76 (MQ=255)
  aTTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGcc                                                                      >  2:1402666/1‑75 (MQ=255)
    tGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGcctg                                                                    <  2:2283552/76‑4 (MQ=255)
          gctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTAtctc                                                               >  1:2860786/1‑76 (MQ=255)
          gctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTAtctc                                                               >  2:1715383/1‑76 (MQ=255)
           ctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATctct                                                              <  1:3852619/76‑1 (MQ=255)
           ctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATctct                                                              <  2:3129712/76‑1 (MQ=255)
            ttGGCGTGCATTTGATAATCGAAAAAATGGTCTGGTTGCATCCGCTGATTATGGTGGGGGGTGGCCTGTATCTctg                                                             <  2:4773426/76‑1 (MQ=255)
            ttGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctg                                                             >  2:1694007/1‑76 (MQ=255)
             tgtccTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctgc                                                            <  1:4539610/73‑1 (MQ=255)
             tGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctgc                                                            <  2:3016873/76‑1 (MQ=255)
              ggCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACTCTGATTATGGTGGGCGGTGGCCTGTATCTctgct                                                           >  1:4605058/1‑76 (MQ=255)
              ggCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctgct                                                           >  2:4508984/1‑76 (MQ=255)
               gCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctgctg                                                          >  1:4203049/1‑76 (MQ=255)
               gCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctgctg                                                          <  2:420502/76‑1 (MQ=255)
                ccTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTgg                                                         >  1:1639966/1‑76 (MQ=255)
                ccTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTgg                                                         >  1:2554408/1‑76 (MQ=255)
                 cTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTggg                                                        >  1:574686/1‑75 (MQ=255)
                 cTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGa                                                        >  1:1880748/1‑76 (MQ=255)
                 cTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGa                                                        >  1:2355708/1‑76 (MQ=255)
                 cTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGa                                                        <  2:316849/76‑1 (MQ=255)
                  tGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGAt                                                       >  2:2361048/1‑76 (MQ=255)
                  tGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGGGGCCTGTATCTCTGCTGGAt                                                       >  2:4421031/1‑76 (MQ=255)
                   gCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATg                                                      <  1:620312/76‑1 (MQ=255)
                   gCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGAGGTGGACTGTATCTCTGCTGGATg                                                      >  1:1174222/1‑76 (MQ=255)
                    cATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATgg                                                     >  2:773925/1‑76 (MQ=255)
                    cATTTGATCATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATgg                                                     >  2:4141152/1‑76 (MQ=255)
                     aTTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATggg                                                    >  1:1034854/1‑76 (MQ=255)
                     aTTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATggg                                                    <  2:1039362/76‑1 (MQ=255)
                        tGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTa                                                 <  1:2060298/76‑1 (MQ=255)
                         gATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTAc                                                <  1:3061513/76‑1 (MQ=255)
                           ttATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCa                                              <  2:669291/76‑1 (MQ=255)
                           ttATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGATGGATGGGTTACCa                                              >  2:1977096/1‑76 (MQ=255)
                              tCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGAt                                           <  2:1129219/76‑1 (MQ=255)
                              tCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGCTACCAGAt                                           >  2:4631613/1‑76 (MQ=255)
                               cGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATg                                          >  2:1239189/1‑76 (MQ=255)
                               cGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATg                                          >  2:2605619/1‑76 (MQ=255)
                                gAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGc                                         >  1:3056059/1‑76 (MQ=255)
                                gAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGc                                         <  1:537873/76‑1 (MQ=255)
                                gAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGc                                         >  1:3107958/1‑76 (MQ=255)
                                 aaaaaaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCt                                        <  2:216013/76‑1 (MQ=255)
                                 aaaaaaTGGCCTGGCTGCATACGCTGATTATGGTGAGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCt                                        <  1:10520/76‑1 (MQ=255)
                                   aaaaTGGCGTGGGGGGGTACGCTGATTATGGTGGGGGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTAc                                      <  2:279261/76‑1 (MQ=255)
                                   aaaaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTAc                                      <  1:689836/76‑1 (MQ=255)
                                   aaaaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTAc                                      <  2:2353715/76‑1 (MQ=255)
                                   aaaaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTAc                                      <  1:1221623/76‑1 (MQ=255)
                                    aaaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGCTGCTACg                                     >  2:704962/1‑76 (MQ=255)
                                    aaaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACg                                     <  1:2450453/76‑1 (MQ=255)
                                     aaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgt                                    <  1:4696055/76‑1 (MQ=255)
                                      aTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtg                                   >  1:2677189/1‑76 (MQ=255)
                                        ggCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtggt                                 <  2:3721350/76‑1 (MQ=255)
                                         gCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtggtg                                >  2:4298706/1‑76 (MQ=255)
                                           cTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCa                              >  2:4271172/1‑76 (MQ=255)
                                           cTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCa                              >  2:2433053/1‑76 (MQ=255)
                                           cTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGCTACCAGATGCTACGTGGTGCa                              >  2:4267647/1‑76 (MQ=255)
                                            tGGGTGCATACGCGGATTAGGGTGGGCGGGGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCAc                             <  2:4776777/76‑1 (MQ=255)
                                            tGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCAc                             >  2:2167123/1‑76 (MQ=255)
                                            tGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCAc                             >  2:646738/1‑76 (MQ=255)
                                             ggCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACt                            <  2:2587623/76‑1 (MQ=255)
                                             ggCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACt                            >  1:787217/1‑76 (MQ=255)
                                              gCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTg                           >  2:3307050/1‑76 (MQ=255)
                                              gCTGCACACGCTGATTATGGTGGGCGGTCGCCTGTACCTCCGCCGGATGGGTTACCAGATCCAACGTGGTGCACTg                           >  2:3958614/1‑76 (MQ=255)
                                                tGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaa                         >  1:464782/1‑76 (MQ=255)
                                                tGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaa                         >  1:819722/1‑76 (MQ=255)
                                                 gCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGCTGGGTTACCAGATGCTACGCGGTGCACCGaaa                        >  2:1777450/1‑76 (MQ=255)
                                                 gCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaaa                        <  2:1474915/76‑1 (MQ=255)
                                                 gCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaaa                        <  2:2538344/76‑1 (MQ=255)
                                                 gCAGACGCGGAGTAGGGGGGGCGGGGGCCTGTATCGCTGCTGGATGGGTTACCAGATGCGACGTGGTGCACTGaaa                        <  2:2559417/76‑1 (MQ=255)
                                                  cATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGgaaa                       >  1:1456242/1‑76 (MQ=255)
                                                  cATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaaaa                       >  2:3722771/1‑76 (MQ=255)
                                                  cATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaaaa                       >  1:2799802/1‑76 (MQ=255)
                                                    tACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaaaaaa                     <  2:536603/76‑1 (MQ=255)
                                                    tACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaaaaaa                     <  2:2197577/76‑1 (MQ=255)
                                                    tACGCTGATTATGGTGGGCGGTGGACTGTATCTCTGCTGGATGGGGTACCAGATGCTACGCGCTGCACTGaaaaaa                     >  1:1543061/1‑76 (MQ=255)
                                                     aCGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAag                    >  2:2045236/1‑76 (MQ=255)
                                                       gCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAagag                  >  1:3201180/1‑76 (MQ=255)
                                                         tGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGGGCACTGAAAAAAGAggg                >  2:4603056/1‑75 (MQ=255)
                                                          gATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCg               <  2:2314185/76‑1 (MQ=255)
                                                          gATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCg               <  1:773925/76‑1 (MQ=255)
                                                          gATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCg               <  1:2999775/76‑1 (MQ=255)
                                                          gATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCg               <  1:3237668/76‑1 (MQ=255)
                                                          gATTATGGGGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCg               >  2:1736867/1‑76 (MQ=255)
                                                           aTTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCgg              >  2:2698806/1‑76 (MQ=255)
                                                           aTTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCgg              <  1:582446/76‑1 (MQ=255)
                                                           aTTATGGTGGGCGGCGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCgg              >  1:3719291/1‑76 (MQ=255)
                                                              aTGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGttt           <  1:168060/76‑1 (MQ=255)
                                                                ggtggGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCATATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCt         <  1:2584118/76‑1 (MQ=255)
                                                                ggtggGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCt         <  1:2045237/76‑1 (MQ=255)
                                                                ggtggGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCt         <  1:4298706/76‑1 (MQ=255)
                                                                ggtggGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCt         <  1:3357514/76‑1 (MQ=255)
                                                                 gtggGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTg        <  2:4203779/76‑1 (MQ=255)
                                                                  tggTCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGc       <  1:152652/76‑1 (MQ=255)
                                                                  ggggCGGTGGCCTGTATCTCTGCTGGATGGGTTACGAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGc       <  2:761217/75‑1 (MQ=255)
                                                                   gggCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGGGTCTGcc      >  2:2862089/1‑75 (MQ=255)
                                                                     gCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCAcc    <  2:1737982/76‑1 (MQ=255)
                                                                      cGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCt   <  2:2674792/76‑1 (MQ=255)
                                                                      cGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGGGTCTGcccc  >  2:455649/1‑75 (MQ=255)
                                                                       ggTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTTTGCCCCtc  >  2:1444662/1‑75 (MQ=255)
                                                                       ggTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTg  <  1:4719811/76‑1 (MQ=255)
                                                                       ggTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGGTTTTTCACCTg  >  1:3307076/1‑76 (MQ=255)
                                                                             |                                                                      
GGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTG  >  NC_000913/4005525‑4005670

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