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.5% 360.8 / ‑3.1 139G79A (GGC→GCC) b4568hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  major base C (60/74);  minor base T (1/1);  total (63/76)
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.39e-01

Only 100 of 147 total aligned reads displayed.

GGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTG  >  NC_000913/4005525‑4005670
                                                                             |                                                                       
ggATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAGTGGCCTGGCTGCATACGCTGATTATGGTGGGCGgggg                                                                          >  2:1055316/1‑76 (MQ=255)
 taTTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGc                                                                        >  2:4443145/2‑76 (MQ=255)
 gATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGc                                                                        >  1:5513163/1‑76 (MQ=255)
   ttGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGcct                                                                      <  1:326480/76‑3 (MQ=255)
   ttGCGCTGCTTGGCCTGCACTTCATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTggtct                                                                       <  1:3743509/76‑4 (MQ=255)
       gtggctTTGCCTGCGTTTGATTATCTGAAAAATAGCCTGGCTGCATACGGTGGTTATGGTTGGGGGTGGTCTGTAt                                                                  <  1:151457/73‑1 (MQ=255)
       gctgctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTAt                                                                  <  1:2337878/76‑1 (MQ=255)
           ctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATctct                                                              <  2:921032/76‑1 (MQ=255)
           ctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATctct                                                              <  2:4403263/76‑1 (MQ=255)
            ttGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctg                                                             >  1:1365560/1‑76 (MQ=255)
            ttGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctg                                                             <  2:4396435/76‑1 (MQ=255)
             tGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctgc                                                            <  1:516190/76‑1 (MQ=255)
              tgtcTGCATTTGATTATGGAAAAAACGGCCTGGTTGCATTGGCTGTTTATGGTGGGCGGTGGCCTGTATCTgtgct                                                           <  1:2584290/73‑1 (MQ=255)
              ggCCTGCATTTGATTATCGAAAAAATGGCCTTGCTGCATACGCTGATTTTGGTGGGAGGAGGTCAGTATCTctgct                                                           >  1:2058007/1‑76 (MQ=255)
              ggCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGGCTGTATCTctgct                                                           >  1:4953461/1‑76 (MQ=255)
              ggCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctgct                                                           >  1:2216912/1‑76 (MQ=255)
              ggCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctgct                                                           <  1:3107284/76‑1 (MQ=255)
               gCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctgctg                                                          >  1:2091070/1‑76 (MQ=255)
               gCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctgctg                                                          <  1:1539058/76‑1 (MQ=255)
               gCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctgctg                                                          <  1:112624/76‑1 (MQ=255)
                ccTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTgg                                                         <  1:3588186/76‑1 (MQ=255)
                ccTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGAGGGGGCCTGTATCTCTGCTgg                                                         >  1:231557/1‑76 (MQ=255)
                 cTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGa                                                        >  1:2013239/1‑76 (MQ=255)
                 cTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGa                                                        >  1:2902237/1‑76 (MQ=255)
                   gcgtttGATTATCGAAAAAATGGCCGGGCGGCATACGCTGAGTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATg                                                      <  2:5672023/73‑1 (MQ=255)
                   gCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATg                                                      <  1:1660124/76‑1 (MQ=255)
                   gCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATg                                                      >  2:4627864/1‑76 (MQ=255)
                    cATTTGATTATCGAAAAAATGGCCTGGCTGCGTACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATgg                                                     <  2:3717612/76‑1 (MQ=255)
                    cATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATgg                                                     >  1:4515983/1‑76 (MQ=255)
                    cATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATgg                                                     >  1:4163172/1‑76 (MQ=255)
                    cATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATgg                                                     <  1:52877/76‑1 (MQ=255)
                     aTTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATggg                                                    >  2:3646856/1‑76 (MQ=255)
                      tttGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGt                                                   <  1:1676416/76‑1 (MQ=255)
                      tttGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGt                                                   <  2:2859832/76‑1 (MQ=255)
                      tttGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGt                                                   >  2:1839148/1‑76 (MQ=255)
                       ttGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGtt                                                  >  1:276002/1‑76 (MQ=255)
                        tGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCTTGTATCTCTGCTGGATGGGTTa                                                 <  1:4739121/76‑1 (MQ=255)
                          aTTATCGAAACAATGTCCTGGCTACATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATAGGTTAcc                                               <  1:2989207/76‑1 (MQ=255)
                          aTTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTAcc                                               <  1:3427160/76‑1 (MQ=255)
                           ttATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCa                                              <  1:336513/76‑1 (MQ=255)
                           ttATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCa                                              >  1:3632674/1‑76 (MQ=255)
                            tagcGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAg                                             <  2:4548563/73‑1 (MQ=255)
                            tATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctg                                                             >  2:5609437/1‑60 (MQ=255)
                            tATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGAGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTAcccg                                             >  2:259618/1‑74 (MQ=255)
                             aTCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGa                                            >  1:4948024/1‑76 (MQ=255)
                              tCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGAt                                           <  1:1087362/76‑1 (MQ=255)
                                gcaaaaaTGGCCTGGCAGCATACGCTGATTATGGTGGGCGGTGGCCTGTCTCTCTGATGGCTGGGTTATCAGATGc                                         >  1:4368325/3‑76 (MQ=255)
                                gAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGc                                         >  1:1189601/1‑76 (MQ=255)
                                gAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGc                                         >  2:3712155/1‑76 (MQ=255)
                                gAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGc                                         >  1:4351605/1‑76 (MQ=255)
                                gAAAAAATGGCCTGGATGCATACGCTGATTATGGCGGGCGGTGGCCTCTATCTCTGCCGGCTGGGGTGCCACATGc                                         >  1:489729/1‑76 (MQ=255)
                                 aaaaaaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCGGGATGGGTTACCAGATGCt                                        <  2:3193902/76‑1 (MQ=255)
                                  aaaaaTGGCCGGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGCTTACCTTATGCtc                                       <  1:1770644/76‑2 (MQ=255)
                                   aaaaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTAc                                      <  1:1201913/76‑1 (MQ=255)
                                   aaaaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTAc                                      <  2:285340/76‑1 (MQ=255)
                                    aaaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACg                                     >  2:780105/1‑76 (MQ=255)
                                    aaaTGGCCGGGCTGGATACGGTGATTATGGTGGGGGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACg                                     <  2:2049007/76‑1 (MQ=255)
                                     aaTGGCCTGGCTGCATCGTGTGATTGTGGTGGGCGGTGGCCTGTATCGCTGCTGGTTGGGTTCCCAGATGCTCCgt                                    <  1:4501451/76‑1 (MQ=255)
                                        ggCGTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtggt                                 <  2:2128392/76‑1 (MQ=255)
                                         gCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtggtg                                <  2:3989460/76‑1 (MQ=255)
                                          ccTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGc                               <  1:2556524/76‑1 (MQ=255)
                                          ccTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGc                               >  2:3383259/1‑76 (MQ=255)
                                          ccTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGCGGTGc                               >  1:404660/1‑76 (MQ=255)
                                           cTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCa                              >  1:3131226/1‑76 (MQ=255)
                                           cTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCa                              >  2:4260011/1‑76 (MQ=255)
                                            tGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCAc                             <  2:4163171/76‑1 (MQ=255)
                                            tGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCAc                             >  1:5425237/1‑76 (MQ=255)
                                            tGGCTGCATACGCTGATTATCGAGGGCGGTGGCCTGTATCTCTGCTGGATGGCTTACCACATGCTACGTGGTGccc                             >  2:1202142/1‑74 (MQ=255)
                                              gCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTg                           >  1:1010609/1‑76 (MQ=255)
                                               cTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGa                          <  1:2585864/76‑1 (MQ=255)
                                                tGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaa                         <  1:5664896/76‑1 (MQ=255)
                                                 gCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaaa                        <  2:1981268/76‑1 (MQ=255)
                                                 gCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaaa                        <  2:2014625/76‑1 (MQ=255)
                                                  cATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaaaa                       >  1:5670254/1‑76 (MQ=255)
                                                   aTACGCTGCTTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaaaaa                      <  1:2239235/76‑1 (MQ=255)
                                                   aTACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaaaaa                      >  2:1832027/1‑76 (MQ=255)
                                                    tACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaaaaaa                     <  1:5492563/76‑1 (MQ=255)
                                                    tACGCTGATTATGGTGGGCGGGGGCCTGTATCTCTGCTGGGTGGGGTACCAGATGCTACGTGGTGTACTGgaaaaa                     >  2:5828658/1‑76 (MQ=255)
                                                    tACGCTGATTATGGTGAGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaaaa                       <  2:1579051/74‑1 (MQ=255)
                                                    tACGCTGATTATGGTGAGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaaaa                       >  1:1579051/1‑74 (MQ=255)
                                                     tcGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAag                    <  2:3127696/75‑1 (MQ=255)
                                                     aCGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAag                    >  1:2757596/1‑76 (MQ=255)
                                                       gCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAagag                  <  1:175633/76‑1 (MQ=255)
                                                       gCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAagag                  >  1:3831386/1‑76 (MQ=255)
                                                         tGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGc                <  2:5387316/76‑1 (MQ=255)
                                                          gATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCg               <  1:1839148/76‑1 (MQ=255)
                                                          gATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGCGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCg               <  2:1560508/76‑1 (MQ=255)
                                                            ttATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGt             <  1:618036/76‑1 (MQ=255)
                                                            ttATGGTGGGCGGTGGACTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCggg             >  2:2643445/1‑75 (MQ=255)
                                                             tATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCgggt            >  2:3032979/1‑74 (MQ=255)
                                                              aTGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGttt           <  1:143505/76‑1 (MQ=255)
                                                                ggtggGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCt         <  1:606146/76‑1 (MQ=255)
                                                                 gtggGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTg        <  2:3892194/76‑1 (MQ=255)
                                                                  tggGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGGTTCttt       >  2:3650440/1‑74 (MQ=255)
                                                                   gggCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCa      <  2:2741827/76‑1 (MQ=255)
                                                                    ggCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAAAGGCGGGTTTTCCAc     >  2:2764351/1‑76 (MQ=255)
                                                                     gCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAagag                  >  2:2360758/1‑62 (MQ=255)
                                                                     gCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAagag                  <  1:2360758/62‑1 (MQ=255)
                                                                      cGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCt   <  2:3775872/76‑1 (MQ=255)
                                                                       ggTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCAcccc  >  2:5843889/1‑74 (MQ=255)
                                                                             |                                                                       
GGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTG  >  NC_000913/4005525‑4005670

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